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ABSTRACT 

Condensation  of  R-113  was  studied  using  an  evaporator/condenser  test  platform. 
The  condenser  section  contained  four  horizontal  tubes  (nominal  outer  diameter  15.9 
mm)  forming  a  vertical  in-line  column  with  a  pitch-to-diameter  ratio  of  2.25  and  a  con- 
densing length  of  1.2  m.  The  condenser  tubes  could  be  operated  either  individually  (ie. 
as  a  single  tube  apparatus)  or  as  a  small  tube  bundle.  This  allowed  investigation  of  the 
effects  of  condensate  inundation  on  different  types  of  condenser  tubes.  Tubes  tested  were 
smooth  copper  tubes,  copper/nickel  KORODENSE  roped  tubes,  KORODENSE  tubes 
wrapped  with  wire,  and  copper/nickel  finned  tubes  (26  fpi).  The  outside  heat  transfer 
coefficient,  h0,  was  calculated  by  experimentally  determining  the  overall  heat  transfer 
coefficient,  Ue ,  and  then  using  a  modified  Wilson  Plot  procedure.  Great  care  was  taken 
to  ensure  the  results  were  not  vitiated  by  the  presence  of  noncondensibles. 

Results  obtained  with  the  smooth  copper  tubes  are  in  agreement  with  published  data 
and  verify  satisfactory  operation  of  the  test  platform.  Furthermore,  problems  associated 
with  the  apparatus  encountered  by  previous  workers  were  successfully  overcome.  In 
comparison  to  the  top  smooth  tubes,  the  copper/nickel  top  KORODENSE  tube  yielded 
about  a  22%  increase  in  h0.  When  different  diameters  of  wire  were  wrapped  around  the 
copper/nickel  KORODENSE  tubes,  an  optimum  pitch-to-wire  diameter  of  7  was  found 
yielding  almost  a  90%  increase  in  h0,  compared  to  the  top  smooth  tube.  This  marked 
increase  is  presumably  due  to  surface  tension  effects  thinning  the  condensate  film.  The 
copper/nickel  finned  tubes  gave  the  best  enhancement  (approximately  a  7  fold  increase). 
Comparison  of  the  data  obtained  from  the  top  copper/nickel  finned  tubes  agreed  well 
with  the  model  of  Beatty  and  Katz. 

The  effect  of  condensate  inundation  was  to  reduce  h0  for  the  lower  tubes  compared 
to  the  top  tube  in  the  bundle.  Comparison  of  h0  for  the  second  tube  compared  to  the  top 
tube  showed  that  the  effects  of  condensate  inundation  are  reduced  most  by  wrapping  the 
KORODENSE  tubes  with  fine  wire. 
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NOMENCLATURE 

Standard  Symbols 

A  area, m2 

b  fin  spacing,  m 

CB  parameter  defined  in  Eq  (2-11) 

Cp  specific  heat,  J/kg  K 

D  diameter,  m 

e  fin  height,  m 

g  gravitational  acceleration,  mjs2 

h  heat  transfer  coefficient,  Wjm2  K 

hb  heat  transfer  coefficient  in  flooded  portion  of  finned  tube,  Wjm2  K 

hh  heat  transfer  coefficient  in  unflooded  region  of  finned  tube,  Wjm2  K 

hf  heat  transfer  coefficient  for  fins  in  unflooded  region,  Wjm2  K 

k  thermal  conductivity,  Wjm  K 

L  length,  m 

L  average  condensing  length  defined  by  Eq  (2-4) 

m  mass  flow  rate,  kg/s 

N  number  of  tubes  in  a  vertical  column 

Nu  Nusselt  number 

Pr  Prandtl  number 

q  heat  transfer  rate,  W 

q"  heat  flux,  Wjm2 

R  radius,  m 

Rf  fouling  resistance,  K/W 

Rw  wall  resistance,  K/W 

Re  Reynolds  number 

Ref  film  Reynolds  number 

T  temperature,  K 

t  fin  thickness,  m 

U  overall  heat  transfer  coefficient,  Wjm2  K 
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V  velocity,  m/s 

Greek  Symbols 

S  tape  thickness,  m 

tj  fin  efficiency 

£  aspect  ratio  of  condensate  surface,  Eq  (2-13) 

0m  rotation  angle  normal  to  fin  surface 

H  dynamic  viscosity,  kg/m  s 

p  density,  kgjmz 

a  surface  tension,  N/m 

4>f  flooding  angle  defined  in  Eq  (2-7) 

Subscripts 

c  coolant 

eq  equivalent 

ef  effective 

f  film 

fs  fin  sides 

ft  fin  tips 

i  inside  or  inlet 

o  outside  or  outlet 

r  root 

sat  saturation 
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I.     INTRODUCTION 

A.     BACKGROUND 

In  the  last  few  years,  budgetary  considerations  have  led  to  a  reduction  in  the  funding 
level  available  for  new  vessel  construction  for  the  US  Navy.  In  construction  of  ships, 
weight  can  often  be  equated  to  cost.  This  has  led  the  Navy  to  look  at  various  ways  by 
which  the  weight  of  auxiliary  systems  can  be  reduced. 

Of  particular  interest  to  the  Navy  has  been  the  issue  of  reduced  weight  and  therefore 
cost  of  refrigeration  systems.  Two  major  components  of  refrigeration  systems  are  the 
evaporator  and  condenser.  If  the  weight  of  these  two  components  can  be  reduced,  then 
the  overall  weight  and  presumably  cost,  can  also  be  reduced. 

Condenser  weight  can  be  significantly  reduced  by  the  proper  selection  of  materials 
for  condenser  construction.  A  major  weight  savings  comes  from  changing  the  material 
which  is  used  for  construction  of  the  condenser  tubes.  As  an  example,  the  DDG-51  was 
originally  intended  to  have  a  refrigeration  condenser  made  of  titanium  finned  tubes.  This 
would  provide  a  substantial  weight  savings  and  a  reduced  cost  (presumably).  However, 
the  system  was  not  installed  on  the  current  flight  of  ships.  Rather,  traditional 
copper/nickel  finned  tubes  were  utilized  for  the  refrigeration  condenser,  l 

In  re-designing  refrigeration  condensers,  changes  in  condenser  tube  material  based 
solely  on  weight  consideration  should  not  be  the  only  consideration.  The  reason  for  this 
is  that  a  lightweight  material  may  not  have  heat  transfer  characteristics  which  are  com- 
parable to  copper  or  even  copper/nickel.  Therefore,  condenser  tube  designs  have  to  be 
developed  to  maintain  or  even  enhance  heat  transfer  or  else  the  condenser  will  have  to 
be  enlarged.  This  enlargement  is  necessary  in  order  to  provide  the  requisite  surface  area 
to  maintain  the  required  level  of  heat  transfer. 

Copper/nickel  finned  condenser  tubes  have  been  used  in  refrigeration  condensers 
onboard  US  Navy  ships  for  many  years.  However,  the  use  of  titanium,  as  was  planned 
for  the  DDG-51  class  ships,  holds  the  promise  of  providing  substantial  weight  savings 
(not  to  mention  improved  corrosion  resistance  on  the  tube  side).  Unfortunately, 
titanium  has  a  thermal  conductivity  which  is  substantially  less  than  that  of  copper  or 
copper/nickel.  The  question  which  must  therefore  be  addressed  is  whether  titanium 


1  According  to  Mr.  John  Goodhue,  NAVSEA  code  SEA  56Y15,  the  reason  the  units  were 
not  installed  was  that  they  could  not  be  ready  to  support  proposed  construction  dates  lor  flight  one. 


condenser  tubes  can  be  constructed  in  such  a  manner  as  to  maintain  heat  transfer  (due 
to  enhancement)  despite  the  lower  value  of  thermal  conductivity. 

A  major  problem  along  these  lines  is  how  to  measure  the  performance  of  various 
types  of  condenser  tubes.  Typically,  this  is  done  with  single  tubes  or  small  tube  bundles. 
These  tubes  tend  to  be  relatively  short.  Often,  the  experimental  apparatus  bares  little 
resemblance  to  a  shipboard  condenser.  Thus,  any  results  become  suspect  when  they  are 
extended  to  the  real  situation. 

In  order  to  examine  heat  transfer  characteristics  of  various  designs  of  condenser 
tubes,  the  best  means  is  to  construct  a  test  platform  composed  of  an  evaporator  and 
condenser  with  a  size  comparable  to  that  of  a  shipboard  system.  Unfortunately,  this  is 
seldom  feasible  due  to  costs  involved  and  the  size  of  the  plants.  However,  it  is  possible 
to  construct  a  reasonably  large  scale  test  platform  in  which  tube  performance  can  be 
evaluated.  This  evaluation  must  include  both  single  tube  performance  as  well  as  tube 
bundle  operation.  The  reason  for  this  is  that  condensate  forming  on  one  tube  falls  onto 
lower  tubes.  This  effect,  known  as  condensate  inundation,  can  result  in  decreased  heat 
transfer  on  lower  tubes  due  to  a  thickening  of  the  condensate  layer.  Hence,  another 
aspect  which  must  be  considered  is  whether  condenser  tube  geometry  can  be  changed  in 
such  a  manner  as  to  reduce  the  adverse  effects  of  condensate  inundation.  This  could 
help  make  up  for  decreased  thermal  conductivity  associated  with  lightweight  materials 
such  as  titanium. 

The  construction  of  a  suitable  condenser/evaporator  test  platform  was  begun  at  the 
Naval  Postgraduate  School  by  Zebrowski  [Ref.  1].  The  construction  and  initial  operation 
of  the  test  platform  was  completed  by  Mabrey  [Ref.  2]  in  1988.  Mabrey  conducted  ex- 
periments using  R-113  and  R-114  on  smooth  copper  tubes  as  well  as  copper/nickel 
finned  tubes  (26  fpi).  During  the  course  of  his  experiments,  an  orange-peel  or  dimpled 
appearance  was  observed  on  the  smooth  copper  tubes.  Additionally,  the  copper/nickel 
finned  tubes  showed  a  buildup  of  "material"  between  the  fins.  This  material  was  believed 
to  be  ice  although  this  was  not  proved.  The  cause  of  the  orange-peel  appearance  was 
also  never  identified. 

B.     OBJECTIVES. 

This  thesis  concentrates  on  the  condenser  section  of  the  test  platform  with  specific 
interest  in  the  condensation  of  R-113.  R-113  was  chosen  for  this  study  because  of 
convienence  ,  ie.,  it  is  a  liquid  at  atmospheric  pressure.  The  condenser  section  consists 
of  four  horizontal  tubes  mounted  in  a  vertical  in-line  column.  The  nominal  diameter  of 


the  tubes  is  16  mm  with  a  condensing  length  of  1.2  m.  While  the  condensing  medium 
used  in  this  study  was  R-113,  the  condensation  of  many  different  types  of  fluids  can  be 
evaluated.  In  addition,  the  condenser  can  be  used  to  study  heat  transfer  on  almost  any 
type  of  condenser  tube  subject  to  the  nominal  outer  diameter  and  condensing  length  re- 
strictions indicated. 

The  main  objective  of  this  thesis  was  to  resolve  the  contamination  problem  found 
by  Mabrey  [Ref.  2],  and  to  make  the  test  platform  operational.  It  was  felt  that  such  a 
system,  once  operational,  could  be  used  to  study  heat  transfer  using  condenser  tubes 
made  of  various  materials  and  geometries.  In  addition,  since  the  tubes  are  arranged  in 
a  vertical  in-line  column,  the  effect  of  condensate  inundation  could  also  be  studied. 

C.    SPECIFIC  OBJECTIVES 

The  specific  objectives  of  this  study  were  as  follows: 

1.  Resolve  the  contamination  problem  reported  by  Mabrey  [Ref.  2] 

2.  Establish  baseline  condensation  data  by  examining  the  condensation  of  R-113  on 
smooth  copper  tubes. 

3.  Examine  the  heat  transfer  characteristics  of  R-113  on  various  tube  geometries 

4.  Investigate  the  effects  of  condensate  inundation 


II.     LITERATURE  SURVEY 

A.  INTRODUCTION 

Condensation  occurs  when  vapor  molecules  come  into  contact  with  a  surface  which 
has  a  temperature  below  the  vapor  saturation  temperature.  This  vapor,  coming  into 
contact  with  the  cooler  surface,  forms  a  liquid  which  can  accumulate  in  two  ways.  If  the 
liquid  wets  the  entire  surface,  a  continuous  film  of  liquid  is  formed.  On  the  other  hand, 
if  the  liquid  fails  to  wet  the  entire  surface,  droplets  are  formed.  These  two  types  of 
condensation  are  referred  to  as  filmwise  and  dropwise  condensation  respectively. 
Dropwise  condensation  yields  larger  heat  transfer  coefficients  but  is  difficult  to  maintain 
[Ref.  3].  All  condensers  today  are  therefore  designed  to  operate  with  film  condensation 
conditions. 

In  1916,  Nusselt  [Ref.  4]  reasoned  that  as  vapor  condenses  on  a  cooler  surface,  such 
as  a  horizontal  tube,  a  continuous  film  of  condensate  will  form.  This  film  grows  in 
thickness  as  the  condensate  flows  around  the  tube  under  the  influence  of  gravity.  Hence, 
the  thickness  of  the  condensate  is  not  uniform  around  a  horizontal  tube.  Since  the 
condensate  acts  as  a  resistance  to  heat  transfer,  Nusselt  argued  that  the  local  heat 
transfer  coefficient  will  be  a  maximum  at  the  top  of  the  tube  where  the  film  thickness  is 
least  and  decreases  around  the  tube  as  the  condensate  thickness  increases. 

B.  CONDENSATION  OF  REFRIGERANTS  ON  SMOOTH  TUBES 

Goto  et  al.  [Ref.  5]  examined  film  condensation  of  various  refrigerants,  including 
R-113,  on  smooth  copper  tubes.  Their  results  showed  values  for  the  vapor-side  heat 
transfer  coefficient  on  the  order  of  1.2  to  1.5  kWjm7  K  for  a  wall-film  temperature  dif- 
ference of  3  to  25  K.  Furthermore,  their  results  were  about  10%  below  that  predicted 
by  Nusselt  theory. 

.  Masuda  and  Rose  [Ref.  6,7]  reported  that  condensation  of  R-113  on  smooth  copper 
tubes  yielded  an  outside  heat  transfer  coefficient  of  about  1500  W\m2K  for  wall-film 
temperature  differences  of  12  to  20  K.  Their  results  showed  good  agreement  with 
Nusselt  theory. 

Marto  et  al.  [Ref.  8]  showed  that  for  condensation  of  R-113  on  smooth  copper 
tubes,  the  vapor-side  heat  transfer  coefficient  was  on  the  order  of  1000  Wjm2  K  for  a 
temperature  drop  across  the  film  of  20  to  25  K.  Their  values  were  also  in  good  agreement 
with  Nusselt  theorv. 


C.     ENHANCEMENT  OF  FILM  CONDENSATION 

The  most  commonly  used  procedure  for  enhancement  of  film  condensation  is 
through  the  use  of  integral-fin  tubes.  With  this  type  of  condenser  tube,  two  distinct  re- 
gions exist;  namely,  a  flooded  and  an  unflooded  region.  Typically,  as  fluid  condenses 
on  a  finned  surface,  it  is  pulled  from  the  fin  tip  to  the  fin  root  by  surface  tension  forces. 
Gravity  then  causes  the  fluid  between  fins  to  drain  to  the  bottom  of  the  tube  where  it 
drains  off.  Hence,  condensation  on  a  finned  tube  involves  a  combination  of  surface 
tension  and  gravity  forces.  It  is  the  added  surface  tension  which  tends  to  enhance  heat 
transfer  over  smooth  tubes  by  thinning  the  condensate  film  at  the  fin  tip. 

Since  the  present  study  uses  R-113  as  the  working  fluid,  the  remainder  of  this 
chapter  will  focus  on  previous  experiments  which  have  used  low  surface  tension  fluids 
such  as  refrigerants.  For  additional  information  involving  other  fluids,  the  reader  is  re- 
ferred to  recent  review  articles  by  Marto  [Ref.  9,10],  Webb  [Ref.  11]  and  Sukhatme  [Ref. 
12]. 

1.     Condensation  on  Single  Enhanced  Tubes:  Theoretical  Models 

The  first  theoretical  model  of  condensation  on  finned  tubes  was  published  by 
Beatty  and  Katz  [Ref.  13].  Their  theoretical  model  is  given  by: 
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Here  r\f  is  the  fin  efficiency,  Af  is  the  fin  surface  area,  Au  is  the  unfinned  surface  area,  and 
A,f  is  the  effective  total  surface  area.   The  above  equation  can  be  simplified  to: 
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where  the  equivalent  diameter  is  given  by: 
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The  equation  for  the  heat  transfer  coefficient  is  a  Nusselt-type  condensation  equation 
based  on  the  equivalent  diameter  for  a  finned  tube.  This  equation  predicted  the  exper- 
imental results  obtained  by  Beatty  and  Katz  to  within  ±  20%. 

There  are  several  important  limitations  to  the  Beatty  and  Katz  model  [Ref.  13]. 
It  does  not  distinguish  between  fin  tips  and  fin  flanks,  essentially  assuming  that  the  fins 
have  no  thickness.  Therefore  the  model  is  not  suitable  for  trapezoidal  fins.  A  second,  and 
perhaps  more  significant  limitation  is  that  the  Beatty  and  Katz  model  fails  to  take  into 
account  surface  tension  effects.  Instead,  the  model  assumes  that  the  flow  of  the 
condensate  film  is  simply  a  function  of  gravity.  While  this  may  not  pose  a  serious  limi- 
tation when  low  surface  tension  fluids  such  as  refrigerants  are  examined,  this  assumption 
can  become  critical  when  condensation  of  high  surface  tension  fluids  such  as  steam  are 
investigated. 

Another  limitation  of  the  Beatty  and  Katz  model  is  that  it  does  not  take  into 
account  condensate  flooding.  Rudy  and  Webb  [Ref.  14]  argued  that  the  flooded  portion 
of  the  tube  was  ineffective  in  condensation.  They  therefore  proposed  that  the  Beatty  and 
Katz  model  be  modified  by: 
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where  $,,  the  flooding  angle,  is  defined  as: 
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where  R0  is  the  outside  radius,  \j/  is  the  fin  tip  half  angle  and  o  is  the  fluid  surface  tension. 
It  is  interesting  to  note  that  this  modification  resulted  in  an  underprediction  of  Rudy  and 
Webb's  R-ll  condensation  data  by  about  35%. 

Webb  et  al.  [Ref.  15]  next  argued  that  surface  tension  drainage  must  be  consid- 
ered in  the  model.  They  derived  the  equation: 
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where  e  is  the  fin  height,  t  is  the  fin  thickness  and  b  is  the  width  of  the  space  between 
the  fins.  Substitution  into  the  Beatty  and  Katz  expression  yielded  agreement  with  ex- 
perimental data  of  +  10%. 

Owen  et  al.  [Ref.  16]  argued  that  the  Rudy  and  Webb  [Ref.  15]  model  had  to 
be  corrected  in  order  to  account  for  heat  transfer  in  the  flooded  part  of  the  tube.  Owen 
et  al.  proposed  the  equation: 

h0^hBK-f  +  hb(l--jf-)  (2-9) 


where  h0  is  the  average  outside  heat  transfer  coefficient  and  hb  is  the  heat  transfer  coef- 
ficient in  the  flooded  region  of  the  tube.  This  value  was  estimated  using  a  one- 
dimensional  radial  conduction  model.  This  corrected  model  of  Owen  et  al.  resulted  in 
agreement  of  +  30%  with  earlier  experimental  data. 

In  1985,  Webb  et  al.  [Ref.  17]  proposed  a  model  that  included  heat  transfer  in 
the  flooded  region  as  well  as  surface  tension  drainage  on  the  fin  flanks  and  gravity 
drainage  on  the  interfin  surface.  Their  model  equation  for  the  average  heat  transfer  co- 
efficient is: 
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In  this  expression,  hh  is  the  heat  transfer  coefficient  for  the  unflooded  root  surface  be- 
tween the  fins  and  hf  is  the  heat  transfer  coefficient  on  the  fin  flanks  in  the  unflooded 
area,  and  hb  is  the  heat  transfer  coefficient  in  the  condensate  flooded  region.  The  term 
CB  is  the  fraction  of  finned  tube  surface  flooded  by  condensate,  or, 

4>f 

where  <j)f  is  the  flooding  angle  given  in  equation  (2-7). 

In  the  Webb  et  al.  [Ref.  17]  expression  for  the  average  heat  transfer  coefficient, 
hh,  hj  and  hb  are  defined  mathematically  as  follows: 
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where  Ref  is  the  condensate  film  Reynolds  number. 
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where  Sm  is  the  length  of  the  convex  surface  and  0m  is  the  rotation  angle  normal  to  the 
fin  surface  [Ref.  18].  The  term  I  is  a  parameter  that  characterizes  the  aspect  ratio  of  the 
condensate  surface. 


Finally,  hb  is  defined  as: 


h  =  ^T-  (2-14) 


This  complex  model  of  Webb  et  al.  was  able  to  predict  heat  transfer  coefficients  to 
within  ±  20%  of  experimental  data  obtained  with  R-ll. 

Honda  and  Nozu  [Ref.  19]  and  Honda  et  al.  [Ref.  20]  have  proposed  the  most 
comprehensive  model  to  date.  This  model  treats  the  heat  transfer  problem  as  one  in- 
volving vapor-to-coolant  heat  transfer  through  a  finned  tube  wall.  In  this  model, 
trapezoidal  fins  were  divided  into  flooded  (f)  and  unflooded  (u)  regions.  Equation  (2-7) 
was  used  to  determine  the  flooding  angle.  Condensate  flow  in  between  the  fins  was  as- 
sumed to  be  driven  by  gravity  whereas  for  the  fin  surface,  condensate  flow  was  assumed 
to  be  driven  by  both  surface  tension  forces  and  gravity.  The  condensate  film  thickness 
on  the  fin  surface  was  arbitrarily  divided  into  a  thick  region,  where  no  heat  transfer  oc- 
curred, and  a  thin  film  region.  The  equation  of  the  condensate  film  thickness  was  solved 
numerically. 

Honda  and  Nozu  [Ref.  19]  showed  that  the  average  Nusselt  number  could  be 
written  in  terms  of  relevant  parameters  for  the  flooded  (f)  and  unflooded  (u)  regions  as: 
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Here,  Twu  and  Twf  are  the  dimensionless  average  wall  temperatures  at  the  fin  root  and  4>f 
is  defined  as: 
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It  should  be  pointed  out  that  Twu  and  7"^  are  determined  by  solving  a  circumferential  wall 
heat  conduction  equation.  This  wall  temperature  variation  must  be  considered  since  heat 
transfer  rates  through  flooded  and  unflooded  regions  are  very  different. 

Honda  and  Nozu  [Ref.  19]  used  their  model  to  predict  average  heat  transfer 
coefficients.  These  predictions  were  compared  to  experimental  data  on  eleven  different 
fluids  and  twenty-two  tubes.  The  predicted  values  matched  most  experimental  values  to 
within  ±  20%. 

2.     Condensation  on  Single  Enhanced  Tubes:  Experimental  Results 

The  first  experiments  conducted  on  enhanced  tubes  were  performed  by  Katz  and 
colleagues  [Ref.  21,22,  13  ].  These  authors  looked  at  condensation  of  R-12  and  R-22  on 
smooth  and  finned  tubes.  Their  results  showed  a  significant  enhancement  in  heat  transfer 
for  finned  tubes. 

In  1971,  Karkhu  and  Borovkov  [Ref.  23]  examined  condensation  of  R-113  on 
five  finned  tubes.  Their  results  showed  that  the  rate  of  condensation  for  these  tubes  with 
trapezoidal  fins  could  be  increased  by  50  to  100%  over  smooth  tubes. 

Arai  et  al.  [Ref.  24]  and  Kisaragi  et  al.  [Ref.  25]  studied  the  condensation  of  R-12 
and  R-113  respectively  on  finned  tubes  with  spiral  grooves  cut  along  the  length  of  the 
tube.  Arai  et  al.  showed  that  grooves  with  a  pitch  of  0.7  mm  gave  the  best  enhancement 
in  heat  transfer.  Kisaragi  and  co-workers  also  compared  grooved  finned  tubes  with  tubes 
on  which  the  fins  had  been  covered  with  a  porous  coating.  In  comparison  to  the  finned 
tube,  a  two  fold  enhancement  was  seen  with  the  grooved  finned  tube  while  a  four  fold 
enhancement  was  obtained  with  the  finned  porous  coated  tube. 

Carnavos  [Ref.  26]  performed  one  of  the  first  detailed  studies  of  condensation 
of  R-ll  on  enhanced  copper  tubes.  In  this  study,  eleven  different  tubes  were  tested.  The 
tubes  included  finned  tubes  as  well  as  finned  tubes  which  had  an  axial  spiral  groove  cut 
along  their  length.  Various  sizes  and  shapes  of  fins  were  included.  Enhancements  of  4 
to  6  were  seen  when  compared  to  smooth  tubes. 

Rudy  [Ref.  27]  and  Webb  et  al.  [Ref.  15]  reported  experimental  data  obtained 
for  R-ll  condensing  on  copper  tubes  with  different  fin  geometries.  Their  results  showed 
that  increasing  the  fin  density  from  748  to  1378  fpm  resulted  in  about  a  70%  increase  in 
heat  transfer  presumably  due  to  an  increased  surface  area. 

Honda  et  al.  [Ref.  28]  studied  condensation  of  R-113  on  copper  finned  tubes  of 
various  fin  density.  Their  results  showed  a  6  to  9  fold  enhancement  for  R-113 
condensation  for  finned  tubes  when  compared  to  smooth  tubes.  A  similar  study  by 
Kabov  [Ref.  29]  examined  the  condensation  of  R-12  and  R-22  on  tubes  with  different  fin 


densities  His  results  also  showed  that  the  degree  of  enhancement  is  determined  by  fin 
spacing. 

Masuda  and  Rose  [Ref.  6]  examined  heat  transfer  for  the  condensation  of  R-113 
on  low  integral-fin  copper  tubes.  Fourteen  different  fin  spacings  (0.25  to  20  mm)  were 
examined.  The  optimal  fin  spacing  was  0.5  mm  which  resulted  in  a  7.3  fold  enhancement 
in  the  vapor-side  heat  transfer  coefficient. 

Lin  and  Berghmans  [Ref.  30]  examined  condensation  of  R-ll,  R-113,  and  R-114 
on  an  enhanced  tube.  This  enhanced  tube,  known  commercially  as  "EVERFIN",  is  made 
of  trapezoidal  fins  (1378  fpm)  with  a  spacing  between  the  fins  of  0.064  mm.2  This  tube 
gave  an  8  to  10  fold  increase  in  the  vapor-side  heat  transfer  coefficient  when  compared 
to  a  smooth  tube. 

Condensation  of  R-ll  on  a  variety  of  enhanced  copper  tubes  has  been  reported 
recently  by  Sukhatme  et  al.  [Ref.  31].  All  tubes  had  trapezoidal  shaped  integral  fins.  The 
fm  density,  fin  height  and  fm-tip  half  angle  were  systematically  varied  in  order  to  deter- 
mine the  optimal  values  which  would  result  in  maximum  enhancement.  The  maximum 
enhancement  obtained  by  these  authors  was  about  10.3  relative  to  smooth  tubes.  This 
enhancement  corresponded  to  a  fin  spacing  of  0.35  mm,  a  fin  height  of  1.22  mm  and  a 
fin-tip  half  angle  of  10  degrees.  In  addition,  a  set  of  specially  enhanced  tubes  were  tested 
by  these  authors.  These  tubes  had  axial  grooves  cut  into  them  with  the  trapezoidal  fins 
present.  An  enhancement  of  12.3  was  obtained  with  these  specially  enhanced  tubes 
when  compared  to  smooth  tubes. 

Marto  et  al.  [Ref.  8]  studied  condensation  of  R-113  on  a  variety  of  copper 
integral-fin  tubes  with  rectangularly  shaped  fins.  The  optimum  fin  spacing  was  found 
to  lie  between  0.25  and  0.5  mm.  The  vapor-side  heat  transfer  coefficient  was  enhanced 
on  the  order  of  4  to  7  fold. 

Michael  et  al.  [Ref.  32]  continued  the  studies  of  Marto  et  al.  [Ref.  8]  by  looking 
at  the  effect  of  root  diameter  on  the  condensation  of  R-113  on  integral-fin  tubes.  They 
found  that  the  root  diameter  had  only  a  negligible  effect  on  heat  transfer  enhancement. 
In  addition,  these  investigators  found  that  the  optimal  fin  spacing  was  about  0.5  mm. 
They  obtained  an  enhancement  of  5  to  6. 


2  While  the  authors  quote  a  fin  spacing  of  0.064  m,  this  value  seems  too  small.  The  schematic 
diagram  of  the  fins  together  with  other  dimensions  published  by  these  authors  tend  to  suggest  that 
the  correct  fm  spacing  is  0.64  mm.  This  value  is  more  in  line  with  values  quoted  by  other  authors. 
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In  comparing  enhancements  reported  by  various  authors,  a  word  of  caution  is 
required.  In  calculating  values  for  the  vapor-side  heat  transfer  coefficients,  the  surface 
area  to  be  used  in  this  calculation  must  be  defined.  This  area  can  be  based  on  outside 
or  envelope  diameter,  root  diameter  or  total  fin  area.  The  great  majority  of  authors 
define  this  area  by  the  outside  diameter  of  the  tube  as  opposed  to  the  root  diameter.  In 
the  various  sources  cited  above,  this  has  been  the  case.  That  is,  the  surface  area  has  been 
defined  by  the  outside  tube  diameter. 

While  it  is  clear  that  integral-fin  tubes  can  result  in  an  enhancement  of  the 
vapor-side  heat  transfer  coefficient,  other  alternatives  exist.  As  an  example,  Fujii  et  al. 
[Ref.  33]  studied  the  condensation  of  R-ll  on  smooth  tubes  which  were  wrapped  with 
fine  wire.  These  authors  found  an  optimal  pitch-to-wire  diameter  ratio  of  2  which  re- 
sulted in  an  enhancement  of  2.5  in  the  vapor-side  heat  transfer  coefficient  when  com- 
pared to  copper  smooth  tubes.  These  authors  argued  that  the  enhancement  in  the  heat 
transfer  coefficient  was  due  to  a  thinning  of  the  condensate  film.  The  reason  that  the 
condensate  film  thickness  is  reduced  is  that  the  presence  of  the  wire  creates  a  low  pres- 
sure region  in  the  film  near  the  base  of  the  wires.  This  low  pressure  region  is  due  to  the 
concave  shape  of  the  condensate  film  next  to  the  wire.  As  a  result,  condensate  film  is 
pulled  along  the  tube  wall  towards  the  wire.  This  results  in  a  thinning  of  the  condensate 
film  between  the  wires. 

3.     Summary  of  Single  Enhanced  Tube  Condensation 

It  is  clear  from  the  preceding  discussion  that  theoretical  models,  despite  their 
degree  of  sophistication,  still  only  match  experimental  data  to  within  ±  20%.  This 
suggests  that  the  process  of  condensation  on  enhanced  tubes  is  extremely  complex  and 
far  from  being  fully  understood. 

In  terms  of  experimental  data  involving  condensation  of  refrigerants  on  en- 
hanced tubes,  several  comments  need  to  be  made.  The  experiments  reported  here  have 
shown  that  integral-fin  tubes  can  give  enhancements  of  the  vapor-side  heat  transfer  co- 
efficient on  the  order  of  7  to  12.  The  magnitude  of  this  enhancement  depends  upon  such 
factors  as  fm  shape,  fin  spacing,  and  fin  height.  Trapezoidal  shaped  fins  with  a  spacing 
of  0.5  mm  give  the  highest  enhancements. 

Unfortunately,  the  above  experiments  were  all  performed  on  enhanced  tubes 
made  of  copper.  In  Naval  refrigeration  systems,  where  seawater  is  used  as  the  coolant, 
corrosion  is  of  major  concern.  This  necessitates  the  need  for  copper/nickel  or  titanium 
condenser  tubes  in.  order  to  minimize  corrosion.  However,  as  can  be  seen  from  the 
foregoing   summary    of  the   literature,   essentially   no   work   has   been   done   on   the 
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condensation  of  refrigerants  on  copper/nickel  or  titanium  finned  tubes.  Since  both 
copper/nickel  and  titanium  condenser  tubes  have  a  lower  thermal  conductivity  than 
copper,  one  must  wonder  what  effect  this  will  have  on  the  vapor-side  heat  transfer  co- 
efficient.  Clearly,  experiments  using  these  tubes  must  be  carried  out. 
4.     Condensation  on  Enhanced  Tube  Bundles:  Theoretical  Studies 

Condensation  in  tube  bundles  presents  an  entirely  different  situation  than 
condensation  on  single  tubes.  In  a  tube  bundle,  two  conflicting  factors  come  into  play. 
First,  condensate  forming  on  an  upper  tube  drains  off  this  tube  onto  lower  tubes  within 
the  bundle.  This  increases  the  thickness  of  the  condensate  layer  on  lower  tubes  which 
further  increases  the  resistance  of  these  tubes  to  heat  transfer.  This  effect  is  referred  to 
as  condensate  inundation. 

Condensate  falling  from  one  tube  to  another  can  have  a  second  effect.  As 
droplets  fall  from  one  tube  to  the  next,  ripples  and  splashing  (turbulence)  are  induced  in 
the  condensate  film.  This  can  actually  enhance  heat  transfer  on  the  lower  tubes. 

Vapor  velocity  can  also  be  a  significant  factor  in  tube  bundle  condensation.  As 
the  vapor  travels  into  the  tube  bundle  and  condenses,  the  local  velocity  decreases.  This 
decrease  in  local  velocity  leads  to  a  decrease  in  the  heat  transfer  coefficient  similar  to  an 
inundation  effect  since  vapor  shear  effects  that  thin  the  condensate  film  are  reduced. 

In  analyzing  condensation  in  horizontal  tube  bundles,  Nusselt  [Ref.  4]  assumed 
that  all  condensate  flowed  from  one  tube  to  another  as  a  continuous  laminar  sheet.  He 
showed  that  the  ratio  of  the  average  heat  transfer  coefficient  for  a  vertical  column  of 
horizontal  tubes  with  respect  to  the  top  tube  could  be  expressed  as: 

-TL  =  N    4  (2-17) 

where  N  is  the  number  of  tubes  in  the  vertical  column.    A  local  coefficient  can  also  be 
calculated  for  the  Nth  tube  using  Nusselt  theory.  This  coefficient  is  given  by: 

h  —  — 

'  =A'4  -(iV-'l)4  (2-18) 


Kern  [Ref.  34]  proposed  a  less  conservative  approach.  He  argued  that  droplets 
of  condensate  falling  from  one  tube  to  another  causes  ripples  which  disturb  the 
condensate  film,  thereby  diminishing  the  effect  of  condensate  inundation.  His  math- 
ematical model  took  the  form: 
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r        — l 

rl-N  6 

-?-=N  (2-19) 

The  local  value  of  h  for  the  Nth  tube  can  be  expressed  as: 

h  -  - 

-^  =  A'6  -(.v-l)6  (2-20) 

In  1972,  Eissenberg  [Ref.  35]  proposed  that  condensate  may  not  drain  only  in 
a  vertical  direction,  but  rather,  may  be  diverted  sideways,  presumably  due  to  local  vapor 
flow  conditions.  This  so-called  "side  drainage"  model  predicts  a  less  severe  effect  of 
inundation  yields: 

h  -~ 

-p-  =  0.60  +  /V    4  (2-21) 

At  the  present  time,  numerous  investigations  have  been  carried  out  in  an  at- 
tempt to  better  understand  the  phenomenon  of  condensate  inundation.  However,  the 
data  are  scattered  and  as  a  consequence,  no  successful  theoretical  model  of  condensate 
inundation  exists  [Ref.  3].  However,  for  design  purposes,  Butterworth  [Ref.  36]  has  re- 
commended that  the  Kern  model  be  used. 

5.     Condensation  in  Enhanced  Tube  Bundles:  Experimental  Studies 

Young  and  Wohlenberg  [Ref.  37]  examined  the  condensation  of  R-12  on  a 
vertical  column  of  smooth  horizontal  tubes.  These  investigators  found  that  the  effect  of 
the  number  of  tubes  in  the  vertical  column  is  less  than  that  predicted  by  Nusselt. 

Katz  et  al.  [Ref.  38]  investigated  the  condensation  of  R-12  on  a  vertical  column 
of  six  horizontal  finned  tubes.  Their  results  also  showed  a  less  severe  inundation  effect 
than  predicted  by  the  model  of  Nusselt. 

Smirnov  and  Lukanov  [Ref.  39]  used  R-ll  to  study  condensation  on  20  hori- 
zontal rows  of  finned  tubes  arranged  in  a  triangular  fashion.  Their  results  showed  a 
dramatic  reduction  in  heat  transfer  for  the  first  five  tubes.  However,  after  this,  the  heat 
transfer  coefficient  remained  almost  constant  for  the  rest  of  the  tubes  in  the  column. 

Honda  et  al.  [Ref.  40]  have  studied  the  condensation  of  R-113  on  a  vertical 
column  of  horizontal  low  finned  tubes.  An  inundation  tube  was  used  to  simulate  addi- 
tional rows  of  tubes.  These  authors  described  four  modes  of  condensate  flow;  namely, 
droplet  mode,  column  mode,  column  and  sheet  mode,  and  sheet  mode.  Unfortunately, 
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these  authors  did  not  calculate  the  heat  transfer  coefficient  so  that  no  comparison  with 
theoretical  models  is  possible. 

Webb  and  Muravvski  [Ref.  41]  examined  the  condensation  of  R-ll  on  a  column 
of  enhanced  condenser  tubes.  The  tubes  tested  were  a  standard  1024  fpi  tube,  a 
GEWA-SC  tube,  a  Turbo-C  tube  and  a  TRED-D  tube.  The  GEWA-SC  tube  is  a  1024 
fpi  tube  but  the  fins  have  a  Y-shaped  cross-section.  The  Turbo_c  and  TRED-D  tubes 
are  similar  to  the  1024  fpi  tube  except  that  the  fms  have  a  saw-tooth  shape.  Webb  and 
Murawski  found  that  a  standard  1024  fpm  tube  showed  virtually  no  row  (inundation) 
effect.  However,  when  the  experiments  were  conducted  on  Turbo-C,  GEWA-SC  and 
Tred-D  tubes,  row  effects  (ie.,  inundation  effects)  were  seen.  The  inundation  effect  was 
greatest  for  the  Tred-D  tube  which  also  had  the  lowest  single  tube  performance.  These 
investigators  found  values  for  the  exponent  in  equations  (2-17,  2-19,  and  2-21)  on  the 
order  of  0.12  to  0.26.  Recall  that  Nusselt's  theory  used  a  value  of  0.25. 

More  recently,  Honda  et  al.  [Ref.  42]  studied  the  condensation  of  R-113  on  an 
in-line  horizontal  tube  bundle.  The  tubes  tested  consisted  of  either  flat-sided  (rectangu- 
lar) annular  fins  or  three  dimensional  trapezoidal  fins.  These  authors  found  that  the  ef- 
fects of  condensate  inundation  were  greatest  for  the  three  dimensional  fin  tubes. 
Unfortunately,  these  authors  also  had  vapor  velocities  on  the  order  of  3.4  m/s  or  greater. 
Therefore,  it  is  difficult  to  compare  their  results  directly  to  the  theoretical  models  of 
Nusselt,  Kern,  or  Eissenberg. 

An  interesting  facet  of  the  experiments  of  Honda  et  al.  [Ref.  42]  deserves  men- 
tion. These  investigators  compared  the  effects  of  condensate  inundation  on  the  vapor- 
side  heat  transfer  coefficient  for  smooth  and  finned  copper  tubes.  Their  results  showed 
that  despite  condensate  inundation,  the  heat  transfer  coefficient  for  the  finned  tubes  was 
still  substantially  above  that  of  the  unfinned  tubes.  In  addition,  the  decrease  in  the  heat 
transfer  coefficient  with  increasing  row  number  was  less  with  the  finned  tubes  suggesting 
that  the  effects  of  condensate  inundation  were  less  severe  with  finned  tubes  than  with 
unfinned  tubes.  Presumably,  this  result  is  due  to  the  fact  that  condensate  draining  off 
of  a  finned  tube  drains  in  columns  rather  than  continuous  sheets.  Lower  tubes  will 
therefore  not  be  completely  covered  with  condensate  but  instead  will  have  unflooded 
regions. 

Honda  and  Nozu  [Ref.  43]  examined  the  effects  of  condensate  inundation  on 
smooth  and  integral-fin  copper  tubes  using  a  theoretical  model  derived  recently  by 
Honda  et  al.  [Ref.  44]  This  model  was  a  single  tube  model  which  was  extended  to  include 
the  effects  of  condensate  inundation  by  realizing  that  condensate  did  not  flow  off  tubes 
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as  sheets  but  rather  as  columns.  The  model  took  account  of  the  fact  that  regions  which 
received  drainage  of  condensate  had  different  heat  transfer  coefficients  than  regions 
which  did  not.  The  results  support  the  more  recent  experimental  results  of  Honda  et  al. 
[Rcf.  42].  The  effect  of  inundation  was  less  severe  for  the  finned  tubes  than  for  the 
smooth  tubes.  The  authors  suggest  that  this  is  due  to  condensate  being  channeled  on  the 
finned  tubes  whereas  for  the  smooth  tubes,  the  entire  tube  surface  becomes  covered  with 
a  condensate  film  of  varying  thickness. 

6.     Summary  of  Enhanced  Tube  Bundle  Condensation 

Experimental  data  has  shown  that  the  Nusselt  model  for  the  effects  of 
condensate  inundation  may  be  too  conservative.  Rather,  most  experimental  data  seems 
to  fall  between  the  Kern  and  Eissenberg  predictions. 

Other  experimental  data  suggests  that  the  effects  of  condensate  inundation  are 
less  for  finned  tubes  as  opposed  to  smooth  tubes.  Presumably,  this  is  due  to  the  manner 
in  which  condensate  flows  from  one  finned  tube  to  another.  The  presence  of  the  seems 
acts  to  channel  the  flow  of  condensate  leaving  large  areas  relatively  free  of  flooding  from 
inundation. 

As  a  final  comment,  the  use  of  wire  wrapped  tubes  provides  an  interesting  al- 
ternative to  the  more  conservative  finned  tube  approach.  However,  the  effects  of 
refrigerant  condensate  inundation  on  wire  wrapped  tubes  has  yet  to  be  examined. 
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III.     EXPERIMENTAL  APPARATUS 

A.    TEST  PLATFORM 

Details  of  the  design  and  construction  of  the  condenser/ evaporator  test  platform  are 
provided  by  Zebrowski  [Ref.  1]  and  Mabrey  [Ref.  2].  Therefore,  only  a  brief  summary 
will  be  presented  here. 

The  test  platform  is  shown  schematically  in  Figure  1  while  a  photograph  of  the  unit 
is  shown  in  Figure  2.  The  unit  is  made  from  rolled  stainless  steel  plate  (6.35  mm  thick) 
and  designed  to  withstand  an  absolute  pressure  of  308  kPa  (approximately  3  bar). 
1.     Condenser  Section 

The  top  section  is  the  condenser.  It  is  1.3  m  in  length  with  an  internal  diameter 
of  0.61  m.  The  end  plates  are  made  of  stainless  steel  plate  (6.35  mm  thickness)  with  a 
diameter  of  0.71  m.  They  are  connected  to  the  condenser  shell  using  a  bolted  flange  as- 
sembly. A  rubber  gasket,  3.22  mm  thick,  provides  a  vacuum  tight  seal.  The  condenser 
unit  has  five  viewports  which  allow  the  condensation  process  to  be  viewed. 

Four  test  condenser  tubes  are  contained  in  the  condenser  forming  a  vertical  in- 
line column  with  a  pitch-to-diameter  ratio  of  2.25.  These  are  held  in  place  at  each  end 
by  a  nylon  block  bolted  to  the  end  plate.  A  rubber  gasket  provides  a  vacuum  tight  seal. 
A  stainless  steel  plate  is  Fitted  over  the  ends  of  the  condenser  tubes  and  bolted  to  the 
nylon  block.  O-rings  are  placed  in  between  the  nylon  block  and  steel  plate  in  order  to 
provide  a  vacuum  tight  seal. 

In  order  to  provide  a  means  of  pressure  control  during  testing  of  the  condenser 
tubes,  the  test  platform  was  fitted  with  five  auxiliary  condenser  tubes.  These  auxiliary 
condenser  tubes  are  constructed  of  9.53  mm  diameter  copper  tube,  wound  into  coils. 
Each  coil  is  supported  by  a  stainless  steel  rod  welded  to  one  of  the  end  plates.  The  aux- 
iliary condenser  coils  enter  and  exit  through  this  same  condenser  end  plate. 

A  stainless  steel  shroud  assembly,  shown  in  Figure  3,  is  installed  in  the 
condenser  unit.  This  shroud  surrounds  the  test  condenser  tubes  as  well  as  the  auxiliary 
condenser  coils.  A  window,  installed  in  the  side  of  the  shroud,  allows  visualization  of  the 
test  condenser  tubes  through  the  viewports.  The  purpose  of  this  shroud  is  to  funnel  va- 
por from  the  evaporator  below,  around  the  outside  of  the  shroud  and  then  vertically 
down  past  the  condenser  tubes.  The  velocity  of  the  vapor  past  the  condenser  tubes  was 
0.1  m/s. 


16 


THERMOCOUPLES    i 


THERMOCOUPLES  2 


INSTRUMENTED  1 
CONDENSER 
TUBES 


=4 


AUXILIARY 

CONDENSER 

COILS 


127  mm  VIEWPORTS 
-♦ 1.22  m 1- 


■  IEWPORTS 


0.6  m  DIA  (INSIDE) 


THERMOCOUPLE  #465         1 14  mm 


INSTRUMENTED 

CONDENSER 

TUBES 


102  mm  VIEWPORT 

mm  OUTER  SUPPORT   RING 


16mm  INSTRUMENTED  HEATER  TUBES 

»  L     ®/     203  mm  VIEWPORT 

©  /        305 -mm  OUTER  SUPPORT   RING 

16  mm  SIMULATION  HEATER  TUBES 


Figure  1.      Schematic  diagram  of  test  platform. 
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Figure  2.      Photograph  of  the  evaporator/condenser  test  platform. 
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Figure  3.      Schematic  diagram  of  condenser  showing  shroud. 
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2.  Evaporator  Section 

The  flooded  evaporator  is  formed  from  rolled  stainless  steel  and  has  an  internal 
diameter  of  0.61  m  with  a  width  of  0.279  m.  Two  viewports  allow  visualization  of  the 
boiling  process.  This  unit  is  connected  to  the  condenser  by  a  collar  assembly  0.28  m  in 
diameter  and  0.203  m  in  length.  The  volume  of  liquid  contained  in  the  evaporator  is 
approximately  0.06  m\  just  sufficient  to  cover  the  boiler  tubes. 

Vapor  is  produced  in  the  evaporator  using  a  combination  of  three  groups  of 
heaters.  The  first  group  is  the  main  bundle  and  is  composed  of  thirty-five  tubes  in  a 
triangular  pitch  (pitch  to  diameter  ratio  of  1 .2).  Fifteen  of  these  have  heaters  with  a 
nominal  rating  of  1  kW  each.  Five  of  these  heated  tubes  are  instrumented  with 
thermocouples  to  enable  bundle  boiling  data  to  be  obtained.  In  addition,  unheated 
dummy  tubes  are  used  to  guide  the  two  phase  flow  within  the  bundle. 

The  second  group  of  tubes  is  composed  of  five  simulation  heaters,  each  with  a 
nominal  rating  of  4  kW.  These  are  placed  below  the  main  bundle  in  order  to  simulate 
further  rows  of  tubes. 

The  final  group  consists  of  four  auxiliary  heaters,  each  with  a  nominal  rating 
of  4  kW.  Two  of  these  are  placed  on  each  side  of  the  evaporator  in  order  to  keep  the 
liquid  pool  at  saturation  conditions;  this  avoids  any  subcooling  of  the  liquid.  Power  to 
all  three  sets  of  heaters  is  provided  by  three  separate  STACO  240V,  23.5kVA  rheostat 
contollers. 

For  further  information  on  the  design  and  construction  of  the  test  platform,  the 
reader  is  referred  to  Zebrowski  [Ref.  1],  Mabrey  [Ref.  2],  Murphy  [Ref.  45],  and  Chilman 
[Ref.  46]. 

3.  Coolant  System 

The  coolant  is  a  54/46  by  weight  mixture  of  ethylene  glycol  and  water.  The 
coolant  system  consists  of  two  main  coolant  pumps.  These  pumps  are  0.5  hp  constant 
speed  pumps  which  take  a  suction  from  the  main  coolant  sump.  One  pump  is  used  to 
supply  five  rotameters.  A  ball  valve,  located  at  the  entrance  to  each  rotameter,  is  used 
to  control  the  flow  rate.  Flow  to  each  test  condenser  tube  is  supplied  by  an  individual 
rotameter.  The  fifth  rotameter  is  used  to  supply  one  of  the  auxiliary  condenser  coils.  A 
needle  valve  is  used  to  provide  fine  control  of  flow  through  the  coil.  This  coil  is  normally 
used  to  control  saturation  pressure  during  experimental  runs. 

The  second  coolant  pump  is  used  to  supply  coolant  flow  to  the  remaining  four 
auxiliary  condenser  coils.  Flow  through  each  coil  is  controlled  by  a  needle  valve  located 
at  the  entrance  to  each  coil.  It  should  be  noted  that  the  four  auxiliary  cooling  coils  are 


20 


seldom  needed  during  the  condenser  experiments;  their  primary  purpose  is  for  use  during 
evaporator  experiments,  (see  Chilman  [Ref.  46]). 

In  addition  to  supplying  the  four  auxiliary  condenser  coils,  the  second  coolant 
pump  also  provides  coolant  flow  to  a  refrigerant  storage  tank  and  to  a  cold  trap  located 
at  the  top  of  the  condenser.  The  storage  tank  can  be  used  to  store  refrigerant  during 
condenser  or  evaporator  tube  changeout.  The  cold  trap  is  placed  in-line  between  the 
condenser  and  a  vacuum  pump.  Vapor  from  the  condenser  is  drawn  through  the  cold 
trap  which  is  cooled  by  the  ethylene  glycol/water  mixture.  This  process  is  used  to  remove 
noncondensibles  from  the  condenser  while  preventing  loss  of  refrigerant  to  the  atmos- 
phere. 

Flow  of  the  coolant  mixture  from  the  auxiliary  condenser  cooling  coils,  storage 
tank,  and  cold  trap  are  returned  to  the  sump  tank.  Coolant  flow  from  the  test  condenser 
tubes  empties  into  a  common  header  which  then  empties  into  the  same  coolant  sump 
tank. 

In  order  to  keep  the  coolant  mixture  temperature  approximately  constant,  it  is 
cooled  by  an  eight  ton  refrigeration  system  located  outside  the  building.  This  system 
continually  recirculates  the  coolant  using  a  0.75  hp  pump.  Coolant  temperature  can  be 
maintained  at  any  temperature  from  ambient  temperature  (approximately  20  °Q  down 
to  -20  °C. 

4.     Instrumentation 

Prior  to  entering  each  condenser  tube,  the  coolant  flows  through  a  mixing 
chamber.  The  inlet  temperature  is  measured  using  a  single  copper/constantan  type-T 
thermocouple.  Coolant  flows  from  each  condenser  tube  through  an  outlet  mixing 
chamber.  Outlet  coolant  temperatures  are  then  measured  using  two  thermocouples  for 
each  tube.  An  average  outlet  temperature  is  then  taken  for  each  tube.  All  thermocouples 
had  been  previously  calibrated  by  Mabrcy  [Ref.  2]. 

Refrigerant  temperatures  were  measured  at  several  different  locations  in  the 
evaporator/condenser  test  platform  using  similar  single  thermocouples.  Vapor  temper- 
ature was  measured  in  two  separate  locations  at  the  top  of  the  condenser  and  one  lo- 
cation just  above  the  liquid  level.  Refrigerant  liquid  temperature  was  measured  at  both 
the  top  and  bottom  of  the  tube  bundle.  The  thermocouples  were  inserted  into  stainless 
steel  housings  (with  copper  tips)  which  were  installed  in  the  condenser/evaporator  test 
platform. 
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B.     CLEANING  AND  DECONTAMINATION  PROCEDURES 

As  mentioned  earlier,  there  was  a  contamination  problem  with  the  test  platform 
during  the  work  of  Mabrey  [Ref.  2].  Therefore,  the  experimental  apparatus  was  initially 
disassembled  and  cleaned  in  accordance  with  the  following  procedure. 

Prior  to  removing  the  evaporator,  all  heaters  and  thermocouples  were  disconnected. 
The  lagging  around  the  mating  flange  was  then  removed.  A  Porta-lift  was  connected  to 
the  evaporator  mating  flange  in  order  to  support  the  weight  of  the  evaporator  prior  to 
removing  the  bolts  holding  the  evaporator  and  condenser  together.  After  supporting  the 
evaporator,  the  bolts  were  removed  and  the  evaporator  lowered  from  the  condenser  unit. 

Referring  to  Figure  2,  the  condenser  unit  is  capped  at  each  end  by  a  steel  end  plate. 
The  left  end  plate  supports  the  auxiliary  condenser  units  and  is  therefore  removed  after 
the  right  end  plate.  The  test  condenser  tubes  must  be  removed  prior  to  the  disassembly 
of  the  condenser.  The  assembly  which  holds  the  condenser  tubes  in  place  is  composed 
of  a  nylon  plate  bolted  to  the  condenser  end  plate.  A  gasket  insures  a  tight  seal  between 
the  nylon  and  condenser  end  plates.  A  stainless  steel  plate  is  bolted  against  the  nylon 
plate  with  O-rings  in  between  which  ensure  an  air  tight  seal  between  the  environment 
and  the  inside  of  the  condenser. 

The  plates  holding  the  test  condenser  tubes  in  place  must  be  removed  prior  to  re- 
moval of  the  condenser  end  plates.  The  stainless  steel  plates  at  both  ends  are  removed 
by  unbolting  the  retaining  nuts.  The  nuts  holding  the  nylon  plate  in  place  are  then  re- 
moved. A  rubber  mallet  is  used  to  tap  out  the  test  condenser  tubes  from  the  left  end  of 
the  platform.  This  causes  the  right  end  nylon  plate  to  move  away  from  the  right 
condenser  end  plate.  Screwdrivers  can  then  be  used  to  pry  the  nylon  end  plate  away  from 
the  condenser. 

The  condenser  end  plates  can  be  removed  after  the  test  condenser  tubes  have  been 
removed.  As  a  word  of  caution,  these  end  plates  weigh  in  excess  of  100  pounds. 
Therefore,  they  must  be  supported  during  the  removal  process.  In  these  studies,  a 
Porta-lift  was  used  for  this  purpose.  The  right  condenser  end  plate  was  removed  first, 
followed  by  the  left  end  plate.  Prior  to  removal  of  the  left  end  plate,  the  shroud  assembly 
had  to  be  removed.  This  was  facilitated  by  first  removing  the  thermocouple  probes  which 
penetrate  the  top  of  the  condenser  unit.  The  shroud  was  then  pulled  out  of  the 
condenser. 

If  the  auxiliary  condenser  tubes  are  to  be  removed,  then  the  the  left  condenser  end 
plate  will  have  to  be  removed.  This  is  because  the  five  auxiliary  condenser  tubes  are  at- 
tached to  the  left  end  plate.   The  location  of  the  test  platform  in  relation  to  the  coolant 
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sump  is  such  that  this  end  can  not  be  removed  with  the  evaporator/condenser  test  plat- 
form in  its  present  location.  Rather,  if  the  left  head  must  be  removed,  the  entire  test 
platform  must  be  rotated  180  degrees  about  a  vertical  axis. 

Initial  inspection  of  the  condenser  unit  showed  only  one  small  area  of  discoloration 
which  surrounded  the  middle  viewport.  No  other  obvious  signs  of  contamination  were 
evident.  The  condenser  and  evaporator  were  both  cleaned  using  "3-D  Supreme  Cleaning 
Agent",  a  commercially  available  product  which  leaves  no  residue.  All  components  were 
rinsed  with  fresh  water,  dried,  then  given  a  final  cleaning  with  acetone.  The  auxiliary 
condenser  coils  were  cleaned  and  decontaminated  in  the  same  manner. 

The  copper  condenser  tubes  also  appeared  somewhat  discolored.  These  were  cleaned 
using  "Shower  Power",  a  commercially  available  cleaning  agent.  Prior  to  installation  into 
the  condenser  unit,  they  were  given  a  final  cleaning  with  acetone  to  ensure  removal  of 
any  residue  left  over  from  the  cleaning  procedure. 

During  the  reassembly,  all  gaskets  were  replaced.  Prior  to  replacement  of  the 
evaporator  unit,  a  test  blank  was  installed  on  the  condenser  and  a  vacuum  test  success- 
fully performed.  After  this  test  was  completed,  the  test  blank  was  removed  and  the 
evaporator  reinstalled.  All  penetrations  were  reconnected,  the  boiling  tube  bundle 
reinstalled  and  another  vacuum  test  performed. 

Prior  to  charging  the  system  with  R-113,  the  ethylene  glycol/water  coolant  pumps 
were  energized  and  flow  started  through  the  test  condenser  tubes.  For  this  initial  test, 
smooth  copper  tubes  were  used.  The  temperature  of  the  coolant  was  -15°C.  Since 
Mabrey  [Ref.  2]  had  noted  that  during  condensation,  a  thickening  between  the  fins  of 
copper/nickel  finned  tubes  was  occurring  which  suggested  ice  formation,  it  was  felt  that 
water  vapor  might  be  present  in  the  test  platform  since  it  had  been  open  to  the  atmos- 
phere. Therefore,  by  running  the  coolant  mixture  at  this  temperature  through  the 
condenser  tubes,  if  no  condensation  was  seen,  this  would  rule  out  the  presence  of  water 
vapor.  During  this  intitial  run,  no  condensation  was  seen  on  any  tubes.  In  addition,  the 
orange  peel  appearance  noted  by  Mabrey  on  the  smooth  copper  tubes  was  not  in  evi- 
dence. 

After  checking  the  evacuated  system  for  possible  contamination,  the  test  platform 
was  charged  with  fresh  R-113.  Coolant  flow  was  again  started.  This  time  an  orange  peel, 
or  dimpled  appearance  could  be  seen  on  the  top  condenser  tube.  This  suggested  some 
type  of  contamination  layer. 
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Power  (~  20  kW)  was  applied  to  the  heaters  while  the  coolant  temperature  was 
raised  to  ~  0°C.  The  orange  peel  appearance  on  the  top  tube  was  still  present.  Meas- 
urements of  the  vapor  temperature  at  the  top  and  bottom  of  the  condenser  showed  ap- 
proximately a  20  °C  differential,  signifying  the  presence  of  significant  amounts  of 
noncondensible  gases.  These  gases  enter  with  the  refrigerant  and  are  released  into  the 
condenser  when  the  liquid  is  boiled. 

A  vacuum  pump,  connected  to  the  top  of  the  condenser  was  used  to  remove  these 
noncondensible  gases.  When  this  was  done,  the  temperature  differential  between  the  top 
and  bottom  of  the  condenser  decreased  to  less  than  0.5°  C.  The  orange  peel  appearance 
quickly  disappeared  and  all  four  condenser  tubes  demonstrated  good  filmwise 
condensation.  Based  on  these  observations,  it  was  standard  practice  to  degas  the  system 
for  at  least  four  hours  after  the  system  was  open  to  the  atmosphere  or  fresh  R-113  was 
added  to  the  evaporator. 
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IV.     EXPERIMENTAL  PROCEDURES 

A.     GENERAL  PROCEDURES 

The  experimental  procedure  was  identical  for  all  experiments.  The  coolant  pumps 
were  energized  and  flow  started  through  one  of  the  auxiliary  condensers.  Flow  was  also 
started  through  the  four  test  condenser  tubes  at  a  nominal  coolant  velocity  of  0.2  m/s. 
The  auxiliary  heaters  were  set  at  3  kW,  the  instrumented  tube  bundle  set  at  3.5  kW,  and 
the  simulation  heaters  set  at  1.5  kW. 

The  vapor  temperature  at  the  bottom  and  the  top  of  the  condenser  was  allowed  to 
rise  to  47.5  °C  which  corresponds  to  the  saturation  temperature  of  R-113  at  atmospheric 
pressure.  The  rate  of  temperature  rise  was  controlled  to  approximately  0.3  degrees  per 
minute  by  adjusting  the  rate  of  coolant  flow  through  the  auxiliary  condenser.  During  the 
heat  up  period,  a  vacuum  pump,  connected  to  the  top  of  the  condenser  via  the  cold  trap, 
was  run  continuously  in  order  to  remove  any  noncondensiblcs.  The  absence  of 
noncondensibles  was  indicated  when  the  vapor  temperature  difference  between  the  top 
and  bottom  of  the  condenser  unit  was  less  than  0.5  °C. 

When  the  vapor  temperature  was  about  45  °C,  coolant  velocity  was  adjusted  to  be- 
tween 1  and  1.2  m/s  for  the  four  condenser  tubes.  The  heater  power  levels  were  then 
readjusted  in  order  to  maintain  a  relatively  constant  condenser  pressure  while  leaving 
one  auxiliary  condenser  coil  on-line.  At  lower  coolant  velocities,  condenser  pressure  was 
maintained  using  the  auxiliary  condenser  coils.  The  power  settings  for  the  heaters  were 
not  altered  during  the  course  of  the  experiment. 

During  all  experimental  runs,  the  vapor  temperature  at  the  top  and  bottom  of  the 
condenser  unit  was  maintained  between  46.5  and  48.5  °C  which  corresponds  to  a 
condenser  pressure  of  approximately  one  atmosphere.  Pressure  control,  was  accom- 
plished by  adjusting  coolant  flow  through  the  auxiliary  condenser  tubes. 
I.     Tube  Bundle 

Once  the  vapor  temperature  was  between  46.5  and  48.5  °C  ,  coolant  flow  was 
adjusted  to  the  desired  coolant  velocity  for  the  four  test  condenser  tubes.  Here,  the  de- 
sired coolant  velocity  is  an  input  to  the  computer  program.  The  program  then  calculates 
the  required  flowmeter  setting  using  the  calibration  curves  (Appendix  A).  The  initial 
coolant  velocity  used  was  0.2  m/s.  After  setting  the  flow  rate,  the  temperature  difference 
between  the  inlet  and  outlet  of  all  four  condenser  tubes  was  monitored  using  the  com- 
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puter  program,  DRPCON7.  Approximately  ten  minutes  was  required  for  this  temper- 
ature differential  to  stabilize.  Once  stable,  a  data  set  was  taken.  This  procedure  was 
repeated  in  increments  of  0.1  m/s  up  to  a  maximum  of  1.2  m/s.  The  procedure  was  then 
repeated  by  decreasing  the  flow  in  0.1  m/s  increments  down  to  0.2  m/s.  Two  data  sets 
were  obtained  at  each  coolant  velocity.  It  should  be  mentioned  that  for  the 
copper/nickel  finned  tubes,  the  range  of  flow  settings  was  changed.  The  minimal  flow 
setting  was  0.6  m/s  while  the  maximum  flow  setting  was  2.0  m/s.  The  reasons  for  this 
change  are  discussed  in  "3.  Enhancement  of  Outside  Heat  Transfer  Coefficient."  on 
page  30. 

2.    Single  Tube  Experiments 

After  completion  of  the  tube  bundle  runs,  the  experiment  was  repeated  for  each 
individual  tube  starting  with  the  lowest  tube  (tube  D).  Flow  through  the  tubes  above  the 
test  tube  was  stopped  and  flow  through  the  tube  of  interest  was  set  for  a  velocity  of  0.2 
m/sec.  The  vapor  pressure  was  again  maintained  using  an  auxiliary  condenser  coil. 
During  this  time,  the  upper  tubes  were  continuously  observed  for  condensate  formation. 
No  single  tube  data  was  obtained  until  condensate  formation  on  the  upper  tubes  had 
ceased.  This  typically  required  twenty  to  thirty  minutes. 

Once  the  coolant  temperature  rise  had  stabilized,  a  data  set  was  obtained.. Fol- 
lowing this,  the  coolant  flow  rate  was  increased  in  0.1  m/s  increments  up  to  the  maxi- 
mum flow  velocity.  A  data  set  was  obtained  at  each  flow  velocity.  Throughout  the 
experiment,  tubes  above  the  test  tube  were  visually  monitored  to  ensure  that  no 
condensate  formation  occurred.  The  flow  velocity  was  then  decreased  in  0.1  m/s  incre- 
ments down  to  the  minimum  coolant  flow  velocity,  with  a  data  set  obtained  at  each  flow 
velocity.  This  procedure  was  repeated  for  each  test  condenser  tube  in  the  bundle. 

B.     SPECIFIC  EXPERIMENTS 
1.     Baseline  Experiments 

The  purpose  of  the  baseline  experiments  was  to  validate  the  performance  of  the 
condenser/evaporator  test  platform  by  comparing  data  to  that  of  Mabrey  [Ref.  2].  These 
initial  baseline  experiments  were  performed  using  smooth  copper  tubes  (inside  diameter 
of  13.26  mm  ,  outside  diameter  of  15.88  mm,  thermal  conductivity  of  386  W/m  K).  The 
condensing  length  of  the  tubes  was  1.219  m.  The  tubes  were  supplied  by  Wolverine  Tube 
Co. 
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2.     Enhancement  of  Inside  Heat  Transfer  Coefficient 

The  main  purpose  of  this  study  was  to  obtain  accurate  values  of  the  outside  heat 

transfer  coefficient,  h0.  In  order  to  obtain  accurate  values  for  h0,  the  resistance  to  heat 

transfer  on  the  inside  of  the  tube  -7——  needs  to  be  minimized.   The  reason  that  this  re- 

h,A, 

sistance  must  be  minimized  is  that  if  the  inside  resistance  is  a  large  part  of  the  overall 
vapor-to-coolant  resistance,  then  h„  will  be  become  inaccurate.  In  order  to  minimize  the 
inside  resistance,  one  must  enhance  the  inside  heat  transfer  coefficient. 

The  inside  heat  transfer  coefficient  is  influenced  by  the  type  of  coolant  flow 
through  the  tubes.  In  the  test  platform,  the  coolant  is  a  54/46  percent  (by  weight) 
ethylene  glycol/water  mixture.  At  temperatures  below  0  °C,  this  fluid  is  extremely 
viscous  with  Reynolds  numbers  in  the  range  of  300  to  2800,  suggesting  laminar  flow. 
When  the  coolant  flow  is  laminar,  the  inside  heat  transfer  coefficient  is  low.  However, 
if  this  laminar  flow  can  be  mixed,  then  the  inside  resistance  can  be  decreased  and  the 
heat  transfer  coefficient  increased.  This  can  be  accomplished  through  the  use  of  radial 
mixing  elements.  Radial  mixing  elements  have  the  effect  of  mixing  the  thermal  boundary 
layer  (formed  on  the  inside  of  the  tube)  with  the  older  core  flow  thereby  improving  heat 
transfer.  It  should  be  noted  that  a  penalty  in  ing  such  elements  is  an  increased  pres- 
sure drop  which  must  always  be  kept  in  mind  in  practical  applications.  However,  for  the 
purposes  of  this  study,  increased  pressure  drop  was  not  important. 

In  this  thesis,  two  types  of  mixing  elements  were  investigated.  The  first  was  a 
twisted  tape  insert  made  of  copper  [Ref.  47,48].  This  tape  insert  had  a  thickness  of  0.559 
mm  with  a  pitch  for  a  1 80  degree  twist  of  three  times  the  inner  tube  diameter. 

The  second  type  of  insert  used  in  this  study  was  a  commercially  available  prod- 
uct referred  to  as  HEATEX  radial  mixing  elements. 3  These  elements  are  manufactured 
from  stainless  steel  and  consist  of  a  number  of  small  flexible  loops  (petals)  of  steel  wire 
attached  to  a  central  wire  core.  These  flexible  petals  are  placed  in  the  tube  such  that  they 
face  the  oncoming  flow.  Their  outside  diameter  is  slightly  larger  than  the  inside  diameter 
of  the  tube  and  are  therefore  kept  in  position  with  a  "pinch"  fit.  The  petals  disturb  the 
thermal  boundary  layer  and  promote  fluid  mixing.  Figure  4  shows  the  two  types  of  in- 
serts. It  should  be  noted  that  two  different  sizes  of  HEATEX  were  used  since  condenser 
tubes  of  two  different  inside  diameters  were  examined. 


3  manufactured  by  CAL  GAVIN  Birmingham,  England. 
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Figure  4.      Twisted  tape  insert  and  HEATEX  radial  mixing  element. 
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As  part  of  the  baseline  experiments,  the  smooth  copper  tubes  were  also  tested 
with  no  insert  installed.  Experiments  were  repeated  on  the  smooth  copper  tubes  using 
the  twisted  tape  inserts  as  well  as  the  HEATEX  elements.  This  data  was  then  used  as  a 
baseline  for  comparison  with  other  types  of  condenser  tubes.  It  should  be  mentioned 
that  when  no  insert  is  present,  the  inside  tube  resistance  to  heat  transfer  is  very  large 
relative  to  the  outside  resistance.  In  some  cases  (typically  when  the  coolant  flow  rate  is 
low),  this  can  produce  a  negative  value  for  the  outside  heat  transfer  coefficient.  This  is 
simply  the  inaccuracy  in  the  calculation  of  hc  when  the  inside  resistance  dominates  and 
clearly  demonstrates  the  need  for  an  insert  as  mentioned  above. 
3.     Enhancement  of  Outside  Heat  Transfer  Coefficient. 

The  overall  heat  transfer  coefficient,  Ue,  is  influenced  by  the  outside  as  well  as 
the  inside  heat  transfer  coefficient.  The  means  by  which  the  outside  heat  transfer  coeffi- 
cient is  enhanced  is  by  altering  the  outside  geometry  of  the  condenser  tube. 

Experiments  to  investigate  enhancement  of  the  outside  heat  transfer  coefficient 
utilized  three  types  of  condenser  tubes.  The  first  type  of  tube  tested  was  a  copper/nickel 
finned  tube  (26  fins/inch).  This  tube  has  an  inside  diameter  of  10.16  mm  with  a  thermal 
conductivity  of  43  W/m  K.  The  tubes  were  manufactured  by  Wolverine  Tube  Co.  These 
tubes  were  tested  as  a  tube  bundle  as  well  as  individually.  Experiments  were  performed 
with  no  insert  as  well  as  with  the  twisted  tape  and  HEATEX  inserts.  These  tubes  are 
shown  in  Figure  5.  One  point  must  be  mentioned  with  regard  to  the  finned  tubes.  The 
use  of  fins  on  the  outside  of  the  tube  enhances  the  heat  transfer  considerably.  Conse- 
quently, the  inside  resistance  again  becomes  a  large  fraction  of  the  overall  resistance, 
even  when  inserts  are  used. 

One  additional  point  must  be  discussed  with  regard  to  the  copper/nickel  finned 
tubes.  Because  the  size  of  the  holes  in  the  nylon  end  plates  is  fixed,  the  outside  diameter 
of  any  tubes  to  be  tested  must  be  the  same.  In  order  to  have  these  finned  tubes  the  same 
outside  diameter,  the  inside  diameter  had  to  be  reduced.  Since  the  inside  tube  diameter 
was  smaller  for  the  finned  tubes,  higher  coolant  velocities  were  obtained. 

The  second  series  of  experiments  used  a  copper/nickel  KORODENSE  tube  also 
manufactured  by  Wolverine  Tube  Co.  This  tube,  shown  in  Figure  5,  has  a  shallow  spiral 
groove  along  the  entire  condensing  length  with  a  pitch  of  10  mm.  As  with  the  finned 
copper/nickel  tubes,  experiments  were  performed  on  these  KORODENSE  tubes  with  no 
insert  as  well  as  with  the  twisted  tape  and  HEATEX  inserts.  The  tubes  were  tested  in  a 
tube  bundle  as  well  as  individually. 
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Fujii  and  co-workers  [Ref.  33]  have  previously  shown  that  the  condensation  of 
refrigerants  can  be  enhanced  by  wrapping  smooth  condenser  tubes  with  a  small  diameter 
wire.  Therefore,  a  third  series  of  experiments  were  performed  in  which  copper/nickel 
KORODENSE  tubes  were  tightly  wound  with  wire  as  shown  in  Figure  5.  Since  these 
authors  showed  that  there  was  an  optimum  wire  size  which  gave  the  greatest  enhance- 
ment of  heat  transfer  (0.5  mm  at  a  pitch  of  3  mm),  the  KORODENSE  tubes  used  in  this 
study  were  wrapped  with  one  of  three  different  sizes  of  wire.  The  sizes  used  were  0.029" 
(0.736  mm)  diameter  stainless  steel  wire,  0.049"  (1.244  mm)  diameter  stainless  steel  wire, 
and  0.0675"  (1.715  mm)  diameter  copper  wire.  The  wire  lay  in  the  shallow  grooves  of  the 
KORODENSE  tubes  and  was  soldered  to  the  tubes  at  the  ends  only.  Although  the  wire 
was  tightly  wrapped  around  the  tubes,  it  was  not  physically  attached  except  at  the  ends. 
Experiments  were  carried  out  on  each  tube  individually.  Each  tube  was  tested  in  all  four 
locations  within  the  bundle.  This  was  done  to  see  if  position  within  the  bundle  (i.e.  wake 
from  tubes  above)  had  any  effect.  Only  the  HEATEX  inserts  were  used  with  these  tubes. 

Once  an  optimal  wire  diameter  was  determined,  a  fourth  series  of  experiments 
was  conducted.  This  series  used  the  same  four  KORODENSE  tubes,  only  this  time  each 
was  wrapped  with  the  same  optimal  diameter  wire.  The  tests  were  done  with  the 
HEATEX  elements  installed.  Experiments  were  conducted  on  the  tubes  individually  in 
order  to  compare  with  the  previously  obtained  single  tube  data.  In  addition,  these  wire 
wrapped  tubes  were  tested  as  a  tube  bundle  in  order  to  investigate  whether  or  not  the 
wire  wrap  had  any  effect  on  condensate  innundation. 
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Figure  5. 


Copper/nickel  finned  tube,  Korodense  tube,  Wire  wrapped  tube. 
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V.     DATA  REDUCTION 

A.     DATA  ACQUISITION 

The  data  acquisition  system  consisted  of  a  Hewlett-Packard  series  9000  computer 
coupled  with  a  Hewlett-Packard  model  3479A  Data  Acquisition  system.  The  computer 
program  used  in  this  study  was  originally  written  by  Mabrey  [Ref.  2]  but  has  been  sub- 
stantially modified  in  the  course  of  this  thesis.  The  data  acquisition  thermocouple 
channel  assignments  are  shown  in  Table  1. 


Table  1.     THERMCOUPLE  CHANNEL  ASSIGNMENTS. 

Channel 

Temperature  Measured 

0 

Vapor  (condenser  top) 

1 

Vapor  (condenser  top) 

— 2 

Vapor  (top  of  liquid) 

3 

Liquid  (top  of  bundle) 

4 

Liquid  (top  of  bundle) 

5 

Liquid  (bottom  of  pool) 

6 

Inlet  tube  1 

7 

Inlet  Tube  2 

8 

Inlet  Tube  3 

9 

Inlet  Tube  4 

10,  11 

Outlet  Tube  1 

12,  13 

Outlet  Tube  2 

14,  15 

Outlet  Tube  3 

16,  17 

Outlet  Tube  4 

The  computer  program  allows  the  user  to  acquire  data,  reprocess  data,  or  plot  data. 
During  data  acquisition,  the  user  selects  the  fluid  velocity  desired.  Based  on  coolant  inlet 
temperature,  the  coolant  fluid  properties  and  mass  flow  rates  are  calculated  (see  Ap- 
pendix C).  The  program  then  uses  equations  derived  from  the  flowmeter  calibrations 
(Appendix  A)  in  order  to  calculate  the  appropriate  flowmeter  settings. 

Once  the  flowmeters  are  set  for  the  required  coolant  velocity,  the  computer  system 
displays  the  vapor  temperature  at  the  top  and  bottom  of  the  condenser  as  well  as  the 
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coolant  temperature  rise  for  each  condenser  tube.  When  these  temperature  rises  are 
stable  (approximately  10  minutes),  the  user  can  terminate  the  display  and  acquire  the 
data  for  that  coolant  flow  velocity.  While  the  temperature  rises  were  stabilizing,  coolant 
flow  through  an  auxiliary  condenser  tube  was  adjusted  in  order  to  maintain  constant 
pressure. 

B.     DATA  REDUCTION 

The  main  function  of  the  data  reduction  section  of  the  computer  program  is  to  back 
out  the  outside  heat  transfer  coefficient  from  an  experimentally  determined  value  of  U0 
and  an  assumed  inside  heat  transfer  coefficient  correlation.  This  is  accomplished  in  the 
following  manner. 

The  data  acquisition  system  measures  the  emf  from  the  thermocouples  listed  in  Ta- 
ble 1.  These  values  of  emf  are  converted  to  temperatures  using  correlations  obtained 
from  the  temperature  calibrations  both  contained  in  Appendix  B.  The  volume  flow  rate 
of  the  ethylene  glycol/water  mixture  is  calculated  based  on  the  desired  coolant  velocity. 
The  density,  dynamic  viscosity,  thermal  conductivity,  heat  capacity,  and  Prandtl  number 
of  the  coolant  are  calculated  using  inlet  temperatures  (see  Appendix  C  and  Appendix 
E).  From  these  quantities,  the  coolant  mass  flow  rate  and  fluid  velocity  are  calculated. 
The  coolant  Reynolds  number  for  each  condenser  tube  is  also  calculated.  The  heat 
transferred  to  the  coolant  is  then  found  from 

q  =  mCpAT  (5-1) 

where  m  is  the  mass  flow  rate  and  AT  is  the  coolant  temperature  rise.  The  outside  heat 
flux,  q"0  is  calculated  by  dividing  q  by  the  surface  area  of  the  tube  based  on  the  outside 
tube  diameter. 

The  overall  heat  transfer  coefficient,  UD,  is  given  by: 

U°  =  LMTD  {5  ~  2) 

where  LMTD  is  the  log  mean  temperature  difference  defined  as: 

AT 


LMTD 


In 


T     —  T 

1  sat  *  c, 

~T     —  T 

1  sat  *  c„ 


(5-3) 
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The  inside  heat  transfer  coefficient  (given  as  a  dimensionless  Nusselt  number,  Nuc) 
can  then  be  calculated  using  one  of  three  correlations.  The  choice  of  which  correlation 
to  use  depends  upon  the  type  of  insert  installed  inside  the  condenser  tube.  If  no  insert 
is  installed,  then  Nuc  can  be  calculated  from  a  correlation  provided  by  Hausen  [Ref. 
49]4  for  developing  laminar  flow  in  a  tube  at  constant  wall  temperature. 


Nuc  =  3.66  + 


0.0668^-i  WeiVc 


L  1  +  0.04 


-  )ReDjePre 


0.666 


(5-4) 


If  a  twisted  tape  insert  is  installed,  a  correlation  given  by  Hong  and  Bergles  [Ref. 
47]  can  be  used  to  calculate  Nuc  where 


Nue  =  5.172 


1  +  5.484jcl0"3Pr°-7 


Res 


1.25' 


0.5 


(5-5) 


Here,  y  is  the  twist  ratio  (pitch  over  diameter).  Re,  is  the  coolant  Reynolds  number, 
corrected  for  tape  thickness.  Res  is  given  by: 


Res  = 


Am 


7tfi(dj  -  AS) 


(5-6) 


where  S  is  the  tape  thickness.  This  correlation  is  valid  for  Re,  up  to  about  2300. 

A  number  of  correlations  for  determining  the  inside  heat  transfer  coefficient  have 
been  published  when  using  a  wire  matrix  insert  [Ref.  50,51,52,  53  ].  However,  these 
correlations  are  specific  for  the  type  of  mixing  element  used,  and  are  especially  depend- 
ent upon  the  loop  density  of  the  wire  matrix. 

The  HEATEX  mixing  elements  supplied  by  CAL  GAVIN  (nominal  diameters  of 
13.26  and  10.16  mm)  did  not  match  the  loop  densities  used  by  previous  investigators. 
However,  CAL  GAVIN  did  supply  experimental  data  for  the  heat  transfer  coefficient 
versus  fluid  velocity  at  three  different  temperatures  for  both  sizes  of  elements.  To  con- 
vert this  data  into  a  usable  correlation,  the  following  procedure  was  developed. 


4  Calculations  were  performed  for  coolant  flow  through  the  smooth  copper  tubes  at  coolant 
velocities  of  0.2  and  1.2  m/s.  These  calculations  confirmed  the  coolant  flow  was  not  fully  devel- 
oped. 
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The  correlations  of  Oliver  and  Aldington  [Ref.  51,  52  ]  took  the  form  of 

Nu  =  cRemPr"  (5  -  7) 


where  c,  m,  and  n  are  arbitrary  constants  which  are  dependent  upon  the  type  of  insert 
used.  It  was  assumed  that  a  similar  type  of  equation  could  be  used  in  the  present  study 
by  fitting  the  constants  to  the  data  (once  non-dimensionalized)  provided  by  CAL 
GAVIN. 

The  correlations  of  Oliver  and  Aldington  [Ref.  51,  52  ]  had  Prandtl  number  raised 
to  a  power  of  0.46.  This  was  primarily  to  take  into  account  temperature  effects.  There- 
fore, the  same  value  was  assumed  in  the  present  study.  Values  of  Nu/ZV0-46  were  plotted 
as  functions  of  Re  on  a  log-log  basis.  Figure  6  shows  the  data  obtained  for  the  13.26 
mm  elements  while  Figure  7shows  the  data  obtained  for  the  10.16  mm  elements.  Based 
on  the  appearance  of  the  non-dimensionalized  data,  a  least  squares  linear  regression  was 
performed  on  each  data  set.  The  figures  yield  the  values  of  the  constants  m  and  c  (slope 
and  intercept  respectively)  and  provide  the  following  correlations  for  the  13.26  and  10.16 
mm  elements  respectively: 


and 


Nuf  =  0.226fo°-65Pr0-46  (5-8) 


Nuc  =  0.065/?e°-76/V046  (5  -  9) 


Once  Nuc  was  obtained  from  equations  (5-4,  5-5  5-8  or  5-9),  the  inside  heat  transfer 
coefficient  was  calculated  from: 

Nu./:- 
hi  =  -fS~  (5-10) 
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HEATEX  CORRELATION  (13.26  MM) 

NU  =  0.226(RE«*.65)(PR**.46) 


to 
d 

*  "5- 

a. 


oa- 


rf 


.A 


•a 


t 1 1 — i — i — r 


.«... 


LEGEND 
o  =  T  =  0  C 
o  =  T  =  -10C 
a  =  T  =  -20  C 


t 1  i 1       I      I — I — r 


Itf 


l(f 
RE 


104 


Figure  6.      HITRAN  correlation  for  13.26  mm  elements 
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■b. 


HEATEX  CORRELATION  (10.15  MM) 

NU  =  0.067(RE**.76)(PR**.46) 
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Figure  7.      HITRAN  correlation  for  10.16  mm  elements 
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The  overall  heat  transfer  coefficient,  U0  is  found  from  the  summation  of  resistances 
to  heat  transfer  across  the  whole  tube: 

-nh-  =  -rr  +  Rf+R»  +  -rr  (5-ii) 

Mo     Mo      J  hA 

where  Rf  and  R„  are  the  fouling  and  wall  resistances  respectively.  In  order  to  find  he 
equation  (5-11)  can  be  rearranged  to  give: 

K  =  ~\ ! 1    (A    \  (5_12) 

■jj-  -R^-RfA,,-  —  \-S-J 

In  all  experiments,  the  tubes  were  meticulously  cleaned  and  therefore,  the  fouling  re- 
sistance, Rf,  was  neglected. 

In  order  to  better  assess  the  calculated  experimental  value  of  the  outside  heat 
transfer  coefficient,  some  form  of  comparison  with  prediction  must  be  made.  In  this 
case,  the  experimentally  obtained  value  of  h„  was  compared  with  Nusselt  theory.  How- 
ever, because  the  wall  temperatures  are  not  known  a  priori,  an  iteration  must  be  used 
to  determine  this  predicted  value  of  h0.  An  initial  value  of  h0  was  estimated  and  based 
on  this  estimate,  a  wall  temperature  was  calculated  using  the  equation: 

TWaIl  ~  Tsat         7  (5  —  13) 

no 

The  condensate  film  properties  are  then  calculated  using  a  reference  temperature  given 

by: 

Tfdm  =  —  Tsa[  +  —  Twa!t  (5  —  14) 

Once  these  condensate  fluid  properties  are  known  (Appendix  D),  the  theoretical  outside 
heat  transfer  coefficient  can  be  calculated  from  Nusselt  theory  [Ref  4]: 


h0  =  0.655 


F  gkfPfhfg  n 


(5-15) 


This  calculated  value  is  then  compared  to  the  initial  guess.  The  error  between  the  initial 
guess  and  the  calculated  value  are  compared.  If  the  error  is  greater  than  0.1%,  the 
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"guessed"  value  for  h0  is  updated  and  the  iteration  continued.  This  process  is  repeated 
until  the  error  between  the  calculated  and  guessed  values  is  less  than  0.1%. 

C.     MODIFIED  WILSON  PLOT 

The  outside  and  inside  heat  transfer  coefficients  could  be  determined  directly  by 
measuring  the  vapor  temperature,  the  tube  wall  temperature,  and  the  coolant  temper- 
ature rise.  Of  these,  the  most  difficult  temperature  to  measure  directly  is  the  tube  wall 
temperature.  This  would  require  the  use  of  instrumented  tubes  with  thermocouples  em- 
bedded in  the  tube  wall.  Unfortunately,  the  manufacture  of  such  instrumented  tubes  is 
both  costly  and  time  consuming,  especially  when  large  numbers  of  tubes  are  to  be  tested. 

An  alternative  to  the  use  of  instrumented  tubes  is  to  employ  the  "Modified  Wilson 
Plot"  procedure.  This  procedure  has  been  outlined  in  detail  by  Marto  [Ref.  10].  The 
modified  Wilson  plot  procedure  relies  upon  the  overall  heat  transfer  coefficient  being 
reliably  measured  from  experimental  data.  A  summation  of  heat  transfer  resistances 
from  vapor-to-coolant  is  then  expressed  as: 


11  1    /  ^o 

■jr—r-  +V.+M0  +  T  \-r 


(5-16) 


As  mentioned  above,  the  fouling  resistance  is  assumed  to  be  negligible.  This  summation 
of  resistances  can  be  rearranged  into  a  linear  equation: 


1 


RWA0 


(5-17) 


In  looking  at  the  correlations  given  above  for  hc  and  h„  it  can  be  seen  that  they  all 
employ  a  leading  coefficient.  As  an  example,  Nusselt  theory,  (equation  5-15)  can  be  ex- 
pressed in  the  form: 


h0  =  o.F 


(5-18) 


where: 


F  = 


SkjP/hfg 

v-jd0q"o 


(5-19) 


In  a  similar  manner,  the  equations  for  the  inside  heat  transfer  coefficients  (equations  5-4, 
5-5,  5-8,  and  5-9)  can  also  be  transformed  into  the  form: 
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ht  =  CO. 


(5-20) 


where  the  form  of  Q.  would  depend  on  which  type  of  tube  insert  was  being  used.  As  an 
example,  using  equation  5-5,  Q  would  become: 


10"3Prc°-7(  -y- 


1.25" 


0.5 


(5-21) 


The  linearized  version  of  the  summation  equation  (equation  5-17)  can  then  be  rear- 
ranged to  give: 


7  :     ~  Rw^o 


Ac?         1 
F  =  — - —  +  — 
CQAi       a 


This  is  a  linear  equation  with  two  unknowns,  namely,  C  and  a  .   Letting 

1 


Y  = 


—  RWA0 


and 


X  = 


A0F 

An 


we  can  obtain  a  simplified  linear  equation: 

Y=mX+ b 
where 

1 


m  = 


C 


and 


*-i 


(5-22) 


(5-23) 


(5-24) 


(5-25) 


(5  -  26) 


(5-27) 


Since  Q.  and  F  are  both  temperature  dependent,  an  iterative  solution  procedure  must  be 
used.  This  is  accomplished  by  applying  a  least  squares  fit  to  the  linear  equation. 

It  should  be  pointed  out  that  the  accuracy  of  the  modified  Wilson  plot  procedure 
relies  heavily  on  the  number  of  data  points  taken  as  well  as  the  range  of  coolant  veloci- 
ties used. 
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VI.     RESULTS  AND  DISCUSSION 

A  summary  of  the  experimental  runs  performed  in  this  thesis  is  presented  in 
Table  2.  The  abbreviation  "SM"  is  for  the  smooth  copper  tubes,  "KD"  is  for  the 
KORODENSE  tubes,  and  "CF"  is  for  the  copper/nickel  finned  tubes. 

Table  2.     SUMMARY  OF  EXPERIMENTAL  RUNS 


Data  Set 

Tube  type 

Insert 

Comments 

SM01 

Copper 

Tape 

Test  run 

SM02 

Copper 

Tape 

Test  run 

SM03 

Copper 

HEATEX 

Test  HEATEX  correlation 

SM05 

Copper 

HEATEX 

Tube  A  insert  bad 

SM06 

Copper 

HEATEX 

SM07 

Copper 

Tape 

KD01 

KORODENSE 

Tape 

1st  test  of  KORODENSE 

KD02 

KORODENSE 

None 

KD03 

KORODENSE 

HEATEX 

KD04 

KORODENSE 

Tape 

KD05 

KORODENSE 

HEATEX 

Wire  wrapped 

KD06 

KORODENSE 

HEATEX 

Wire  wrapped 

KD07 

KORODENSE 

HEATEX 

Wire  wrapped 

KD08 

KORODENSE 

HEATEX 

Wire  wrapped 

KD09 

KORODENSE 

HEATEX 

0.049"  Wire,  bad  insert 

KD10 

KORODENSE 

HEATEX 

0.049"  Wire 

KD11 

KORODENSE 

HEATEX 

0.049"  Wire 

KD12 

KORODENSE 

HEATEX 

0.049"  Wire 

CF01 

Cu/Ni  fin 

Tape 

dropwise  condensation 

CF02 

Cu/Ni  fin 

Tape 

Test  run 

CF03 

Cu/Ni  fin 

Tape 

CF04 

Cu/Ni  fin 

HEATEX 

CF05 

Cu/Ni  fin 

None 
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A.    DECONTAMINATION 

Prior  to  charging  the  system  with  R-l  13,  the  system  was  evacuated  to  approximately 
15  inHg  and  left  for  24  hours.  For  this  initial  test,  smooth  copper  tubes  were  used.  The 
ethylene  glycol  coolant  pumps  were  energized  and  flow  started  through  the  test 
condenser  tubes.  The  temperature  of  the  ethylene  glycol  was  -15°C.  No  condensation 
was  seen  on  any  tubes.  In  addition,  the  orange  peel  appearance  noted  by  Mabrey  [Ref. 
2]  was  not  observed. 

After  checking  the  evacuated  system  for  possible  contamination,  the  test  platform 
was  charged  with  fresh  R-l  13.  Coolant  flow  was  again  started.  No  orange  peel,  or 
dimpled  appearance  could  be  seen  on  the  top  condenser  tube. 

Power  (~  20  kW)  was  applied  to  the  heaters  while  the  coolant  temperature  was 
raised  to  ~  0°C.  An  orange  peel  appearance  was  noted  for  the  top  tube.  Measurements 
of  the  vapor  temperature  at  the  top  and  bottom  of  the  condenser  showed  approximately 
a  20  °C  differential. 

A  vacuum  pump,  connected  to  the  top  of  the  condenser,  was  used  to  remove  any 
noncondensible  gases.  When  this  was  done,  the  temperature  differential  decreased  to  less 
than  0.5° C.  The  orange  peel  appearance  quickly  disappeared  from  the  top  condenser 
tube.  This  suggests  that  at  least  part  of  the  contamination  problem  reported  bv  Mabrev 
was  due  to  the  presence  of  noncondensibles. 

In  a  subsequent  series  of  experiments  conducted  on  a  set  of  titanium  finned  tubes, 
a  foreign  substance  appeared  to  build  up  on  the  tubes  during  the  experimental  run.  This 
occurred  despite  continuous  operation  of  the  vacuum  pump.  The  tubes  were  removed 
from  the  condenser  and  inspected.  The  foreign  material  had  an  oily  feel  to  it  and  an  odor 
similar  to  that  of  ethylene  glycol. 

The  titanium  finned  tubes  had  smooth  extensions  (6"  in  length)  attached  to  their 
ends  using  epoxy.  These  extensions  penetrated  the  nylon  block  and  steel  plate  and  pro- 
vided the  smooth  surface  necessary  for  the  O-rings  to  give  an  airtight  seal  between  the 
condenser  and  the  atmosphere.  When  these  tubes  were  removed  from  the  condenser,  it 
was  noted  that  the  extensions  on  two  of  the  tubes  had  come  loose.  The  condenser  was 
opened  for  inspection.  Upon  removal  of  the  end  plate  and  shroud  assembly,  pools  of 
ethylene  glycol  were  seen  inside  the  condenser  shell  and  on  the  auxiliary  condenser  tubes. 
In  addition,  ethylene  glycol  was  found  around  the  shroud  assembly.  The  entire 
condenser  and  evaporator  assembly  were  therefore  dismantled  and  cleaned  in  accordance 
with  the  procedures  outlined  in  "Experimental  Apparatus". 
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After  reassembly  of  the  condenser  and  evaporator,  the  smooth  copper  tubes  were 
cleaned  and  reinstalled.  Operation  of  the  test  platform  failed  to  show  any  orange-peel 
or  dimpled  appearance  on  the  tubes.  In  addition,  all  tubes  demonstrated  good  filmwise 
condensation. 

It  should  be  mentioned  that  in  subsequent  tests,  copper/nickel  finned  tubes  were 
installed  in  the  condenser.  However,  unlike  the  titanium  finned  tubes,  the  smooth  ex- 
tensions on  these  tubes  were  soldered.  During  system  operation,  no  foreign  material 
was  noted  between  the  fins.  Condensate  retention  between  the  fins  was  negligible.  The 
tubes  demonstrated  good  film  wise  condensation. 

In  summary,  it  is  believed  that  the  contamination  problems  noted  by  Mabrev  [Ref. 
2j  were  due  to  two  problems.  The  first  problem  was  the  presence  of  noncondensibles 
during  operation  of  the  condenser.  This  resulted  in  the  relatively  low  values  for  overall 
heat  transfer  as  well  as  the  fact  that  the  top  tube,  which  would  be  most  affected  bv  the 
noncondensibles.  had  the  lowest  heat  transfer.  Secondly,  while  it  can  not  be  proved 
conclusively,  the  author  believes  that  ethylene  glvcol  may  also  have  been  present  in  the 
svstem  for  certain  tubes  and  that  this  is  what  was  the  material  noted  in  between  the  fins 
bv  Mabrev. 

B.     OVERALL  HEAT  TRANSFER  COEFFICIENT 

In  the  following  sections,  the  results  pertaining  to  the  overall  heat  transfer  coeffi- 
cient are  presented.  It  is  important  to  keep  in  mind  that  this  value  is  experimentally  de- 
termined and  does  not  rely  on  any  correlations  other  than  for  the  fluid  properties.  In 
addition,  it  should  be  mentioned  that  in  the  following  graphs,  tube  A  is  the  top  tube  in 
the  bundle,  tube  B  the  second,  tube  C  the  third,  and  tube  D  the  bottom  tube. 

1.     Smooth  Copper  Tubes 

Figure  8  shows  the  overall  heat  transfer  coefficient,  U0,  plotted  as  a  function 
of  coolant  velocity.  The  tubes  did  not  contain  any  inserts.  Note  the  low  values  for  U0. 
These  values  are  almost  constant  for  velocities  over  the  range  of  0.2  to  0.7  m/s.  In  this 
series  of  experiments,  Reynolds  numbers  were  typically  of  the  order  of  300  to  2300  in- 
dicating that  flow  is  laminar.  Without  any  insert  to  promote  mixing,  one  would  there- 
fore expect  relatively  low  values  for  Ue  .  This  is  because,  in  laminar  flow,  the  inside 
resistance  to  heat  transfer  is  relatively  large  and  therefore  leads  to  a  lower  value  for  heat 
transfer.  Finally,  there  is  essentially  no  difference  between  the  tubes  at  low  coolant  ve- 
locities. This  suggests  that  with  these  relatively  low  condensation  rates,  inundation  is  not 
a  significant  factor.  At  higher  flow  rate,  it  can  be  seen  that  the  values  for  heat  transfer 
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start  increasing  and  that  differences  between  the  tubes  become  more  apparent.  As 
coolant  velocity  is  increasing,  flow  is  beginning  to  become  turbulent.  This  reduces  the 
inside  resistance  to  heat  transfer  which  leads  to  an  increase  in  the  overall  heat  transfer. 
As  heat  transfer  increases,  the  effects  of  inundation  become  more  apparent.  This  is  pre- 
sumably due  to  the  outside  resistance  becoming  more  important  as  the  inside  resistance 
becomes  less  important. 

The  results  obtained  when  these  smooth  copper  tubes  were  tested  individually 
is  shown  in  Figure  9.  Here,  all  the  tubes  behave  similarly.  In  addition,  comparison  with 
Figure  8  shows  that  the  values  of  U„  for  the  tubes  when  tested  individually  are  approx- 
imately the  same  as  for  the  top  tube  during  the  bundle  run  as  one  would  expect. 

The  results  for  the  same  run  with  tape  inserts  installed  are  shown  in  Figure  10. 
It  can  be  seen  that  U0  increases  as  coolant  velocity  increases.  The  values  are  clearly 
greater  than  when  no  insert  is  used  (approximately  50%  greater)  indicating  the  advan- 
tage gained  when  enhancing  the  inside  heat  transfer  coefficient  in  situations  where 
coolant  flow  is  laminar.  The  effect  of  condensate  inundation  on  the  overall  heat  transfer 
coefficient  is  also  clearly  seen.  Here,  condensate  inundation  results  in  approximately  a 
15%  decrease  in  U0  when  the  top  tube  is  compared  to  the  second  tube.  Interestingly, 
there  appears  to  be  little  difference  in  the  values  for  U0  for  tubes  B,  C,  and  D  indicating 
that  the  condensate  film  is  of  similar  thickness  for  these  tubes.  The  data  for  a  similar 
run  performed  by  Mabrey  [Ref.  2]  is  shown  in  Figure  11.  The  data  obtained  in  the 
present  study  is  almost  twice  the  values  reported  by  Mabrey.  In  addition,  the  present 
data  shows  that  the  top  tube  has  the  highest  heat  transfer  which  is  to  be  expected  as- 
suming noncondensibles  are  not  present.  This  was  not  the  case  in  Mabrey's  work. 

The  results  for  these  same  tubes  when  operated  individually  are  shown  in 
Figure  12.  All  the  curves  are  essentially  the  same  with  values  for  U0  corresponding  to 
that  obtained  for  the  top  tube  during  the  bundle  run.  The  maximum  variation  is  about 
12%. 

The  next  series  of  experiments  were  performed  on  the  smooth  copper  tubes  with 
the  HEATEX  elements  installed.  Figure  13  shows  the  results  of  these  experiments  when 
the  tubes  were  operated  as  a  bundle.  As  with  the  experiments  performed  with  the 
twisted  tape  inserts,  there  is  an  increase  in  U„  as  coolant  velocity  is  increased.  In  this 
case  the  increase  is  100%  in  comparison  to  the  case  with  no  insert!  Also,  as  in  the  pre- 
vious experiment,  condensate  inundation  appears  to  result  in  about  a  15%  decrease  in 
U0  for  the  second  tube  compared  to  the  top  tube.  There  appears  to  be  little  difference 
in  the  heat  transfer  coefficients  for  tubes  B,  C,  and  D,  although  some  effect  may  be  ar- 
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Figure  8.      Bundle  run  on  smooth  copper  tubes  without  inserts. 
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Figure  9.      Single  tube  run  on  smooth  copper  tubes  without  inserts. 
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Figure  12.      Smooth  copper  tubes  with  tape  insert  operated  individually. 
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gued  at  the  higher  coolant  velocities.  This  may  reflect  the  fact  that,  at  the  higher  coolant 
velocities,  the  inside  resistance  decreases  while  the  outside  resistance  increases  in  im- 
portance. Anything  which  causes  further  increases  in  this  resistance,  such  as  condensate 
inundation,  may  cause  noticeable  effects  in  heat  transfer. 

As  before,  these  same  tubes  with  the  HEATEX  inserts  installed  were  tested  in- 
dividually; the  results  are  shown  in  Figure  14.  The  values  of  U„  for  all  four  tubes  are 
virtually  identical  and  match  the  values  for  the  top  tube  obtained  during  the  bundle  run 
(Figure  13).  The  values  appear  to  be  about  10%  greater  than  for  the  tape  insert  ex- 
periments indicating  that  the  HEATEX  provides  better  tube-side  enhancement. 

While  the  preceding  figures  have  shown  the  effect  of  inserts  and  also  inundation 
on  the  heat  transfer  coefficient,  two  important  questions  must  be  addressed.  First,  what 
is  the  uncertainty  associated  with  the  calculated  values  for  U„.  Secondly,  the  repeat- 
ability of  the  data  must  also  be  examined. 

An  uncertainty  analysis  was  performed  on  U„  as  shown  in  Appendix  F  based 
on  the  procedures  of  Kline  and  McClintock  [Ref.  54].  At  the  low  coolant  velocities,  the 
uncertainty  was  5.2%.  This  uncertainty  was  based  primarily  on  the  uncertainty  in  the 
mass  flow  measurement.  Here,  one  must  include  a  scale  interpolation  term  as  well  as  a 
term  to  account  for  the  time-wise  jitter  for  the  flowmeter.  Since  the  flow  rate  is  low  for 
this  coolant  velocity,  these  two  terms  become  more  significant  than  at  the  higher  flow 
rates.  This  is  especially  true  if  the  magnitude  of  the  time-wise  jitter  is  the  same  for  all 
flow  rates  which  was  the  case  in  this  study. 

At  high  coolant  velocities,  the  uncertainty  in  the  value  for  U0  was  about  5.6%. 
Here,  the  major  factor  driving  this  value  is  the  uncertainty  in  the  determination  of  tem- 
perature. The  reason  for  this  is  that  the  AT  between  the  tube  inlet  and  outlet  is  smaller 
than  for  low  coolant  velocities.  Since  the  relative  uncertainty  for  the  thermocouples  re- 
mains the  same  regardless  of  coolant  velocities,  uncertainty  in  the  temperature  meas- 
urement, specifically  the  LiMTD,  becomes  the  driving  factor  in  determining  the 
uncertainty  in  U0. 

The  uncertainty  analysis  performed  here  was  only  for  the  overall  heat  transfer 
coefficient.  Un.  In  this  analysis,  any  uncertainty  in  the  correlations  used  to  determine  the 
physical  properties  of  the  ethylene  glvcol  coolant  or  R-113  were  neglected.  Therefore, 
the  level  of  uncertainty  reported  here  is  slightly  conservative. 

In  order  to  examine  the  reproducibility  of  the  data,  repeat  experiments  were 
performed  on  the  smooth  copper  tubes  with  the  HEATEX  elements  installed.  It  is  im- 
portant to  mention  that  during  the  time  between  the  two  runs,  the  tubes  were  removed 
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Figure  13.      Smooth  copper  tubes  with  HEATEX  inserts.  Bundle  operation. 
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from  the  condenser,  cleaned  and  reinstalled.  The  results  for  the  duplicate  experiment  are 
shown  in  Figure  15.  When  these  values  are  compared  to  Figure  13,  differences  of  the 
order  of  3%  is  seen  in  the  values  of  U0  for  the  top  tube.  UB  for  the  lower  tubes  are  slightly 
larger  with  a  more  pronounced  difference  between  the  tubes.  However,  these  differences 
all  lie  within  the  6%  uncertainty  band  indicating  that  the  system  gives  good 
reproducibility. 

In  a  similar  manner,  the  results  for  a  repeat  experiment  on  the  smooth  copper 
tubes  tested  individually  with  the  HEATEX  installed  is  shown  in  Figure  16.  When  the 
values  for  U0  shown  here  are  compared  with  those  shown  in  Figure  14,  no  significant 
differences  are  seen.  Again,  this  demonstrates  the  good  reproducibility  obtained  by  this 
system.  However,  a  wrord  of  caution  is  in  order.  In  all  the  experiments  reported  here, 
great  care  has  been  taken  to  control  as  manv  variables  as  possible.  Of  particular  im- 
portance is  to  ensure  that  the  heater  power  levels  are  set  at  the  same  level  for  each  ex- 
periment. This  will  tend  to  keep  vapor  velocity  in  the  condenser  at  similar  levels  for  each 
experiment  which  will  result  in  better  repeatability. 

The  effects  of  having  no  insert  as  opposed  to  the  twisted  tape  or  HEATEX  in- 
sert is  shown  in  Figure  17.  In  this  figure,  the  values  for  U0  for  the  top  tube  only  have 
been  plotted  as  a  function  of  coolant  velocity  for  the  smooth  copper  tubes.  As  previously 
noted,  without  any  insert,  coolant  flow  is  laminar  and  U0  is  relatively  constant  up  to 
about  0.7  m/s.  When  either  the  twisted  tape  or  HEATEX  inserts  are  installed,  U„  in- 
creases markedly  as  coolant  velocity  increases.  Interestingly,  the  HEATEX  results  in 
about  a  10  to  14%  increase  in  U0  when  compared  to  the  twisted  tape  inserts.  Hence, 
from  a  heat  transfer  standpoint,  the  HEATEX  element  is  the  better  insert.  This  is 
probably  because  the  wire  loops  which  make  up  the  insert,  touch  much  of  the  inside 
surface  of  the  tube  thereby  "scrapping"  the  thermal  boundary  layer  from  the  wall.  With 
the  twisted  tape  insert,  mixing  is  provided  only  by  inducing  swirling  which  is  probably 
not  as  an  effective  mechanism. 

2.     Copper/Nickel  KORODENSE  Tubes 

The  results  obtained  for  the  overall  heat  transfer  coefficient,  U0,  for 
KORODENSE  tubes  operated  as  a  bundle  are  shown  in  Figure  18.  This  series  of  ex- 
periments was  performed  without  any  inserts  installed.  Note  the  similarity  between  these 
results  and  those  obtained  under  similar  conditions  with  the  smooth  copper  tubes  (Fig- 
ure S).  Again,  values  of  Ug  appear  relatively  constant  for  coolant  velocities  up  to  0.7 
m/s.  As  coolant  velocity  increases  further,  values  of  Ue  increase  dramatically.  In  addi- 
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Figure  15.      Repeat  run  on  smooth  copper  tubes  with  HEATEX  insert. 
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Figure  16.      Repeat  run  on  smooth  copper  tube  with  HEATEX.  Individual  runs. 
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tion,  an  effect  of  condensate  inundation  begins  to  be  apparent  at  these  higher  coolant 
velocities. 

Results  obtained  when  these  same  tubes  were  operated  individually  are  shown 
in  Figure  19.  Again,  when  there  is  no  inundation,  all  tubes  have  similar  behavior.  In 
addition,  the  values  obtained  here  match  those  obtained  for  the  top  tube  during  the 
bundle  experiment  (see  Figure  18).  Interestingly,  the  relative  magnitude  of  Ug  is  re- 
markably similar  to  that  obtained  for  the  smooth  copper  tubes  (Figure  9).  This  is 
somewhat  surprising.  The  KORODENSE  tubes  are  corrugated  on  both  the  inside  as  well 
as  the  outside.  Therefore,  one  would  have  expected  somewhat  higher  values  for  heat 
transfer  with  the  KORODENSE  than  with  the  smooth  copper  tube. 

Values  of  U0  for  the  KORODENSE  tubes  with  the  twisted  tape  inserts  installed 
are  shown  in  Figure  20.  This  data  was  obtained  during  bundle  operation.  As  coolant 
velocity  increased,  U„  increased  steadily.  This  was  true  for  all  four  tubes.  The  effect  of 
inundation  reduced  U„  by  about  15%  (tube  B  relative  to  tube  A).  There  appeared  to  be 
more  separation  in  the  values  of  U„  for  tubes  B,  C,  and  D  than  was  the  case  for  the 
smooth  copper  tubes  (see  Figure  10)  suggesting  a  greater  inundation  effect.  In  compar- 
ison with  the  data  obtained  for  the  smooth  copper  tubes,  the  KORODENSE  tubes  re- 
sulted in  a  slight  increase  (about  6  to  8%)  in  Uc  This  difference  may  be  due  to  the  fact 
that  the  insert  reduces  the  inside  heat  transfer  resistance  so  that  the  accuracy  of  the 
measurement  is  improved. 

The  results  obtained  for  these  same  tubes  operated  individually  are  shown  in 
Figure  21.  As  expected,  all  four  tubes  behave  similarly  when  no  inundation  is  present. 
In  addition,  the  values  for  U„  match  those  obtained  for  the  top  tube  during  the  bundle 
operation  (Figure  20).  In  comparison  to  the  smooth  copper  tube  data  (Figure  12),  the 
data  for  the  KORODENSE  tubes  appears  somewhat  higher  (about  6  to  8%).  However, 
given  the  known  uncertainty,  this  difference  is  probably  not  significant. 

Results  obtained  for  the  KORODENSE  tubes  during  bundle  operation  with  the 
HEATEX  elements  installed  are  shown  in  Figure  22.  Note  that  Ue  again  increases 
markedly  for  all  four  tube  as  coolant  velocity  increases.  An  inundation  effect  can  be  seen 
for  tubes  B,  C,  and  D  which  is  similar  to  that  seen  in  the  twisted  tape  runs.  Tube  B 
shows  approximately  a  15%  decrease  in  U0  compared  to  tube  A.  This  effect  is  magnified 
at  higher  coolant  velocities  (condensation  rates)  as  one  would  expect.  When  this  data  is 
compared  to  that  obtained  for  the  smooth  copper  tubes  with  HEATEX  installed  (see 
Figure  13),  it  can  be  seen  that  Uc  is  greater  for  the  KORODENSE  tubes.  This  suggests 
that  roping  of  the  tube  provides  some  enhancement  for  the  outside  (and  possibly  the 
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Figure  18.      KORODENSE  tubes  with  no  insert,  bundle  operation. 
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19.      KORODENSE  tubes  operated  individually,  no  inserts  installed. 
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Figure  20.      KORODENSE  tubes  with  tape  inserts,  bundle  operation. 
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Figure  21.      KORODENSE  tubes  with  tape  inserts,  individual  operation. 
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inside)  heat  transfer  coefficients.  While  we  expect  some  small  difference,  it  must  be 
pointed  out  that  the  differences  observed  here  are  at  the  bounds  of  the  calculated  un- 
certainty. 

Results  for  these  tubes  tested  individually  are  shown  in  Figure  23.  As  before, 
the  values  of  Ue  are  similar  for  all  four  tubes  and  match  the  values  obtained  for  the  top 
tube  during  the  bundle  experiments  (see  Figure  22).  In  comparison  to  the  data  obtained 
for  the  smooth  copper  tubes,  the  values  of  Ue  are  somewhat  higher  as  expected,  but  are 
still  at  the  bounds  of  the  calculated  uncertainty. 

The  effect  of  using  no  insert,  twisted  tape  and  HEATEX  on  Ue  can  be  seen  in 
Figure  24.  The  data  are  for  the  top  tube  only.  Note  that,  as  with  the  smooth  copper 
tubes  (Figure  17),  the  HEATEX  results  in  about  a  10  to  14%  increase  in  U0  in  com- 
parison to  the  twisted  tape  insert.  Both  provide  significant  enhancement  in  U0  when 
compared  to  no  insert  (almost  100%  enhancement).  As  with  the  smooth  copper  tubes, 
when  no  insert  is  used,  Ue  appears  relatively  constant  for  coolant  velocities  over  the 
range  of  0.2  to  0.7  m/s  indicating  laminar  flow. 
3.     Wire  Wrapped  KORODENSE  Tubes 

This  series  of  experiments  was  performed  to  determine  an  optimum  wire  diam- 
eter as  discussed  in  "Experimental  Procedure".  All  tubes  were  tested  individually  and  all 
had  HEATEX  mixing  elements  installed.  Figure  25  shows  U0  as  a  function  of  coolant 
velocity  for  the  case  of  an  ordinary  KORODENSE  tube  together  with  tubes  wrapped 
with  0.029",  0.049"  and  0.0675"  diameter  wire,  each  tested  in  the  top  position  (tube  A). 
The  0.049"  diameter  wire  gives  the  biggest  enhancement  in  Ue  (about  25%).  While  the 
other  two  wire  diameters  also  provide  enhancement  in  U0  (about  15%),  it  is  below  that 
of  the  0.049"  wire.  The  no  wire  case  agrees  with  the  previously  obtained  data  (see 
Figure  22). 

Figure  26  shows  similar  data  but  in  this  case  for  tube  B.  Again,  the  0.049"  di- 
ameter wire  provides  the  greatest  enhancement  in  U0  (again  about  25%). 

Figure  27  shows  what  happens  when  the  different  diameter  wires  are  used  on 
tube  C.  In  this  case  an  interesting  fact  can  be  observed.  Here,  both  the  0.049"  and  0.029" 
wire  provide  virtually  the  same  degree  of  enhancement  for  U0  (about  20%)  with  both 
being  better  than  the  0.0675"  wire.  Figure  28  shows  the  results  obtained  for  tube  D. 
Note  that  while  the  0.049"  wire  again  provides  the  greatest  degree  of  enhancement  for 
U„  (about  16%),  the  differences  between  the  three  different  diameter  wires  are  not  as 
great. 
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Figure  22.      KORODENSE  tubes  with  HEATEX  inserts,  bundle  operation. 
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Figure  23.      KORODENSE  tubes  with  HEATEX  inserts,  individual  operation. 
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Figure  25.      Effect  of  wire  diameter  on  U0  for  tube  A. 
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Figure  26.      Effect  of  wire  diameter  on  U0  for  tube  B. 
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Figure  28.      Effect  of  wire  diameter  on  U9  for  tube  D. 
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Figure  27  and  Figure  28  suggest  that  the  position  of  the  wire  wrapped  tube 
within  the  bundle  may  affect  the  overall  heat  transfer  coefficient.  In  order  to  gain  further 
insight  on  this,  the  data  obtained  from  these  experiments  was  replotted  so  that  Ue  for  a 
given  wire  diameter  was  plotted  as  a  function  of  coolant  velocity  for  all  four  tube  posi- 
tions. The  data  for  the  tube  without  any  wire  is  shown  in  Figure  29.  The  data  is  all  very 
similar  and  no  distinct  trends  are  evident.  Note  that  comparison  of  this  data  with 
Figure  23  again  demonstrates  the  reproducibility  of  the  system. 

Figure  30  (0.029"  wire  diameter)  shows  that  tubes  A,  B  and  C  behave  similarly. 
Tube  D  produces  consistently  lower  values  for  U0  than  the  other  tubes.  Figure  31 
(0.049"  wire)  demonstrates  what  appears  to  be  a  positional  effect  on  U0.  The  data  for 
tube  A  and  B  is  virtually  identical  and  lies  above  that  of  tubes  C  and  D  which  are  also 
virtually  identical.  Figure  32  suggests  that  there  is  also  a  positional  effect  for  the 
0.0675"  wire.  The  top  tube  seems  to  provide  the  greatest  enhancement  in  U0.  The  level 
of  enhancement  for  tubes  C  and  D  are  approximately  the  same,  both  being  lower  than 
tubes  A  and  B. 

Figure  33  shows  U0  plotted  as  a  function  of  wire  diameter  for  an  arbitrarily 
chosen  coolant  velocity  of  1.0  m/s.  Note  that  the  strongest  positional  affect  appears  to 
be  for  the  0.049"  wire.  Also,  it  is  clearly  evident  that  the  0.049"  wire  provides  the  greatest 
enhancement  for  heat  transfer. 

The  possible  positional  effect  may  be  due  to  vapor  velocity.  Although  the  ve- 
locity of  the  vapor  was  very  small  in  this  study  (approximately  0. 1  m/s)  there  may  still 
be  an  effect.  As  the  vapor  penetrates  the  tube  bundle,  it  slows  down.  As  its  velocity  is 
reduced,  it  can  no  longer  help  to  strip  condensate  from  the  tubes.  Hence,  heat  transfer 
is  decreased.  This  effect  looks  much  like  an  inundation  effect.  An  interesting  point  is 
why  this  shows  up  in  these  experiments  and  not  with  the  plain  KORODENSE  or 
smooth  copper  tubes.  With  the  wire  wrapped  tubes,  the  condensate  layer  is  presumably 
thin.  Therefore,  anything  which  could  help  to  thin  it  further  would  be  of  importance  in 
improving  heat  transfer.  The  loss  of  this  small  amount  of  vapor  shearing  may  be  im- 
portant in  this  case. 

In  this  study,  the  0.049"  diameter  wire  appears  to  give  the  best  enhancement  in 
heat  transfer  (see  Figure  33).  This  yields  a  pitch  to  wire  diameter  ratio  of  7.2.  Fujii  et 
al.  [Ref.  33]  found  that  for  condensation  of  R-12  on  wire  wrapped  smooth  copper  tubes, 
the  optimal  pitch-to-wire  diameter  ratio  was  on  the  order  of  2.  This  difference  could  be 
due  to  our  use  of  the  KORODENSE  tubes.  For  a  pitch-to-wire  diameter  ratio  of  2,  we 
would  have  had  to  use  wire  which  was  of  the  order  of  5  mm  in  diameter.  This  would  have 
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Figure  30.      Effect  of  0.029"  wire  on  U0  for  all  tubes. 
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Figure  31.      Effect  of  0.049"  wire  on  U0  for  all  tubes. 


73 


UZERO  WIRE  WRAPPED  KORODENSE 

0.0675"  WIRE,  R-113,  HEATEX 


M  - 

CM 

|o 
>§■ 

O 

OS 

w 

©• 

03 


o 
o" 


o 
o 


XT 

o 


+ 


—I — 

0.3 


LEGEND 
a  =  TUBE  A 
o  =  TUBE  B 
a  =  TUBE  C 
+  =  TUBE  D 


— i — 

0.4 


—I — 
0.6 


-1 — 
0.7 


—I — 
0.6 


—1 — 
0.9 


-I — 
1.1 


0.1        0.2 


0.5 


VELOCITY  (M/SEC) 


Figure  32.      Effect  of  0.0675"  wire  on  U0  for  all  tubes. 


1.2 


l.s 


74 


WIRE  WRAPPED  KORODENSE  TUBES 

COOLANT  VELOCITY  =  1  M/S 


o 
©■ 

CM 


O 

K 
W 


LEGEND 
□  =  TUBE  A 
o  =  TUBE  B 
a  =  TUBE  C 
+  =  TUBE  D 


B 


0.0000 


0.0225  0.0460 

WIRE  DIAMETER  (IN) 


0.0670 


Figure  33.      Um  as  a  function  of  wire  diameter  for  the  four  tubes. 


75 


been  unrealistic  for  our  experiments.  It  may  be  that  the  optimal  pitch-to- wire  diameter 
may  in  fact  be  about  2.  Additional  experiments  need  to  be  performed  with  different 
pitch  KORODENSE  tubes. 

4.  KORODENSE  Tubes  Wrapped  with  0.049"  Wire 

Once  the  optimal  wire  diameter  had  been  determined,  two  identical  experiments 
were  performed  on  the  KORODENSE  copper/nickel  tubes  wrapped  with  the  0.049"  di- 
ameter wire.  The  reason  that  two  experiments  were  performed  was  to  check  the  repeat- 
ability of  the  data.  Figure  34  and  Figure  35  show  the  results  from  these  two  experiments 
for  bundle  operation  of  the  tubes.  The  two  experiments  show  excellent  repeatability.  In 
comparing  these  results  with  the  non-wire  wrapped  tubes  (Figure  22),  several  interesting 
points  can  be  made.  The  values  of  U„  for  all  four  tubes  are  greater  for  the  wire  wrapped 
tubes  than  for  the  plain  tubes  (approximately  16%  for  the  top  tube).  In  addition,  the 
inundation  effect  (seen  as  a  lowering  of  Ug)  is  reduced  for  the  wire  wrapped  tubes. 
Finally,  note  that  the  differences  between  the  tubes  are  more  apparent  than  for  the 
non-wire  wrapped  tubes.  This  suggests  that  the  wire  wrapping  may  have  an  effect  upon 
condensate  inundation.  This  will  be  addressed  in  more  detail  later. 

Figure  36  and  Figure  37  show  the  results  obtained  when  the  tubes  wrapped 
with  the  0.049"  diameter  wire  were  operated  individually.  All  the  tubes  behave  similarly 
and  demonstrate  good  reproducibility  between  the  two  runs. 

5.  Copper/Nickel  Finned  Tubes 

Figure  38  shows  the  results  obtained  when  no  insert  was  installed  and  the  tubes 
were  operated  as  a  bundle.  Comparison  of  these  results  with  those  shown  in  Figure  8 
for  the  smooth  tubes  shows  that  the  finned  tubes  yield  higher  heat  transfer  coefficients 
as  expected  (approximately  12%).  Also,  as  with  the  smooth  copper  tubes,  the  top  tube 
again  displays  the  highest  value  of  U0,  due  to  the  fact  that  the  lower  tubes  suffer  from 
condensate  inundation.  The  single  tube  data  for  the  finned  tubes  is  shown  in 
Figure  39.  All  four  tubes  are  virtually  the  same  as  the  top  tube  during  bundle  operation. 

Figure  40  and  Figure  41  show  the  results  obtained  for  the  finned  tubes  with  the 
tape  insert  installed  during  bundle  and  individual  tube  operation  respectively.  In  bundle 
operation,  the  values  of  U0  are  nearly  double  those  obtained  for  both  the  smooth  copper 
tubes  (Figure  10)  and  the  KORODENSE  tubes  (Figure  20)  using  similar  inserts.  As 
with  the  smooth  and  KORODENSE  tubes,  inundation  again  causes  a  reduction  in  Ue . 
This  reduction  for  the  second  tube  appears  to  be  about  10%.  Comparison  of  the  single 
tube  operation  (Figure  41)  with  the  bundle  operation  (Figure  40)  shows  that  all  four 
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Figure  34.      0.049"  uire  wrapped  KORODENSE  tubes,  bundle  operation,  expt.  1. 
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Figure  35.      0.049"  wire  wrapped  KORODENSE  tubes,  bundle  operation,  expt.  2. 
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Figure  36.      0.049"  wire  wrapped  KORODENSE  tubes,  individual  operation. 
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Figure  37. 


0.049"  wire  wrapped  KORODENSE  tubes,  individual  operation. 
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Figure  38.      Copper/nickel  finned  tubes,  no  insert,  bundle  operation. 
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Figure  39.      Copper/nickel  finned  tubes,  no  insert,  individual  operation. 
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tubes  operated  individually  have  approximately  the  same  value  of  U0  as  the  top  tube 
during  the  bundle  operation. 

Data  obtained  from  experiments  on  the  copper/nickel  finned  tubes  with  the 
HEATEX  elements  installed  are  shown  in  Figure  42  and  Figure  43  for  bundle  and  in- 
dividual tube  operation  respectively.  Note  that  in  comparison  with  the  tape  insert  ex- 
periments (Figure  40  and  Figure  41),  use  of  the  HEATEX  inserts  results  in  about  a 
20%  increase  in  U„  which  again  indicates  the  improvement  afforded  by  using  this  type 
of  insert.  Again  the  inundation  effect  results  in  a  decrement  of  10  to  14%.  During  in- 
dividual tube  operation,  all  four  tubes  again  behave  like  the  top  tube  during  the  bundle 
operation. 

As  with  the  smooth  copper  tubes  and  the  KORODENSE  tubes,  the  effects  of 
using  different  types  of  inserts  can  be  summarized.  This  has  been  done  in  Figure  44. 
This  data  is  for  the  top  tube  only.  The  similarity  with  data  from  the  other  tubes  (see 
Figure  17  and  Figure  24)  is  good  with  HEATEX  yielding  approximately  a  20%  increase 
in  Ug  over  the  twisted  tape  insert  and  a  remarkable  200%  increase  over  the  case  with  no 
insert! 

C.     OUTSIDE  HEAT  TRANSFER  COEFFICIENT 

The  following  figures  deal  with  the  outside  heat  transfer  coefficient,  h„.  In  all  cases, 
except  as  specifically  noted,  the  outside  heat  transfer  coefficient  was  determined  using 
the  modified  Wilson  Plot  procedure.  As  before,  tube  A  is  the  top  tube  with  tube  D  being 
the  bottom  tube  in  the  bundle.  In  addition,  no  values  were  obtained  for  h0  for  exper- 
iments which  had  no  insert.  The  reason  for  this  is  that,  for  the  coolant  velocities  used, 
without  inserts,  the  inside  resistance  was  simply  too  large  for  any  accurate  results  to  be 
obtained  for  h0. 

1.     Smooth  Copper  Tubes 

Figure  45  and  Figure  46  show  the  values  obtained  for  h„  for  the  smooth  copper 
tubes  with  the  tape  insert  installed  during  bundle  and  individual  tube  operation  respec- 
tively. Several  features  deserve  comment.  First,  values  of  h0  for  the  top  tube  are  below 
those  predicted  by  Nusselt  theory.  These  differences  range  from  about  15%  at  low 
coolant  velocities  to  about  7%  at  the  higher  coolant  velocities.  In  addition,  note  the 
inundation  effect  (seen  as  a  reduction  in  h0  for  tubes  B,  C,  and  D).  Reductions  in  h0  range 
from  28%  to  19%  as  coolant  velocity  increases  from  0.2  to  1.2  m/s  respectively.  As  was 
the  case  for  the  overall  heat  transfer  coefficient,  when  the  tubes  were  operated  individ- 
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Figure  40.      Copper/nickel  finned  tubes,  tape  insert,  bundle  operation. 
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Figure  41.      Copper/nickel  finned  tubes,  tape  insert,  individual  operation. 
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Figure  42.      Copper/nickel  finned  tubes,  HEATEX  insert,  bundle  operation. 
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Figure  43.      Copper/nickel  finned  tubes,  HEATEX  insert,  individual  operation. 
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Figure  44.      Summary  of  effect  of  tube  inserts  on  Ue 


ually,  all  had  values  for  h0  which  were  virtually  identical  to  those  obtained  for  the  top 
tube  during  bundle  operation. 

Figure  47  and  Figure  48  show  values  for  h0  obtained  when  the  HEATEX  ele- 
ments were  installed  in  the  copper  tubes  during  bundle  and  individual  tube  operation 
respectively.  As  was  the  case  for  the  tape  insert,  values  of  h0  range  from  about  10  to 
14%  below  that  predicted  from  Nusselt  theory.  Again,  there  is  a  substantial  inundation 
effect  which  results  in  a  15  to  18%  decrease  in  h0  when  comparing  the  top  tube  to  the 
second  tube.  Note  that  the  values  for  h0  obtained  for  the  individual  tubes  agree  with 
those  obtained  for  the  top  tube  during  bundle  operation.  Finally,  comparison  of  the 
value  for  hc  obtained  when  the  tape  inserts  were  used  as  opposed  to  the  HEATEX  ele- 
ments reveal  no  effect  on  h0  which  is  to  be  expected  since  we  are  examining  the  outside 
heat  transfer  coefficient. 

It  is  interesting  to  note  that  Goto  et  al.  [Ref.  5]  found  that  values  of  h0  for 
R-113  condensing  on  long,  smooth  copper  tubes  were  about  10%  below  that  predicted 
by  Nusselt  theory.  Marto  et  al.  [Ref.  8  ]  found  that  the  vapor-side  heat  transfer  coeffi- 
cient for  R-113  condensing  on  short  (0.133m)  smooth  copper  tubes  was  in  relatively 
good  agreement  with  Nusselt  theory.  This  result  was  subsequently  confirmed  by  Michael 
et  al.  [Ref.  32].  However,  it  is  important  to  point  out  that  both  the  latter  reports  used 
very  short  condenser  tubes.  The  tubes  used  by  Goto  et  al.  were  about  0.6  m  in  length, 
while  the  tubes  used  in  this  study  had  a  condensing  length  of  1.2  m.  Hence  tube  length 
may  be  an  important  consideration.  The  reason  for  this  is  as  follows.  The  coolant  en- 
ters the  condenser  tube  with  a  given  inlet  temperature.  As  the  coolant  flows  through  the 
tube,  the  coolant  temperature  increases.  The  wall  termperature  is  also  changing  down 
the  length  of  the  tube.  At  the  tube  entrance,  the  temperature  difference  between  the  tube 
wal  and  the  vapor  saturation  temperature  causes  the  formation  of  a  thick  condensate 
layer.  This  layer  is  thinner  at  the  tube  outlet  due  to  the  smaller  wall-vapor  saturation 
temperature  difference.  This  causes  condensate  to  flow  in  the  axial  direction.  These  fac- 
tors taken  together  imply  that  the  temperature  difference  is  a  function  of  the  length  of 
the  tube.  This  factor  is  not  taken  into  account  in  the  Nusselt  analysis.  As  coolant  ve- 
locity increase^  the  difference  betwee  the  coolant  inlet  and  outlet  temperatures  is  re- 
duced. The  wall  temperature  becomes  more  constant  along  the  tube  length  so  that  values 
for  h0  are  in  better  agreement  with  those  determined  experimentally. 

The  effect  of  using  the  twisted  tape  and  HEATEX  insert  on  the  inside  and  out- 
side heat  transfer  coefficients  is  demonstrated  in  Figure  49  and  Figure  50  respectively. 
In  both  cases,  the  values  for  h0  are  virtually  identical  as  expected.  However,  the 
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Figure  45.      Smooth  Copper  tubes,  tape  insert,  bundle  operation. 
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Figure  46.      Smooth  Copper  tubes,  tape  insert,  individual  operation. 
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Figure  47.      Smooth  Copper  tubes,  HEATEX  insert,  bundle  operation. 
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Figure  48.      Smooth  Copper  tubes,  HEATEX  insert,  individual  operation. 
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HEATEX  elements  essentially  double  the  inside  heat  transfer  coefficient.  This  is  what 
provides  the  improvement  in  overall  heat  transfer  coefficient  when  the  HEATEX  ele- 
ments are  used. 

The  effects  of  inundation  were  examined  by  calculating  an  average  and  a  local 
heat  transfer  coefficient.  The  local  h0  is  found  by  dividing  the  value  for  h0  for  a  given 
tube  by  the  value  for  the  top  tube.  Values  obtained  are  then  compared  with  Nusselt 
theory  (Eq.  2-18)  and  the  Kern  model  (Eq.  2-20).  The  average  heat  transfer  coefficient 
is  determined  by  taking  an  average  value  of  h0  for  the  tube  of  interest  plus  the  tubes 
above  it  and  dividing  this  value  by  h0  for  the  top  tube.  Values  obtained  are  then  com- 
pared to  Nusselt  theory  (Eq.  2-17),  the  Kern  model  (Eq.  2-19),  or  the  Eissenberg  model 
(Eq.  2-12).  For  both  the  local  and  average  values,  the  flatter  the  curves,  the  less  the  ef- 
fect of  condensate  inundation. 

The  effects  of  inundation  on  h0  for  the  smooth  copper  tubes  are  shown  in 
Figure  51,  Figure  52,  Figure  53  and  Figure  54  at  three  arbitrarily  chosen  coolant  ve- 
locities. As  expected,  there  is  virtually  no  difference  on  the  local  or  average  values  of  the 
heat  transfer  coefficient  when  the  tape  insert  is  used  as  opposed  to  the  HEATEX  ele- 
ment. In  addition,  the  values  for  the  average  heat  transfer  coefficient  lie  in  between  the 
values  predicted  by  the  Kern  and  Eissenberg  models  (equations  2-19,  and  2-21).  Values 
for  the  local  heat  transfer  coefficient  lie  slightly  above  those  predicted  by  the  Kern  model 
(equation  2-20).  These  results  suggest  that  Nusselt  theory  is  too  conservative  with  re- 
gard to  the  effect  of  condensate  inundation.  It  is  interesting  to  note  that  the  main 
inundation  effect  occurs  for  the  second  tube  with  little  additional  effect  for  subsequent 
tubes  in  the  bundle.  This  again  supports  the  notion  that  the  smooth  copper  tubes  can 
only  support  a  limited  amount  of  condensate 
2.     KORODENSE  Tubes 

Figure  55  and  Figure  56  show  values  for  h0  for  the  KORODENSE  tubes  with 
the  tape  inserts  installed.  The  figures  show  data  obtained  when  the  tubes  were  operated 
as  a  bundle  and  individually  respectively.  The  values  of  h0  for  the  top  tube  are  above  that 
predicted  from  Nusselt  theory  which  is  due  to  the  small  enhancement  provided  by  the 
roping  of  the  tubes.  This  enhancement  is  approximately  15  to  20%  in  comparison  to  the 
smooth  copper  tubes.  As  with  the  smooth  copper  tubes,  there  is  again  a  substantial  ef- 
fect of  inundation  (on  the  order  of  20%  for  the  second  tube  compared  to  the  top  tube). 
There  is  also  more  separation  between  the  values  of  h0  for  tubes  B,  C,  and  D  than  for 
the  smooth  copper  tubes  indicating  that  the  lower  tubes  may  not  be  "saturated"  with 
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Figure  52.      Local  h9  for  HEATEX  element. 
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condensate.  Finally,  note  that  the  relative  values  for  h0  for  the  KORODENSE  tubes  are 
about  10  to  15%  higher  than  the  corresponding  values  for  the  smooth  copper  tubes. 

The  data  for  the  KORODENSE  tubes  operated  individually  show  much  more 
scatter  than  was  seen  for  the  smooth  copper  tubes  especially  at  low  coolant  flow  rates. 
Also,  rather  than  being  virtually  identical  with  those  obtained  for  the  top  tube  during 
bundle  operation,  the  values  of  h0  lay  scattered  between  those  for  the  top  tube  (bundle 
operation)  and  the  values  predicted  by  Nusselt  theory.  Presumably,  this  is  due  to  either 
inaccuracies  or  uncertainty  in  the  inside  correlation  especially  at  low  coolant  (low  rates. 
At  higher  coolant  velocities,  the  scatter  appears  to  be  within  about  10%. 

Values  for  he  for  the  KORODENSE  tubes  obtained  when  the  HEATEX  ele- 
ments were  installed  are  shown  in  Figure  57  and  Figure  58  for  bundle  and  individual 
tube  operation  respectively.  As  with  the  experiments  conducted  with  the  tape  inserts 
installed,  values  of  h0  again  lie  above  those  predicted  by  Nusselt  theory.  The  enhance- 
ment appears  to  be  approximately  15  to  20%.  Again  there  is  a  substantial  inundation 
effect  which  is  of  the  order  of  20%  for  the  second  tube  in  comparison  to  the  top  tube. 
The  data  obtained  for  the  individual  tube  operation  shows  less  scatter.  However,  values 
for  all  the  tubes  again  lie  between  those  predicted  by  Nusselt  theory  and  those  obtained 
for  the  top  tube  during  bundle  operation. 

As  with  the  smooth  copper  tubes,  the  effects  of  changing  inserts  on  the  inside 
and  outside  heat  transfer  coefficients  are  shown  in  Figure  59  and  Figure  60.  The  inside 
heat  transfer  coefficients  are  identical  to  those  of  the  smooth  copper  tubes  (Figure  49 
and  Figure  50).  It  was  felt  that  the  enhancement  provided  by  the  roping  would  be  much 
less  than  the  enhancement  provided  by  the  inserts.  Therefore,  the  same  correlations 
were  used  to  determine  the  inside  heat  transfer  coefficient.  The  values  for  h0  are  slightly 
greater  for  the  KORODENSE  tubes  than  for  the  smooth  copper  tubes  (approximately 
5  to  10%). 

Figure  61  and  Figure  62  show  the  local  heat  transfer  coefficient  for  the  four 
tubes  for  the  tape  insert  and  the  HEATEX  elements.  There  is  virtually  no  difference 
between  the  two.  The  values  lie  above  the  Kern  estimate.  Figure  63  and  Figure  64  show 
the  average  heat  transfer  coefficient  for  the  four  tubes  with  both  inserts.  The  differences 
for  values  using  the  tape  inserts  as  opposed  to  the  HEATEX  elements  are  again  negli- 
gible. The  average  values  for  h0  also  lie  just  above  the  Kern  model.  This  is  in  contrast  to 
the  smooth  copper  tubes  where  the  average  values  lay  near  the  Nusselt  estimate.  This 
suggests  that  the  effect  of  condensate  inundation  is  indeed  less  for  the  KORODENSE 
tubes. 
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Figure  55.      KORODENSE  tubes,  tape  insert,  bundle  operation. 
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Figure  56.      KORODENSE  tubes,  tape  insert,  individual  operation. 
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Figure  61.      Local  h9  for  tape  insert. 
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3.  Wire  Wrapped  KORODENSE  Tubes 

The  effect  on  h0  of  wrapping  wire  around  the  KORODENSE  tubes  is  shown  in 
Figure  65,  Figure  66,  Figure  67  and  Figure  68.  In  viewing  these  figures,  one  should 
remember  that  the  values  for  the  unwrapped  KORODENSE  tubes  were  about  10% 
greater  than  predicted  from  Nusselt  theory  and  about  15%  greater  than  for  the  smooth 
copper  tubes.  From  the  figures,  it  appears  that  the  0.049"  diameter  wire  results  in  the 
greatest  enhancement  in  h0  for  all  four  tube  positions.  However,  the  relative  magnitude 
of  the  enhancement  decreases  as  lower  tubes  are  examined.  Since  this  data  was  obtained 
during  individual  tube  operation,  the  results  are  not  due  to  condensate  inundation. 

Values  for  h0  are  shown  for  each  wire  diameter  in  all  four  tube  positions  in  Fig- 
ure 69,  Figure  70,  Figure  71  and  Figure  72.  When  no  wire  is  used,  the  results  are  the 
same  as  those  shown  in  Figure  58.  For  the  0.029"  wire,  positions  A  and  B  are  essentially 
identical  with  tube  C  showing  the  greatest  enhancement  and  tube  D  the  worst.  The 
magnitude  of  the  deviation  is  about  20%.  The  0.049"  diameter  wire  results  in  the  best 
enhancement  in  h0  regardless  of  position.  However,  note  that  position  D  is  not  over- 
whelmingly greater  than  for  the  0.029"  wire.  The  data  show  that  the  level  of  enhance- 
ment depends  upon  bundle  position.  The  top  tube  positions  behave  similarly  while 
positions  C  and  D  show  reduced  levels  of  enhancement  Finally,  the  0.0675"  diameter 
wire  yields  values  for  h0  which  are  nearly  the  same  although  a  positional  effect  appears 
evident.  An  interesting  comparison  are  the  values  for  h„  for  the  0.049"  and  0.0675"  wire 
in  the  fourth  position.  Here  it  appears  that  the  0.0675"  wire  gives  roughly  the  same  level 
of  enhancement  especially  at  high  coolant  velocities. 

4.  KORODENSE  Tubes  Wrapped  with  0.049"  Wire 

The  effects  on  h„  of  wrapping  0.049"  diameter  wire  on  KORODENSE  tubes  on 
are  shown  in  Figure  73  and  Figure  74.  These  two  experiments  are  identical  runs.  In 
both  cases,  there  is  a  substantial  improvement  in  h„  for  all  tubes  when  compared  to  val- 
ues predicted  from  Nusselt  theory.  It  is  interesting  to  note  that  even  though  tubes  B, 
C,  and  D  are  all  being  inundated,  the  values  for  hc  are  still  well  above  Nusselt  theory 
(32%  for  tube  A,  28%  for  tube  B,  22%  for  tube  C  and  16%  for  tube  D).  In  comparison 
with  the  non-wire  wrapped  KORODENSE  tubes,  wrapping  with  the  0.049"  diameter 
wire  yields  an  increase  of  about  70%  for  tubes  A  and  B  and  about  a  30  to  50%  increase 
for  tubes  C  and  D.  Comparison  with  the  top  smooth  tube  yields  almost  a  90%  increase 
in  he. 

The  results  for  the  single  tube  runs  are  shown  in  Figure  75  and  Figure  76.  Both 
runs  demonstrate  considerable  scatter  in  the  data  of  about  +  22%.   It  is  also  interesting 


112 


HZERO  WIRE  WRAPPED  KORODENSE 


o 
6 

11 

OBt 

A,  K-l 

lio,  ni 

jA1£jA 

P5 

A 

© 

d 

LEGEND 
□  =  NO  WIRE 
o  =  0.029"  WIRE 
a  =  0.049"  WIRE 
+  =  0.0675"  WIRE 

N 

+ 

A; A 

© 

d 

HZERO  (W/M**2  K) 

)0.0             1700.0            191 

A 

A 

A 

A 

A 

A 

A 

6 

1      * 

-b 

9 

+ 

+ 

+ 
0 

+ 

0 

+ 
O 

m 

a 
6 

□ 

0 

o 

d 

D 

D 
□  ': 

a 

D 

1 

0 

D 

1 

a 



a 

0.1  0.2  0.3  0.4 


0.6  0.6  0.7  0.8  0.9 

VELOCITY  (M/SEC) 


1.0  1.1 


1JS  1.3 


Figure  65.      Effect  of  wire  diameter  on  h0  for  tube  A. 
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Figure  66.      Effect  of  wire  diameter  on  h0  for  tube  B. 
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Figure  67.      Effect  of  wire  diameter  on  h0  for  tube  C. 
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Figure  68.      Effect  of  wire  diameter  on  h0  for  tube  D. 
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Figure  69.      Effect  of  no  wire  on  h„  for  all  tubes. 
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Figure  70.      Effect  of  0.029"  wire  on  h0  for  all  tubes. 
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Figure  71.      Effect  of  0.049"  wire  on  h0  for  all  tubes. 
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Figure  72.      Effect  of  0.0675"  wire  on  h0  for  all  tubes. 
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Figure  73.      0.049"  wire  wrapped  KORODENSE  tubes,  bundle  operation,  expt.  1. 
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Figure  74.      0.049"  wire  wrapped  KORODENSE  tubes,  bundle  operation,  expt.  2. 
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to  note  that  the  data  for  these  single  tube  runs  is  somewhat  higher  than  for  the  bundle 
runs  (of  the  order  of  11%).  This  may  be  due  to  uncertainty  in  the  correlations  used  for 
/?,  which  are  reflected  in  h0 

The  values  for  the  average  and  local  vapor-side  heat  transfer  coefficients  are 
shown  in  Figure  77,  Figure  78,  Figure  79  and  Figure  80  for  the  two  experiments  con- 
ducted with  the  0.049"  wire  wrapped  KORODENSE  tubes.  The  average  vapor-side  heat 
transfer  coefficient  is  equal  to  that  predicted  by  the  Eissenberg  model  whereas  the  local 
vapor-side  heat  transfer  coefficient  is  substantially  above  the  Kern  estimate.  Recall  that 
for  the  plain  KORODENSE  tubes,  the  average  heat  transfer  coefficient  was  above  the 
Kern  estimate  but  below  that  predicted  by  the  Eissenberg  model.  The  local  heat  transfer 
coefficient  was  just  above  the  Kern  estimate.  The  differences  between  the  plain  and 
wire-wrapped  KORODENSE  tubes  demonstrate  that  wrapping  the  tubes  with  wire 
substantially  reduces  the  effects  of  condensate  inundation. 
5.     Copper/Nickel  Finned  Tubes 

The  values  for  h0  obtained  from  experiments  on  copper/nickel  finned  tubes  using 
the  tape  inserts  are  shown  in  Figure  81.  The  values  obtained  for  the  top  tube  are  ap- 
proximately 7  times  greater  than  predicted  by  Nusselt  theory.  Condensate  inundation 
caused  about  a  35%  decrease  in  h0  for  the  second  tube. 

The  values  for  h0  obtained  using  the  same  tubes  but  with  the  HEATEX  elements 
are  shown  in  Figure  82.  In  comparison  to  the  values  obtained  for  the  tape  inserts  (see 
Figure  81),  this  data  appears  very  scattered  with  enhancement  over  Nusselt  theory  of 
the  order  of  10  to  15  fold.  This  is  a  large  discrepancy.  It  is  possible  that  the  correlation 
for  the  small  diameter  HEATEX  element  may  be  suspect.  The  data  supplied  by  CAL 
GAVIN  was  collected  on  smooth  copper  tubes  with  a  diameter  similar  to  that  of  our 
copper/nickel  finned  tubes.  This  data  was  collected  over  a  lower  coolant  velocity  range 
than  that  used  in  this  study.  Therefore,  the  correlation  for  this  HEATEX  element  may 
not  be  valid  for  the  finned  tube  over  the  range  of  coolant  velocities  used  here.  This  area 
needs  to  be  studied  further,  perhaps  even  to  gathering  new  experimental  data  which  can 
be  used  to  develop  a  new  correlation. 

The  degree  of  inundation  can  be  seen  in  Figure  83,  Figure  84,  Figure  85  and 
Figure  86.  The  data  for  the  average  vapor-side  heat  transfer  coefficient  for  both 
HEATEX  as  well  as  tape  inserts  lie  between  the  Nusselt  and  Kern  models.  This  is  also 
true  for  the  local  vapor-side  heat  transfer  coefficients. 

An  interesting  comparison  is  to  look  at  the  average  and  local  heat  transfer  val- 
ues for  all  tubes  tested.  Regarding  the  average  heat  transfer  coefficient,  for  the  smooth 
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Figure  75.      0.049"  wire  wrapped  KORODENSE,  individual  runs,  expt.  1. 
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Figure  76.      0.049*  wire  wrapped  KORODENSE,  individual  runs,  expt.  2. 
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Figure  77.      Average  h.  wire  wrapped  KORODENSE,  expt.  1. 
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Figure  78.     Average  h0  wire  wrapped  KORODENSE,  expt.  2. 
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Figure  79.      Local  h,  wire  wrapped  KORODENSE,  expt.  1. 
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Figure  80.      Local  h9  wire  wrapped  KORODENSE,  expt.  2. 
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Figure  81.      h0  for  copper/ nickel  finned  tubes,  tape  insert. 
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Figure  82.      h0  for  copper/nickel  finned  tubes,  HEATEX  insert. 
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copper  tubes,  the  values  lay  between  the  Kern  and  Eissenberg  estimates.  The  plain 
KORODENSE  tubes  had  values  similar  to  those  for  the  smooth  copper  tubes  while  the 
0.049"  wire-wrapped  KORODENSE  had  values  which  matched  estimates  based  on 
Eissenberg's  model.  The  copper/nickel  finned  tubes,  on  the  other  hand,  had  values  which 
were  between  the  Kern  and  Nusselt  estimates.  Since  the  lower  the  values  the  more 
condensate  effect  present,  this  data  would  imply  that  the  effect  of  inundation  is  worse 
on  the  finned  tubes  and  least  on  the  wire-wrapped  tubes.  If  the  data  for  the  local  heat 
transfer  coefficients  is  analyzed  in  a  similar  manner,  the  same  conclusions  are  reached. 

The  wire  wrapped  KORODENSE  tubes  have  the  least  effect  of  condensate 
inundation  presumably  because  the  wire  serves  to  create  a  low  pressure  area  at  the  base 
of  the  wire.  This  causes  condensate  to  be  pulled  along  the  tube  which  enhances  heat 
transfer.  The  values  for  the  finned  tubes  are  not  in  agreement  with  those  of  Honda  et 
al.  [Ref.  42].  These  investigators  found  that  the  degree  of  inudation  effect  was  substan- 
tially less  than  that  reported  here.  There  are  two  apparent  reasons  for  this.  The  pitch- 
to-diameter  ratio  used  by  Honda  et  al.  was  1.37  whereas  in  the  present  study  it  is  2.25. 
This  difference  may  cause  differences  in  the  way  condensate  drains  from  one  tube  onto 
another.  Alternatively,  the  correlations  used  to  determine  the  heat  transfer  coefficient 
may  be  suspect.  Hong  and  Bergles  [Ref.  47]  have  reported  that  their  correlation  for  the 
twisted  tape  has  an  uncertainty  of  about  15%.  Therefore,  it  is  possible  that  the  large 
effect  of  condensate  inundation  is  due  more  to  uncertainty  in  the  correlations  than  to  a 
real  phenomenon. 

In  order  to  assess  the  validity  of  the  data  obtained  on  the  copper/nickel  finned 
tubes,  a  comparison  was  made  with  the  model  of  Beatty  and  Katz  [Ref.  13].  In  order  to 
do  this  comparison,  the  fin  efficiency  had  to  be  determined  using  equations  from 
Incropera  and  DeWitt  [Ref.  55].  It  was  assumed  that  the  fms  could  be  approximated 
as  rectangular  fins  with  adiabatic  tips.  The  fin  efficiency  was  then  approximated  by: 

tanh(mLc) 
""        ,nLc  {6~l) 

According  to  Incropera  and  DeWitt,  the  term  mLc  can  be  approximated  by  [Ref.  55]: 

^-{■$f£  (6-2> 
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Figure  83.      Average  //„  for  Cu/Ni  finned  tubes,  tape  insert. 
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Figure  84.      Average  h0  for  Cu/Ni  Tinned  tubes,  HEATEX  insert. 
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Figure  85.      Local  h0  for  Cu/Ni  finned  tubes,  tape  insert. 
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Figure  86.      Local  h„  for  Cu/Ni  finned  tubes,  HEATEX  insert. 
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where  Le  is  defined  as: 


Lc  =  L  +  j-  (6-3) 


1 

0.25        ,   ,„        Afs         1                    Afi         1               Au          1 

Here,  L  is  defined  as  the  height  of  the  fin  and  t  is  the  fin  thickness.  The  term  Ap  is  de- 
fined as: 

Ap  =  Lct  (6  -  4) 

In  order  to  calculate  a  value  for  the  fin  efficiency,  rjf,  an  iterative  procedure  was 
used.  First,  a  value  for  h0  was  guessed.  Based  on  this  value,  rjf  was  calculated  using 
equation  (6-1).  An  effective  area  was  then  calculated  from  the  following  equation  (Marto 
[Ref.  3]): 

Ae/=  VfAfs  +  rt/Aft  +  AU  (6  -  5) 

where  Ah  is  the  surface  area  of  the  fins  sides,  A/t  is  the  surface  area  of  the  fin  tips  and 
Au  is  the  interfin  surface  area.  Using  these  values,  the  effective  diameter,  D,  was  deter- 
mined from  (Marto  [Ref.  3]): 


(6-6) 


where  L  is  defined  as: 

-       n{Dl  -  D2r) 

L=        AD.  (6-7) 

The  heat  transfer  coefficient  is  then  calculated  from  a  modified  version  of  the  Beatty  and 
Katz  model  (Eq.  2-2)  which  uses  heat  flux  rather  then  the  film  temperature  difference, 
AT. 

(  k2p2ghfg  \j- 

The  results  of  these  calculations  for  the  copper/nickel  finned  tubes  with  the  tape 
inserts  are  shown  in  Figure  87.  In  this  figure,  the  values  of  h„  predicted  from  the  Beatty 
and  Katz  model  have  been  expressed  in  terms  of  the  smooth  tube  area,  that  is  nd„L  . 
Also  note  that  values  for  ha  have  been  calculated  for  copper,  copper/nickel  and  titanium. 
It  is  important  to  note  that  the  raw  data  lies  below  the  predicted  values  for  copper/nickel 
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finned  tubes.  The  data  obtained  after  reprocessing  with  the  modified  Wilson  plot  proce- 
dure lies  somewhat  above  the  prediction.  However,  both  lie  within  ±  20%  which  is 
within  the  range  of  validity  for  the  Beatty  and  Katz  model.  Therefore,  it  would  appear 
that  the  data  obtained  in  the  present  study  for  the  copper/nickel  finned  tubes  with  the 
tape  inserts  are  reasonable. 

The  results  for  the  copper/nickel  finned  tubes  with  the  HEATEX  elements  in- 
stalled are  not  as  good.  These  results  are  shown  in  Figure  88.  This  figure  shows  the 
tremendous  amount  of  scatter  in  the  data  irrespective  of  whether  the  raw  data  or  the 
data  obtained  from  the  modified  Wilson  plot  procedure  is  used.  In  either  case,  it  is  ap- 
parent that  the  data  do  not  match  values  predicted  from  the  model  of  Beatty  and  Katz. 
The  only  difference  in  these  two  experiments  is  the  type  of  insert  used.  The  values  for 
h0  are  based  on  an  experimentally  determined  value  for  U0  and  a  value  for  h,  which  is 
obtained  from  a  correlation.  In  this  study,  the  correlations  for  the  HEATEX  elements 
were  based  on  empirical  data  supplied  from  the  manufacturer.  This  data  was  obtained 
on  smooth  copper  tubes  having  the  same  diameter  as  the  copper/nickel  finned  tubes  used 
in  this  study.  Another  point  which  must  be  considered  is  that  the  range  of  velocity  for 
which  the  empirical  data  was  supplied  was  from  0.2  to  1.2  m/s.  However,  for  the 
copper/nickel  finned  tubes  used  in  this  study,  a  coolant  velocity  of  0.7  to  2  m/s  was  used. 
These  factors  would  suggest  that  the  correlation  used  for  the  small  diameter  HEATEX 
element  is  in  error.  Further  research  will  have  to  verify  whether  or  not  this  is  the  case. 

D.     IMPLICATIONS  OF  PRESENT  WORK 

This  study  used  R-113  as  the  refrigerant  and  a  refrigerated  ethylene  glycol/water 
mixture  as  the  coolant.  However,  Naval  refrigeration  systems  use  R-114  for  the 
refrigerant  and  sea  water  as  the  condensing  medium.  Therefore,  one  must  ask  whether 
the  results  obtained  from  the  present  work  can  be  extended  to  Naval  applications.  In 
order  to  perform  the  following  analysis,  the  dimensions  of  the  smooth  copper  tubes  and 
copper/nickel  finned  tubes  (26  fpi)  used  in  this  study  were  assumed.  The  inside  heat 
transfer  coefficient  was  calculated  using  the  Dittus-Boelter  correlation.  For  the  smooth 
tubes,  h0  was  calculated  using  the  Nusselt  equation,  modified  to  use  heat  flux  rather  than 
the  film  temperature  difference,  AT .  For  the  finned  tube,  h(  was  again  calculated  using 
the  Dittus-Boelter  correlation.  However,  h0  was  calculated  using  the  modified  Beatty  and 
Katz  equation  (Eq.  6-8).  The  fin  efficiency  was  assumed  to  be  1.0.  The  condenser  design 
parameters  were  modeled  after  the  requirements  for  the  DDG-51.  Hence,  the  refriger- 
ation load  was  fixed  at  242.7  tons.  The  number  of  condenser  tubes  was  fixed  at  340 
tubes.  Coolant  inlet  temperature  was  assumed  to  be  65  F  (18.3  C)  while  the  saturation 
temperature  of  R-114  was  107.8  F  (42.1  C).  The  coolant  velocity  was  held  constant  at 
2  m/s. 
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Figure  87.      Experimental  data  (tape  insert)  vs.  Beatty/Katz  model. 
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Figure  88.      Experimental  data  (HEATEX  insert)  vs.  Beatty/Katz  model. 
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Ratios  of  the  surface  areas  for  fumed  to  smooth  tube  and  a  finned  KORODENSE 
tube  to  smooth  tube  were  calculated.  In  order  to  carry  out  this  calculation  for  the 
KORODENSE  tube,  some  type  of  correlation  must  be  used  to  calculate  hr  However, 
as  shown  by  Cunningham  and  Ben  Boudinar  [Ref.  56]  in  their  work,  on  steam 
condensation  on  roped  tubes,  the  enhancement  due  to  the  roping  depends  upon  the  type 
of  groove  cut  into  the  tube.  In  this  analysis,  an  enhancement  of  1.6  for  h,  was  assumed 
for  the  KORODENSE  tube  over  the  finned  tube.  Once  the  area  ratio  is  obtained,  a 
condenser  tube  length  can  be  calculated.  This  is  done  as  follows.  The  area  ratio  can  be 
expressed  in  terms  of  length  by: 


A 


of 


lOJ 


Npid0L5 


(6-9) 


where  A0j  is  the  area  of  the  finned  tube,  A0iS  is  the  area  for  the  smooth  tube,  and  N  is  the 
number  of  tubes.   The  term  af  is  given  by: 


af={Af+A0)N 


(6-10) 


where  Af  and  Ae  are  the  finned  and  unfmned  areas  respectively. 

Table  3  shows  the  results  of  this  calculation.  For  the  finned  tube,  the  area  ratio  of 
1.75  leads  to  a  decrease  in  the  condenser  length  of  about  50%.  On  the  other  hand,  if  a 
roped  KORODENSE  tube  could  be  manufactured  which  had  fins  on  the  outside  surface, 
the  reduction  in  length  of  the  condenser  tubes  would  be  58%.  This  could  lead  to  even 
greater  weight  savings  for  refrigeration  systems. 

Table  3.     SUMMARY  OF  AREA  RATIO  CALCULATIONS 


Smooth 
tube 

Finned 
Tube 

KORODENSE  +  Fins 

*« 

K, 

1.0 

1.75 

1.51 

1.0 

0.497 

0.42 
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VII.     CONCLUSIONS 

The  following  conclusions  can  be  made  based  on  the  results  obtained  in  this  study. 

A.  GENERAL 

1.  The  condenser/evaporator  test  platform  is  fully  operational  and  shows  no  signs  of 
the  contamination  problems  reported  by  Mabrey  [Ref.   2  ]. 

2.  Values  for  outside  heat  transfer  coefficient  obtained  using  smooth  copper  tubes 
appear  to  be  reasonable  and  agree  with  other  published  data. 

B.  OVERALL  HEAT  TRANSFER  COEFFICIENT 

1.  KORODENSE  tubes  provide  enhancements  of  6  to  8%  over  smooth  copper  tubes 
regardless  of  the  type  of  insert  used. 

2.  For  the  KORODENSE  tubes  used  in  this  study,  wire  wrapping  results  in  an  en- 
hancement in  heat  transfer. 

3.  The  optimal  pitch-to-wire  diameter  was  about  7  for  the  KORODENSE  tubes  used 
in  the  present  work. 

4.  Wrapping  the  KORODENSE  tubes  with  0.049"  diameter  wire  resulted  in  an  en- 
hancement of  16%  over  the  non  wire-wrapped  KORODENSE  tubes  and  a  22  to 
24%  enhancement  over  the  smooth  copper  tubes. 

5.  Copper/nickel  finned  tubes  (26  fpi)  yielded  a  doubling  of  the  heat  transfer  in  com- 
parison to  smooth  copper  tubes. 

6.  The  use  of  HEATEX  radial  mixing  elements  yielded  a  10  to  15%  increase  in  com- 
parison to  the  use  of  twisted  tape  inserts. 

C.  OUTSIDE  HEAT  TRANSFER  COEFFICIENT 

1.  The  smooth  tubes  yielded  values  of  h0  which  were  7  to  15%  below  those  estimated 
from  Nusselt  theory. 

2.  The  KORODENSE  tubes  gave  values  of  h0  which  were  about  15%  greater  than 
those  estimated  from  Nusselt  theory  and  approximately  25%  greater  than  values 
obtained  for  the  smooth  copper  tubes. 

3.  Wrapping  the  KORODENSE  tube  with  0.049"  diameter  wire  resulted  in  a  90% 
increase  in  the  outside  heat  transfer  coefficient  in  comparison  to  the  smooth  copper 
tubes. 

4.  The  copper/nickel  finned  tubes  yielded  a  7  fold  increase  in  h0  over  the  smooth 
copper  tubes. 

D.  EFFECTS  OF  CONDENSATE  INUNDATION 

1.  For  smooth  copper  tubes,  condensate  inundation  resulted  is  a  15  to  18%  decrease 
in  h0  for  the  second  tube  compared  to  the  top  tube. 
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2.  KORODENSE  tubes  showed  a  20%  decrease  in  h0  for  this  same  comparison. 

3.  Wrapping  the  KORODENSE  tubes  with  0.049"  diameter  wire  resulted  in  a  decrease 
of  only  5  to  10%  in  h0  for  the  second  tube  compared  to  the  top  tube. 

4.  For  the  copper/nickel  finned  tubes,  condensate  inundation  resulted  in  almost  a 
35%  decrease  in  h0  for  the  second  tube  relative  to  the  top  tube.  This  large  decrease 
is  presumably  due  to  uncertainty  in  the  correlations  used  to  determine  the  heat 
transfer  coefficient. 


143 


VIII.     RECOMMENDATIONS 

Based  on  the  results  of  this  study,  the  following  recommendations  are  made: 

A.  GENERAL  RECOMMENDATIONS 

1.  The  ball  valves  which  control  coolant  flow  to  the  flowmeters  should  be  replaced 
with  needle  valves.  This  will  provide  better  control  of  coolant  flow. 

2.  The  flowmeters  should  be  replaced  with  larger  capacity  meters.  This  will  allow  in- 
vestigation of  heat  transfer  at  higher  coolant  velocities  and  therefore  higher 
condensation  rates.  This  change  will  have  to  be  accompanied  by  an  increase  in  the 
diameter  of  the  coolant  return  header  in  order  to  accommodate  the  increased 
coolant  flow. 

3.  Extra  lines  from  the  top  ends  of  the  condenser  to  the  cold  trap  need  to  be  added. 
This  will  ensure  that  no  isolated  pockets  of  noncondesibles  exist  in  the  condenser. 

B.  EXPERIMENTAL  RECOMMENDATIONS 

1.  -The  correlation  for  the  small  diameter  HEATEX  element  must  be  re-investigated. 

This  may  require  obtaining  additional  data  from  CAL  GAVIN  or  designing  an 
experiment  at  the  Naval  Postgraduate  School  which  can  be  used  to  obtain  the  re- 
quired data. 

2.  The  correlation  used  for  the  twisted  tape  insert  also  needs  to  be  re-evaluated  in 
light  of  the  large  inundation  effects.  As  with  the  HEATEX,  this  may  require  that 
additional  experiments  be  carried  out  at  the  Naval  Postgraduate  School. 

3.  Once  a  reliable  correlation  is  obtained,  the  experiments  conducted  on  the 
copper/nickel  finned  tubes  need  to  be  repeated. 

4.  Experiments  need  to  be  conducted  using  both  titanium  finned  tubes  as  well  as 
titanium  KORODENSE  tubes.  This  data  can  then  be  compared  to  that  obtained 
with  the  copper/nickel  finned  tubes  and  the  copper/nickel  KORODENSE  tubes  to 
determine  whether  titanium  tubes  can  deliver  the  required  performance. 

5.  Smooth  copper  tubes  should  be  wrapped  with  different  wire  diameters  at  different 
pitches  in  order  to  determine  an  optimal  pitch-to-wire  diameter  ratio. 

6.  Experiments  need  to  be  conducted  on  wire  wrapping  of  KORODENSE  tubes  which 
have  a  different  pitch.  These  experiments  should  be  based  on  results  obtained  from 
the  smooth  copper  tube  experiments. 

7.  Experiments  need  to  be  repeated  using  R-114  and  possible  alternate  refrigerants 
that  the  Navy  is  interested  in  such  as  R-24. 

8.  A  KORODENSE  type  tube  which  is  roped  on  the  inside  and  finned  on  the  outside 
needs  to  be  developed  and  tested  since  calculations  have  shown  that  this  type  of 
tube  can  lead  to  smaller  condensers. 
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APPENDIX  A.     FLOWMETER  CALIBRATION 

The  flowmeters  used  to  measure  flow  through  the  main  condenser  tubes  were  cali- 
brated in  the  following  manner.  The  meters  were  disconnected  from  the  condenser  tubes. 
An  empty  55  gallon  drum  was  placed  on  a  scale  (resolution  0.5  lbs).  The  condenser 
pump  was  started  with  flow  through  the  bypass. 

Flowmeters  were  calibrated  at  0  °C  and  24  °C.  Each  flowmeter  was  initially  opened 
to  maximum  flow.  The  flow  was  then  throttled  back  to  10%.  The  weight  of  fluid  in  the 
drum  was  recorded  and  the  timer  started.  The  weight  was  then  recorded  at  the  end  of  the 
calibration  time.  This  procedure  was  repeated  in  10%  increments  up  to  100%  flow. 
Flowmeter  number  2  was  only  calibrated  up  to  90%  flow  since  a  mechanical  stop  pre- 
vented it  from  going  to  100%.  The  calibration  runs  were  then  repeated  from  the  maxi- 
mum flow  back  to  10%  flow  in  10%  increments.  The  mass  flow  rate  was  then  calculated 
in  kg/sec  based  on  an  average  of  the  two  readings.  The  data  for  the  flowmeter  cali- 
brations as  a  function  of  nominal  flow  rate  are  shown  in  Table  4. 


Table  4. 

MASS  FLOW  RATE  (KG/SEC)  AT  0  AND  24  C. 

Flow  Meter  A 

Flow  Meter  B 

Flow  Meter  C 

Flow  Meter  D 

Nominal 

Flow  Rate 

(%) 

o°c 

24°C 

ox 

24°C 

o°c 

24°C       , 

0°C 

24°C 

10 

0.0149 

0.0172 

0.0127 

0.0181 

0.0141 

0.0179 

0.0253 

0.0318 

20 

0.0363 

0.0403 

0.0308 

0.0408 

0.0374 

0.0391 

0.0513 

0.0613 

30 

0.0544 

0.0621 

0.0531 

0.0626 

0.0566 

0.0616 

0.0798 

0.0894 

40 

0.0"94 

0.0851 

0.0739 

0.0844 

0.0793 

0.0899 

0.1129 

0.1216 

50 

0.1003 

0.112 

0.0953 

0.1102 

0.1000 

0.1103 

0.1406 

0.1555 

60 

0.1179 

0.1339 

0.1134 

0.1295 

0.1211 

0.1318 

0.1728 

0.1823 

70 

0.1388 

0.1556 

0.1379 

0.1506 

0.1427 

0.1552 

0.1975 

0.2114 

80 

0.1563 

0.1769 

0.1533 

0.1715 

0.1639 

0.1716 

0.2261 

0.2409 

.    90 

0.1846 

0.2028 

0.1799 

0.1912 

0.1813 

0.1975 

0.2461 

0.2667 

100 

0.2019 

0.2151 

0.2034 

0.2183 

0.2788 

0.2924 
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A  least  squares  linear  regression  was  performed  for  each  flowmeter  at  each  of  the 
two  temperatures  for  which  the  calibration  runs  were  carried  out.  The  linear  regression 
equation  takes  the  following  form: 

m  =  a  +  bN 

where  a  is  the  y-intercept  in  kg/s,  b  is  the  slope  in  kg/s/%  and  N  is  the  flowmeter  setting 
in  %.   The  results  are  shown  in  Table  5. 

Table  5.     FLOWMETER  CALIBRATION  REGRESSION  RESULTS. 


Flow  Meter  A 

Flow  Meter  B 

Flow  Meter  C 

Flow  Meter  D 

o  °c 

24  CC 

o  °c 

24  °C 

o  °c 

24  °C 

o  °c 

24  °C 

Slope 

0.00206 

0.00225 

0.00204 

0.00216 

0.00208 

0.0022 

0.0028 

0.0029 

Y  intercept 

-0.00414 

-0.00329 

-0.00615 

-0.00138 

-0.0036 

-0.00211 

-0.00122 

0.0048 

In  order  to  obtain  flowrates  at  any  intermediate  temperatures,  a  simplified  interpo- 
lation procedure  was  used.  The  results  of  the  calibrations  for  the  four  flowmeteres  are 
shown  in  Figure  89,  Figure  90,  Figure  91,  and  Figure  92. 
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Figure  89.      Calibration  for  flowmeter  A. 
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Figure 


90.      Calibration  for  flowmeter  B. 
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FLOWMETER  C 
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Figure  91. 


Calibration  for  flowmeter  C. 
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Figure  92. 


Calibration  for  flowmeter  D. 
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APPENDIX  B.     EMF  VERSUS  TEMPERATURE 

The  EMF  read  from  the  thermocouples  must  be  converted  into  a  temperature  for 
pruposes  of  calculation.  This  is  accomplished  using  a  1th  order  polynomial  which  has  the 
general  form: 

T=C{0)  +  C{I)VI 

where  V  is  the  EMF  reading  in  volts,  and  I  has  values  from  1  to  7.  The  leading  coeffi- 
cients are  given  in  Table  6. 

Figure  93  shows  the  relationship  between  EMF  and  temperature. 


Table  6.     LEADING   COEFFICIENTS    FOR   CONVERSION   OF 
EMF  TO  TEMPERATURE. 


Leading  Coefficient 

Value 

C(0) 

0.10086 

C(l) 

25727.944 

C(2) 

-767345.8 

C(3) 

78025595.8 

C(4) 

-9247486589 

C(5) 

6.977x10" 

C(6) 

-2.66xl013 

C(7) 

3.9jc1014 
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EMF  VERSUS  TEMPERATURE 
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Figure  93.      EMF  versus  temperature. 
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APPENDIX  C.  PROPERTIES  OF  ETHYLENE  GLYCOL/WATER 

MIXTURE 

The  physical  and  thermodynamic  properties  of  the  ethylene  glycol/water  mixture 
were  determined  using  equations  and  figures  from  Cragoe  [Ref.  56]  and  Gallant  [Ref. 
57]. 

The  density  of  the  ethylene  glycol/water  mixture  is  calculated  in  lines  2490  to  2530 
of  the  computer  program.  Density  as  a  function  of  temperature  for  various  concen- 
trations of  ethylene  glycol  are  shown  in  Figure  94. 

The  kinematic  viscosity  of  the  ethylene  glycol/water  mixture  is  calculated  in  lines 
2370  to  2420.  Figure  95  shows  the  kinematic  viscosity  plotted  as  a  function  of  temper- 
ature for  various  ethylene  glycol  concentrations. 

The  heat  capacity  is  calculated  in  lines  2430  to  2480.  Figure  96is  a  plot  of  heat  ca- 
pacity as  a  function  of  temperature  for  different  ethylene  glycol  concentrations. 

The  thermal  conductivity  is  calculated  in  lines  2560  to  2610  of  the  program.  Values 
of  the  thermal  conductivity  as  a  function  of  temperature  are  plotted  for  various  ethylene 
glycol  concentrations  in  Figure  97. 
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Figure  94.      Density  of  ethylene  glycol/Mater  mixture. 
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Figure  95.      Kinematic  viscosity  of  ethylene  glycol/water  mixture. 
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Figure  96.      Thermal  heat  capacity  of  ethylene  glycol/water  mixture. 
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Figure  97.      Thermal  conductivity  of  ethylene  glycol/water  mixture. 


157 


APPENDIX  D.     PHYSICAL  AND  THERMODYNAMIC  PROPERTIES  OF 

R-113 

The  physical  and  thermodynamic  properties  of  R-113  were  calculated  based  on 
equations  supplied  by  Chapman  [Ref.  58]. 

The  heat  capacity,  Cp,  of  R-113  was  calculated  from  the  equation: 

Cp  =  (929.0  +1.037) 

where  T  is  in  °C  and  Cp  is  in  kJ/kg  K.   Figure  98  shows  Cp  as  a  function  of  temperature 
for  R-113. 

The  dynamic  viscosity  for  R-113  is  calculated  from  the  equation: 

503 
_5  (r-2.15) 

fi=  1.34*10  5x  10 

where  T  is  in  K  and  //  is  given  in  kg/m  s.  The  relationship  between  \i  and  temperature 
is  shown  in  Figure  99. 

The  thermal  conductivity  of  R-113  can  be  calculated  from: 

k  =  1.7308(4.846*10~2  -  7~x  6.57xl0~5) 

where  T  is  in  °F  and  K  is  in  W/m  K.    Figure  100  shows  the  thermal  conductivity  as  a 
function  of  temperature  for  R-113. 

The  density  of  R-113  can  be  calculated  from: 

p  =  1.62jc103  -  Tx  (2.219  +Tx  2.358jc10~3) 

where  T  is  in  °C  and  p  is  in  kgjrtf.  Density  of  R-113  is  shown  as  a  function  of  temper- 
ature in  Figure  101. 

The  latent  heat  of  R-113  can  be  calculated  from  the  equation: 

hfg  =  2326.1(70.55  -  t  x  (0.04838  +  1.26*10~4r)) 

where  T  is  in  °F  and  h^  is  in  J/kg.  Figure  102  shows  the  relationship  between  latent  heat 
and  temperature  for  R-113. 
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Figure  98.      Heat  capacity  versus  temperature  for  R-113. 
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DYNAMIC  VISCOSITY  OF  R-113 
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Figure 


99.      Viscosity  as  a  function  of  temperature  for  R-113. 
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THERMAL  CONDUCTIVITY  OF  R-113 
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Figure  100.      Thermal  Conductivity  of  R-l  13  as  function  of  temperature. 
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DENSITY  OF  R-113 
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Figure  101.      Density  as  a  function  of  temperature  for  R-113. 
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Figure  102.      Latent  heat  versus  temperature  for  R-113. 
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APPENDIX  E.     SAMPLE  CALCULATION 

A  sample  calculation  was  conducted  on  run  1  of  data  set  SM06.  This  run  was  carried 
out  on  smooth  copper  tubes  with  the  HEATEX  insert  installed.  The  sample  calculation 
was  performed  for  the  top  tube  only. 

Temperatures  were  recorded  as  emfs  using  copper/ constantan  thermocouples.  EiMF 
values  were  converted  to  temperatures  using  the  equation  shown  in  Appendix  B. 
Table  7  summarizes  the  results. 

Table  7.     EMF  AND  TEMPERATURES  FOR  RUN  1  DATA  SET  SM06. 


Thermocouple 

EMF  (microvolts) 

Temperature  (C) 

T(0) 

190S 

46.83 

T(l) 

1895 

46.54 

T(2) 

1909 

46.86 

T(3) 

1922 

47.16 

T(4) 

1918 

47.06 

T(5) 

1919 

47.099 

T(6) 

367 

9.44 

T(7) 

369 

9.496 

T(8) 

370 

9.52 

T(9) 

377 

9.7 

T(10) 

893 

22.52 

T(ll) 

893 

22.53 

T(12) 

831 

20.99 

T(13) 

831 

20.99 

T(14) 

798 

20.18 

T(15) 

798 

20.18 

T(16) 

853 

21.54 

T(17) 

853 

21.54 

164 


For  tube  A,  the  temperature  differential  is  calculated  by: 

7(10) +  7(11) 
A7=  '  -7(6)=  13.086  C 

The  average  temperature  is  then  calculated  from: 

Tavg  =  7(6)  +  AT  x  0.5  =  15.98  C 

Based  on  the  ethylene  glycol/water  mixture  physical  properties  (Appendix  C),  and  using 
and  average  temperature  of  15.98  C  and  an  ethylene  glycol  concentration  of  54%,  the 
following  fluid  properties  are  obtained: 

pc=  1071.3  Ag/m3 
Hc  =  0.0051  kglms 
C„e  =  3196.99  JjkgK 

kc  =  0.3742  W\m  k 
Prc  =  43.57 

Based  on  the  calibration  curves  for  flowmeter  A,  the  mass  flowrate  of  coolant  cor- 
responding to  a  setting  of  15.74  can  be  calculated.  A  simplified  interpolation  procedure 
is  used  to  account  for  temperature  differences.  Based  on  this  calculation,  we  obtain  a 
mass  flow  rate  of  0.2979  kg/s  for  this  run. 

The  coolant  velocity  can  be  calculated  from: 


c     pA 


where  the  inside  cross-sectional  area  is  given  by: 


A  n     J1 


and 


Al  =  -2-  (0.013259)2  =  0.000138  rn 


r/  0.02979  A,,m      , 

V<=  (1071.3X0.000138)=  °-20Im/5 
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The  coolant  Reynolds  number,  Ree,  is  calculated  from: 

(0.200(0.013259) 
Re.  =  - — '-  =  559.89 

4.76jc10"6 
The  heat  transferred  to  the  coolant  is  calculated  from: 

q  =  mCPiCAT=  (0.02979)(3196.99)(1 3.085) 

q  =  1246.19  W 
and  the  heat  flux  is  then  calculated  from  : 

where 

A0  =  nd0L  =  rc(0.015875)(1.2192)  =  0.0608  rn 
where  L  is  the  condensing  length  of  the  tube.  Based  on  this  equation,  we  obtain: 

«"--aoi^  =  204949,f'/m2 

In  order  to  calculate  the  inside  and  outside  heat  transfer  coefficients,  the  physical 
and  thermodynamic  properties  of  R-113  must  be  calculated  (Appendix  D).  These  fluid 
properties  are  based  on  a  reference  or  condensate  film  temperature.  In  order  to  obtain 
this  film  temperature,  the  wall  temperature  is  needed.  To  obtain  this,  an  iteration  is 
carried  out  with  an  initial  value  of  hB  assumed  (in  this  case  1000  Wjm2  K).  The  wall 
temperature  is  then  calculated  from: 

T  —  7" 

*  wall        *  sat 


ho 


20494.9 
1000 


Twatt  =  46.925  -      ./":     =  26.43°C 
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The  film  temperature  is  then  calculated  from: 


T  -         sat  2    -j- 

lf~      3      +3    lwall 


or 


46  925        2 
Tf=       ^       +  -J  (26.43)  =  33.26°C 

Based  on  the  equations  of  fluid  properties  for  R-113  (Appendix  D),  the  following  values 
were  obtained  for  the  condensate: 

hA  =  147036.5  J/kg 
k  =  0.0734  W\m  K 
p=  1544.34  Ag/m3 
H  =  6.029x10-*  kgjtns 

Using  these  fluid  properties  together  with  equation  (5-15),  we  obtain: 

h0=  1907.06  Wlm2K 

Using  an  iterative  procedure,  we  adjust  the  initial  guess  of  h0  by  taking  the  average  of 
the  calculated  and  guessed  value  of  h0.  The  wall  and  film  temperatures  are  then  recalcu- 
lated and  new  fluid  properties  for  R-113  are  obtained.  A  new  value  of  h0  is  then  calcu- 
lated. This  procedure  continues  until  the  newly  calculated  and  previosuly  calculated 
values  are  within  0.1%.   For  this  run,  we  obtained: 

h0=\25\  W\rn   K 
The  next  step  is  to  find  a  value  for  the  overall  heat  transfer  coefficient,  U0,  where: 

K.-     ""° 


where,  for  this  run, 


and 


0      LMTD 


LMTD  =  30.48°  C 
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U0=  23Q944g9  =  672.5  Wjm2K 

The  inside  heat  transfer  coefficient,  h„  can  be  obtained  from  a  correlation  for  the 
Nusselt  number  for  single  phase  coolant  flow  in  a  tube.  This  Nusselt  number  is  obtained 
from  a  correlation  (in  this  case)  for  the  HEATEX  element  where: 

Nuc  =  0.226{Rec)°-65(Prc)0A6 

Nuc  =  0.226(559.89)-65(43.75)-46  =  78.41 

and 

Nuckc       (78.41)(0.3742) 


hi         dt  0.013259 


^  =  2213.07  W\rnK 


The  outside  heat  transfer  coefficient  is  calculated  from  equation  (5-12)  where: 

1 


hn  = 


JL-A.    »  A 


The  wall  resistance,  Rm,  is  defmed  as: 

Rm=       2Lk 


0.015875 
In 


0.013259 
2(1.2192)(386) 


Rm=      -vn.rvw™^      =1.91xl()-*KIW 


Therefore, 


L      _  1 

^  -67149-  (0.013°259)(22513.07)  -('-"^XO.OeSOS) 


h0  =  1070.2  W\rn   K 
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Figure  103  is  a  copy  of  the  computer  printout  for  this  run.  The  small  differences 
seen  in  the  computer  calculations  versus  these  sample  calculations  are  the  result  of 
roundoff  error. 
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Data  set  number   =   1 


RESULTS  FOR  TUBE  NO.  1 

Mass  flow  rate  =  2 .  980E-02 

Inlet  temperature  =  9.443E+00 

Saturation  temp  ( Deg  C)  =  4.693E+01 

DELT   temp   Dif.  =  1.307E+01 

Log.  Mean  temp   Dif.  =  3.048E+01 

Heat  flux  =  2.048E+04 

Resistance  of  Tube  Metal=  1.913E-04 

Prandtl  number  =  4.359E+01 

Reynolds  number  =  5.B08E+02 

Reynolds  number  H8rB  S    =  S.745E+02 

Inside  h.t.c.  =  2.21GE+03 

Inside  NUSSULT  NO.  =  7.852E+01 

OVERALL  H.t.c.  ( Uo  )  =  6.720E+02 
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Figure  103.      Computer  printout  of  results  for  run  1  data  set  SM06. 
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APPENDIX  F.     UNCERTAINTY  ANALYSIS  FOR  U0 

Run  number  1  of  data  set  SM06  was  also  chosen  for  the  uncertainty  analysis.  The 
analysis  was  only  conducted  on  U0.  The  measured  and  calculated  parameters  found  in 
the  sample  calculation  (Appendix  E)  were  therefore  used  in  this  section.  In  order  to  as- 
sess the  effects  of  coolant  velocity  on  the  uncertainty  in  the  measurement  of  U0  ,  an 
uncertainty  analysis  was  also  performed  on  run  number  1 1  of  data  set  SM06.  This  run 
had  a  coolant  velocity  of  1.2  m/s  as  opposed  to  0.2  m/s  in  run  number  1. 

The  uncertainty  analysis  utilizes  the  procedure  suggested  by  Kline  and  McClintock 
[Ref.  54].  This  procedure  states  that  if 

R  =  R{x]x2x3 ,*„) 

then  the  uncertainty  in  R,  namely,  SR  is  given  by  the  equation: 

«.   <lB-SxY  +  (l!LSxY  + +  (1E-Sxn 

OXx         7         V   ox2  J  V   oxn 


0.5 


where  x„  is  the  measured  variable  and  5x„  is  the  uncertainty  in  the  measured  variable. 

The  uncertainty  in  the  temperature  measurements  is  due  to  uncertainty  in  the  volt- 
age measured  by  the  thermocouples.  For  the  thermocouples  used  corresponds  to  0.1 
°C.  The  uncertainty  in  in  this  thesis,  this  value  was  assumed  to  be  4  microvolts  or  0.1 
C.  The  uncertainty  in  A7",  Tsat  and  coolant  inlet  and  outlet  temperatures,  Tei  and  Tc0  , 
can  be  calculated  as  follows: 

£(10) +  7(11) 

■*  CO 


2 


Taking  the  partial  derivatives  we  obtain: 


likewise, 


djco    _  1 

dT(\0)  "  2 


°Tco     _   1 
dT{\\)  "  2 
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The  uncertainty  in  the  measurement  of  T„  is  thus: 


STco  = 


8Tr. 


dT{\0) 


ST)  + 


dT,. 


571(11) 


ST 


0.5 


Based  on  the  thermocouple  uncertainty,  ST  =  0.1  C,  we  obtain: 


STC0  =  0.0778  C 


The  coolant  inlet  temperature  is  calculated  from  a  single  thermocouple,  number  6  for 
tube  A.  The  uncertainty  in  this  measurement  is  calculated  as: 


dTn 


571(6) 


=  1 


and, 


STr 


8T,, 


6T(6) 


ST 


0.5 


or, 


STci  =  0.\  C 


Finally,  the  uncertainty  in  AT  is  calculated  from: 

SAT=[(STj  +  (STcff5 


and, 


6AT=  0.1347  C 


Coolant  velocity  is  calculated  from: 


V  = 


m 


c     PA 

In  this  run,  m  was  calculated  to  be  0.02979  kg/s.  The  uncertainty  in  the  reading  of  the 
flowmeter  is  defined  as  the  scale  interpolation.  This  value  corresponds  to  one-half  the 
value  of  the  smallest  marked  increment.  In  this  case,  the  scale  interpolation  factor,  Sf„ 
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is  0.5%  which  corresponds  to  0.00101  kg/s  resolution.  In  addition,  a  time-wise  jitter, 
StJ,  of  0.5%  was  noted  during  the  runs.   Therefore,  the  uncertainty  in  m  is  given  by: 


dm 
m 


8  m 


m 


+ 


5  m 


m 


or 


Srn 
m 


=  0.0479 


The  uncertainty  in  the  cross-sectional  area  of  the  tubes  was  estimated  from: 


where  Sd(  is  given  as  0.1  mm  based  on  tolerances  supplied  from  the  manufacturer. 
Therefore, 


SAt 


=  0.0107=  1.07% 


The  uncertainty  in  the  coolant  velocity  is  now  given  by: 


iff 


a-  hA       "l0-5 

{    m    }       K    A,    ' 


or 


bVr 


=  0.049  =  4.9  % 
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The  uncertainty  in  the  calculation  of  coolant  Reynolds  number  is  given  by: 


SRec 
Rec 


5Ve   2        Sdt 

(-F7£-)2  +  (-7L)2 


nO. 5 


*l 


or  substituting  appropriate  numbers  we  find 

SRec 


Rec 


=  0.0497  =  4.97  % 


The  uncertainty  in  the  calculation  of  the  heat  transferred  to  the  coolant  is  related 
to  the  uncertainty  in  the  measurements  of  coolant  mass  flow  rate  and  coolant  temper- 
ature rise  where: 


bq_ 


[,  dm 
m 


>2  +  (^F>2] 


0.5 


or 


6q 


=  0.0489  =  4.89  % 


The  uncertainty  in  the  heat  flux  is  then: 


<l'o 


n0.5 


where  the  uncertainty  in  A0  is  given  by: 


SA, 


bd„ 


{-^)2  +  {At)2 


no.5 


Here,  the  uncertainty  in  the  length  of  the  tube,  SL  ,  is  assumed  to  be  5  mm  while  5d„  is 
assumed  to  be  0.1  mm  (again  from  given  tolerances).  Substitution  of  the  appropriate 
numbers  yields: 

-A:2-  =  0.0496 
<7  o 
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Finally,  since  Uc  is  given  by: 


Un  = 


0      LMTD 


the  uncertainty  in  the  LMTD  must  be  calculated.   The  LMTD  is  given  by: 

LMTD  =  — - 


In 


T      —  T 

1  sat       l  c,l 

~f     —  T 

1  sat        l  c,o 


In  order  to  calculate  the  uncertainty  in  the  measurement  of  the  LMTD,  the  following 
equations  were  used: 


6LMTD 
8  AT 


1 


In 


T     —  T 

1  sat       J  c,l 

~f     —  T 

1  sat       A  c,o 


8LMTD ATl(Tsat  ~  Tc,o)  -  (Tsat  ~  Tcjfl 

dTsat 


In 


T    -  T  ■ 

1  sat        l  c,i 

~T    -  T 

1  sat       '  c,o 


\* sat       * c,i)\*  sat       *cj) 


8LMTD 
cT.„ 


AT 


In 


T     —  T 

1  sat        l  c,i 

~T     -  T 

1  sat        l  c,o 


\*sat        '  c,o) 


8  LMTD 


AT 


In 


T     -  T     ' 

1  sat        l  c,i 

~f     -  T 

1  sai        l  c,o 


\*  sat       *  c,i) 


The  uncertainty  in  the  LMTD  is  therefore  given  by: 


3  LMTD  = 


8 LMTD   8AJ\2  ,  (  8 LMTD   *-    \2  .  f  8 LMTD   ^    h  ,  /   8 LMTD   <T 
8AT  )   +  {~8T *"'    Hl""^ ST"I    ^\'~Tr-STcJ 


sat 


LM1D    ,r     \2   ,      _cLM_L 
dTc>0     0l")  +V      dTd 


For  the  sample  run  with  a  coolant  velocity  of  0.2  m/s,  the  uncertainty  in  the  measure- 
ment of  the  LMTD  is: 


21  2 


SLMTD  =  [[(2.329)(0.1347)]2  +  [(  -1.0157)(0.0778)]2  4-  [(  -2.9097)(0.0778)]2  +  [(1.894)(0.1)]2] 
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or 

SLMTD  =  0.446  K 
and 

TMrD-  =  "3a48"  =  0-0146=L46% 

The  uncertainty  in  the  overall  heat  transfer  coefficient,  Ue,  is  then  given  by: 


u0 


Wo  ^2 


<7  o 


(  SLMTD_  Y 
\    LMTD    ) 


or 


5  =  [(0.0496)2  +  (0.0146)2]T 


U0 


and 


—2-  =  0.052  =  5.2  % 

Table  7  summarizes  the  results  of  the  uncertainty  analysis  for  the  smooth  copper  tubes 
with  the  HEATEX  insert  for  coolant  velocities  of  0.2  and  1.2  m/s  respectively. 

The  uncertainty  analysis  neglects  any  uncertainty  in  the  physical  properties  of  the 
coolant  or  the  R-113.  For  the  most  part,  this  is  probably  a  reasonable  assumption.  The 
exception  to  this  may  be  the  viscosity  due  to  its  sensitivity  to  temperature.  Hence,  the 
uncertainty  reported  here  is  an  underestimate. 

It  is  curious  that  the  relative  magnitudes  of  the  uncertainty  for  the  low  and  high 
coolant  velocities  are  essentially  the  same.  Two  different  mechanisms  dominate  in  the 
calculated  uncertainty  at  low  and  high  coolant  velocities.  In  terms  of  the  mass  flow  rate, 
the  time-wise  jitter  and  scale  interpolation  terms  are  the  same  for  both  coolant  velocities. 
However,  the  actual  mass  flow  rate  is  less  at  the  lower  flow  rates.  Hence,  the  uncertainty 
in  the  mass  flow  rate  calculation  is  more  significant  at  the  lower  coolant  velocities.  At 
higher  flow  rates  the  uncertainty  in  the  calculation  of  the  LMTD  becomes  the  dominant 
term.  This  is  because  the  coolant  temperature  rise  is  smaller  at  high  coolant  velocities 
but  the  uncertainty  in  the  temperature  measurements  remains  the  same.  The  mass  flow 
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Table  8.     SUMMARY  OR  RESULTS  FOR  SAMPLE  UNCERTAINTY  ANALY- 
SIS. 


Variable 

Vc  =0.2  m/s 

Vt  =1.2m's 

*T« 

o.rc 

o.rc 

&T„ 

0.078°C 

0.078°C 

ST,B, 

0.078°C 

0.078°C 

dAT 

0.135°C 

0.135°C 

dm 
m 

0.048 

0.008 

SA, 

A, 

0.011 

0.011 

sve 

0.049 

0.013 

SRec 
Re, 

0.05 

0.015 

q 

0.05 

0.015 

5A, 

A, 

0.008 

0.008 

ft 

0.05 

0.017 

8LMTD 
LMTD 

0.015 

0.054 

0.052 

0.057 

rate  goes  up  so  that  the  uncertainty  in  the  coollant  mass  flow  rate  calculation  decreases. 
Thus,  at  the  higher  coolant  velocities,  the  uncertainty  in  the  calculation  of  the  LMTD 
is  the  dominant  term. 


177 


APPENDIX  G.     PROGRAM  LISTING 
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1000 

1005 

ube  b 

1015 

1020 

1021 

1022 

1023 

1024 

1025 

1030 

1035 

1040 

1045 

1050 

1055 

1060 

1055 

1070 

1075 

1080 

1085 

1090 

1095 

1100 

1105 

1110 

1115 

1120 

1125 

1130 

1135 

1140 

1145 

1150 

1155 

1160 

1215 

1216 

1217 

1218 

1220 

1225 

1230 

1235 

1240 

1245 


I   FILE: 
!   PURPOSE 
undle . 

CREATED 

UPDATE 1 

UPDATE2 

UPDATE3 

UPDATE4 

UPDATES 
CHANGE  H0D=3 

BEEP 

PRINTER  IS  1 

PRINT  USING 

PRINT  USING 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


DRPC0N7 
This  program 


collects  and  processes  condensat ion  data  for  the  t 


NOVEMBER  2,  1988 

APRIL  25  1991  (S.  MEMORY) 

May  17  1991  (S.  MEMORY) 

JULY  26,  1991  (R.W.  MAZZONE) 

13  AUGUST  199KR.W.  MAZZONE) 

10  OCTOBER  1991  (R.W.  MAZZONE) 

AND  USING  CALCULATED  CI  FROM  WILSON 


"4X, 
"BX, 
"6X, 
"6X, 
"6X, 
"BX, 
"BX, 
"6X, 


USING 
USING 
USING 
USING 
USING 
USING 
PRINTER  IS  701 
INPUT  Icall 
IF  Icall-0  THEN 
Icall-1  THEN 
Icall=2  THEN 
Icall=3  THEN 
Icall=4  THEN 
Icall=5  THEN 
Icall-B  THEN 


"SELECT  OPTION 

"0  TAKING  DATA  OR  REPROCESING  PREVIOUS  DATA 

"1  PLOTTING  H  US  DELTA-T 

"2  PLOTTING  HRAT  US  N 

"3  PLOTTING  WILSON 

"4  PURGE  FILES 

"5  XYREAD""" 

"6  NUSSELT  ESTIMATE 


IF 

IF 

IF 

IF 

IF 

IF 

END 

SUB 

COM 

COM 

COM 

COM 

COM 

DIM 

DATA 

DATA 

DATA 

DATA 

READ 

DATA 


CALL  Mam 
CALL  Plot2 
CALL  Plotl 
CALL  Plot3 
CALL  Purge 
CALL  Xyread 
CALL  Nusselt 


Main 

/Cc/  C(7) 

/Fid/  Ift 

/Nus/  Iin.Tsat ,Qdpl , Hnus , K f ,Rhof , Hf g , Muf , Do, I  tube 

/Will/  Doa(4),Dia(4),Kna(4),Iact  ,Droot(4) 

/Wil2/  Delta, I  sat ,Nset s , Hod,Cia( 3,3 ) , Alpaa( 3 ) 

Enf  <  17>,TpU),T<  17  ) ,  Ho(  4  )  ,Qdp(  4  )  ,Uo(  4  )  ,Pc(  4  ) 

1.0,1.0,1.0,1.0 

5.172,5.172,5.172,5.172 

226 
063 


,5, 
226,0. 
,063,0, 


0.226,0.226,0 

0.063,0.063,0 

Cia( *  ) 

0.10086091,25727.94369,-767345.8295,78025595.81 
DATA  -9247486589, 6. 97688E11 , -2 . 661 92E13 , 3 . 94078E1 4 
READ  C(«  ) 

DATA  0.01 5875 ,0.01 4000 ,0.015798,0.0159,0.0 
DATA  0.013259,0.010160,0.012776,0.01346,0.0 
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1250  DATA  386.0,42.975,42.975,21.9,0.0 

1251  DATA  0.0,0.014,0.0, .0141 ,0.0 
1255  DATA  0.0005588,3 

12G0  READ  Doa( *  ),Dia< *  ),Kma<*  ),Droot<»  ), Delta, Hod   !  Hod=H/Di 

12B5  L=1.2192   !  Condensing  length 

1270  Jset=0 

1271  Okaccpt=l 
1275  BEEP 

1280  INPUT  "ENTER  MONTH,  DATE  AND  TIME  < MM: DD  =  HH : MM= SS  ) " ,Dt g$ 

1285  OUTPUT  709; "TD" ;Dtg$ 

1290  BEEP 

1295  INPUT  "SELECT  OPTION  (0=DAQ,  1-FILE >" , Im 

1300  Ihard=l 

1305  BEEP 

1310  INPUT  "WANT  A  HARDCOPY  PRINTOUT  < 1=DEF=YES , 0=N  ) " , Ihard 

1315  BEEP 

1320  INPUT  "SELECT  ( 0=R-1 14 , 1=STEAM , 2-R-l 13 , 3=EG  ) " , I f t 

1325  Iin-1 

1330  Isat=2 

1335  BEEP 

1340  INPUT  "SELECT  SAT  TEMP  MODE  ( 0=LIQ , 1=VAP , 2  =  ( LIQ+UAP  )/2=DEF  ) " , Isat 

1345  IF  Ihard=l  THEN  PRINTER  IS  701 

1350  IF  lm=0  THEN 

1355  BEEP 

1360  INPUT  "GIUE  A  NAME  FOR  THE  NEW  DATA  FILE" ,File$ 

1365  CREATE  BDAT  File$,20 

1370  ASSIGN  ©File  TO  File$ 

1375  BEEP 

1380  INPUT  "ENTER  TUBE  CODE  (0=SMOOTH  CU,1=FINNED  CU ,2=KOR0DENSE , 3=TIT  FIN 
)",Itube 

1381  BEEP 

1382  INPUT  "0=NO  INSERT,  1=TWISTED  TAPE,  2=HITRAN1,  3=HITRAN2 ", Insert 
1385  BEEP 

1390  INPUT  "ENTER  EG  CONCENTRATION  ( WT  PERCENT  )" ,Egrat 

1395  ENTER  709;Dtg$ 

1400  OUTPUT  »File;Dtg$ 

1405  OUTPUT  @File;Itube,Egrat,Ddl,Dd2,Dd3,Dd4,Dd5 

1410  Iact=10 

1415  BEEP 

1420  INPUT  "SELECT  ( 0=TOP , l=SECOND 10=BUNDLE  =  DEFAULT  ) " , Iact 

1425  .      PRINT 

1430  PRINT  USING  "10X,""FILE  NAME:    " "  ,  12A" ; Fi le$ 

1435  PRINT 

1440  ELSE 

1445  BEEP 

1450  INPUT  "ENTER  NAME  OF  EXISTING  FILE",File$ 

1455  BEEP 

1460  INPUT  "ENTER  NUMBER  OF  DATA  SETS  STORED" ,Nset s 

1465  Iwil-1 

1470  BEEP 
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1475  INPUT  "WANT  TO  CALL  WILSON?  ( 1  =  DEFAULT=YES , 0=NO  ) " , I wi 1 

1480  IF  Iwil-1  THEN 

1485  BEEP 

1490  INPUT  "WHICH  TUBE  ( 0-TOP , l  =  SECOND , . . , 10=BUNDLE  )"  ,  lac t 

1491  INPUT  "0=NO  INSERT,  1=TWISTED  TAPE,  2=HITRAN1,  3=HITRAN2 ", Insert 

1493  CALL  Wilson 

1494  IF  Iact=10  THEN 

1495  FOR  lcount=0  TO  3 

14 9B  Cia< Insert, I count  )*Cia( Insert ,0 ) 

1497  Alpaa< Icount  )=Alpaa(0  ) 

1498  NEXT  Icount 

1499  END  IF 

1500  IF  Ihard=l  THEN  PRINTER  IS  701 
1505  PRINT 

1510  PRINT  USING  "10X,""FILE  NAME:    " " , 12A" ; Fi le$ 

1515  PRINT 

1520  IF  Iact>9  THEN 

1521  FOR  11=0  TO  3 

1522  PRINT  USING  " 10X ," "RESULTS  FOR  TUBE  NO . " " , 3D" ; 1 1  +  1 

1523  PRINT 

1525  PRINT  USING  " 10X ,"" INSIDE  COEFFICIENTS:  "", DD . 3D" ; Cia< Insert , II 
) 

1530  PRINT  USING  " 10X ," "ALPHAS :  " "  ,  4(  IX  ,  MZ . 3DE , 2X  )" ; Alp 
aa<  II  ) 

1531  NEXT  II 
1535  ELSE 

1540  PRINT  USING  " 10X ,"" INSIDE  COEFFICIENT  FOR  TUBE  " " , D , " " :  "",DD.3 
D" ; Iact  +  1 ,Cia( Insert , Iact  ) 

1545  PRINT  USING  "10X, ""ALPHA  FOR  TUBE  "",D,"":              "",DD.3 
D" ; Iact  +  1 ,Alpaa( Iact  ) 

1550  END  IF 
1555  ELSE 
15B0  BEEP 

1551  INPUT  "0=TOP,  l=SECOND,  10=BUNDLE " , lac t 

15G2  INPUT  "0=SMOOTH  CU ,  1-FIN  CU,  2=KOROD,  3=FIN  TITANIUM" , I  tube 

15B3  INPUT  "0=NO  INSERT,  1=TWISTED  TAPE,  2=HITRAN1,  3=HITRAN2 ", Insert 

15B5  INPUT  "ENTER  CI  VALUE  "  ,Cia( Insert , II  ) 

1570  END  IF 

1575  ASSIGN  ©File  TO  File$ 

1580  ENTER  @File;Dtg$ 

1585  ENTER  ©File; Itube,Egrat , Ddl , Dd2 , Dd3 , Dd4 , Dd5 

1586  INPUT  "0=SMOOTH  CU ,  1-CUNI  FIN,  2=KORODENSE,  3=FIN  TITAN  ",Itube 
1590  END  IF 

1595  Iout=l 

1600  BEEP 

1605  INPUT  "WANT  TO  CREATE  AN  OUTPUT  FILE?  ( 1=DEF  =  YES , 0=N0  ) " , lout 

1610  IF  Iout=l  THEN 

1615  BEEP 

1620  INPUT  "ENTER  A  NAME  FOR  OUTPUT  FILE",Fout$ 

1625  CREATE  BOAT  Fout$,5 
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1G30       ASSIGN  @Fout  TO  Fout$ 

1635  END  IF 
1640  ! 

1645  Do  =  Doa< Itube  ) 

1646  Dol=Doa( Itube ) 
1650  Di=Dia( Itube ) 
1655  Km=Ki*ia<  Itube) 

1660  Ax=PI*DiA2/4   !  Cross-sectional  area 

1665  Ao=PI*Do*L 

1670  !   Rn=Do*LOG<Do/Di )/(2*Km) 

1671  IF  Itube=l  THEN 

1672  Dol=Droot( Itube ) 

1673  END  IF 

1674  IF  Itube=3  THEN 

1675  Dol=Droot< Itube  ) 

1676  END  IF 

1677  Ao=PI*Do»L 

1679  Rn-LOG(Dol/Di >/<2*L*Kn> 

1680  Al=PI*Dol*L 

1682  ! 

1683  IF  In-0  THEN 
1685       PRINTER  IS  1 

1690  BEEP 

1691  INPUT  "WHAT  IS  REQUIRED  NOMINAL  COOLANT  VELOCITY  IN  m/s  < MAX .  -  1.4  M 
/s  )?" ,Uegst 

1694  OUTPUT  709; "AR  AF6  AL9  UR5" 

1695  FOR  1=6  TO  9 

1696  OUTPUT  709; "AS  SA" 

1697  Vsun=0 

1698  FOR  J=l  TO  5 

1699  ENTER  709;E 

1700  Vsum-Vsum+E 

1701  NEXT  J 

1702  Emf< I  >=Vsum/5 

1703  T( I  )=FNTvsv(Emf( I  )  ) 

1704  NEXT  I 

1705  Tavgl  =  (T(6)  +  T(7)  +  T(8)  +  T(9))/4 

1706  IF  Tavgl<0  OR  Tavgl>24  THEN 

1707  PRINT  "COOLANT  TEMP.  NOT  IN  RANGE  0  TO  24  C" 

1708  GOTO  1691 

1709  END  IF 

1711  Rhoegl=FNRhoeg(Tavgl,Egrat ) 

1712  Mdotl=Rhoegl*Ax*Vegst 

1713  Pca0=(Mdotl+. 00414)/. 00206 

1714  Pca24=(Mdotl+. 00329)/. 00225 

1715  Pcb0=(Mdotl+. 00615)/. 00204 

1716  Pcb24=(Mdotl+. 00138)/. 002 16 

1717  Pcc0=(Mdotl+. 0036)/. 00208 

1718  Pcc24=(Mdotl+. 00301  )/. 00223 

1719  Pcd0=(Mdotl+. 00122 )/. 0028 
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PRINT  USING  "BX,""%  OF  METER  A  = 

"",DDD. 

D" 

;Pc<0) 

PRINT  USING  "6X,""%  OF  METER  B  = 

■■    " "  ,DDD, 

,D" 

;Pc<  1  ) 

PRINT  USING  ,,EX,""%  OF  METER  C  = 

"",DDD. 

D" 

;Pc<2  ) 

PRINT  USING  "GX,""%  OF  METER  D  = 

■  "" ,DDD. 

,D" 

;Pc(3) 

PRINT 

PRINT  "HIT  CONTINUE  WHEN  READY" 

PAUSE 
I 

PRINTER  IS  701 

Dtld=47.5    !  Desired  temperatur 

e  of  lie 

jui 

d 

Ido-1 

ON  KEY  0,15  RECOVER  1734 

PRINT  USING  "4X, ""SELECT  OPTION 

ii  ti  ti 

PRINT  USING  "6X,""0=TAKE  DATA  "" 

" 

PRINT  USING  "GX,""1=SET  Tsat  (DEFAULT)  ' 

■  ii  ii 

PRINT  USING  "4X," "NOTE:  KEY  0  = 

ESCAPE"' 

i  H 

Ido  =  desired  option 

BEEP 

INPUT  Ido 

1720  Pcd24=(Mdotl-. 0048  )/. 0029 

1721  Pc(0)=Pca0+(Tavgl*(Pca24-Pca0)  )/24 

1722  Pc< 1  )=Pcb0+(Tavgl*(Pcb24-Pcb0)  )/24 

1723  Pc(2  )=Pcc0+(Tavgl*(Pcc24-Pcc0)  )/24 

1724  Pc<3)=Pcd0+(Tavgl*(Pcd24-Pcd0)  )/24 

1725  PRINT  "SET  FLOWMETER  READINGS  CORRESPONDING  TO: 
172G 
1727 
1728 
1729 
1730 
1731 
1733 

1734  Repeat : 

1735  ! 
173B 
1737 
1738 
1739 
1740 
1741 
1742 
1743  ! 
1744 
1745 
174G  !      Set  default  value  for  input 

1747  IF  Ido.M  THEN  Ido=l 

1748  !      Take  data  option 

1749  IF  ldo=0  THEN  1797 

1750  !      Loop  to  check  saturation  temperature 

1751  IF  Ido-1  THEN 

1752  INPUT  "ENTER  DESIRED  Tsat  ( DEFAULT=47 . 5  C  -  R-113)",Dtld 

1753  !         INPUT  "ENTER  DESIRED  Tsat  (DEFAULT=2.2  C  -  R-114  )"  ,Dt Id 

1754  Nn=l 

1755  Nrs=Nn  MOD  15 
175B  Nn=Nn+l 

1757  IF  Nrs=l  THEN 

1758  PRINT  USING  "4X,""  DTsat     Tldl      Tvab     Tvat      Dell 
el2     Del3     Del4 

1759  END  IF 

17B0  !        Read  thermocouple  voltages 

17G1  OUTPUT  709; "AR  AF0  AL17  UR5" 

1762  Nend=17 

17B3  FOR  1=0  TO  5 

17B4  OUTPUT  709; "AS  SA" 

17G5  Usum=0 

176B  FOR  J=l  TO  5 

1767  ENTER  709 ;E 

1768  Usum-Usum+E 

1769  NEXT  J 
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1770  Emf  < I  )=Vsum/5 

1771  T< I  )=FNTvsv(Emf< I  )  ) 

1772  NEXT  I 

1773  FOR  1=6  TO  Nend 

1774  OUTPUT  709; "AS  SA" 

1775  Vsum=0 

177B  FOR  J=l  TO  20 

1777  ENTER  709; E 

1778  Usum=Vsum+E 

1779  NEXT  J 

1780  Emf < I  )=Usum/20 

1781  T< I  >  =  FNTvsv(EmfU  )  ) 

1782  NEXT  I 

1783  !  Compute  average  temperature  of  liquid 

1784  Tliq-(T(3)  +  T(4)  )/2 

1785  Tvap  =  <T<0)  +  T<  1>+T(2  )  )/3 

1786  Tvat  =  (T(0)  +  T(  1  )  )/2 

1787  Tliql=T(5) 

1788  Dell  =  <(T<10)  +  T(ll  ))/2)-T(6  ) 

1789  Del2  =  <  (T(  12  )  +  T(  13)  )/2  )-T<7> 

1790  Del3  =  <  <T(  14  )  +  T<  15)  )/2  )-T(8) 

1791  Del4  =  <  <T<  IB  )  +  T(  17))/2  )-T(9> 

1792  PRINT  USING  "4X , 8( NOD . DO , 3X  ) " ; Dt Id, T( 3  ) , T< 2  ) , Tvat , Del  1 ,Del2 , Del3,De 
14 

1793  WAIT  2 

1794  GOTO  1755 

1795  END  IF 

1796  i  TAKE  DATA  IF  lm=0  LOOP 

1797  PRINTER  IS  701 

1798  OUTPUT  709; "AR  AF0  AL17  UR5" 

1799  Nend=17    !  INCREASE  TO  20  IF  FIVE  TUBES  IN  BUNDLE 

1800  FOR  1=0  TO  5 

1801  OUTPUT  709; "AS  SA" 

1802  Vsum=0 

1803  FOR  J=l  TO  5 

1804  ENTER  709;E 

1805  <Jsum=Vsum+E 

1806  NEXT  J 

1807  Emf < I >=Usum/5 

1808  NEXT    I 

1809  FOR  1=6  TO  Nend 

1810  OUTPUT  709; "AS  SA" 

1811  Usum=0 

1812  FOR  J=l  TO  20 

1813  ENTER  709;E 

1814  Vsum=V5um+E 

1815  NEXT  J 

1816  Emf  (I  )=k>sum/20 
1820  NEXT    I 

1825  ELSE 


184 


1830 
1835 
1840 
1845 
1850 
1855 
I860 
1865 
1870 
1875 
1880 
1881 
1885 
1890 
1895 
1896 
1900 
1905 
1906 
1907 
1908 
1925 
1926 
1930 
1935 
1940 
1945 
1950 
1955 
I960 
1965 
1970 
1975 
1980 
1981 
1985 
1990 
1991 
1992 
1994 
1995 
1996 
2000 
2005 
2010 
2015 
2020 
2025 
2030 
2035 


ENTER  @File;Pc(*),Enf<#),Tp(*) 

END  IF 

DATA  ANALYSIS 

Nend=17 

FOR  1=0  TO  Nend 

T( I  )=FNTvsv(Emf( I ) ) 
NEXT  I 

Tvap  =  (T(0)  +  T(  1  )  +  T(2  )  )/3 
Tliq-(T(3)+T(4))/2 
Tliql=T(5) 
IF  lset=0  THEN 

Tsat=Tl iql 
END  IF 
IF  Isat=l  THEN 

Tsat  =Tvap 
END  IF 
IF  Isat-2  THEN 

Tsat=<Tvap+Tliq  )/2 
END  IF 
Jset=Jset+l 

IF  Okaccpt=0  THEN  Jset=Jset-l 
PRINT 

PRINT  USING  "10X,""Data  set  number 
PRINT 
lbeg=0 
Iend=3 

IF  Iact<10  THEN 
Ibeg=Iact 
Iend=Iact 
END  IF 
FOR  I=Ibeg  TO  lend 

Grad=FNGrad<Emf < 1+6 )  ) 

Tout  =  (T< 2*1  +  10  )  +  T( 2*1  +  11  ))/2 

Delt-Toui-T(I+G) 

Tavg=T( 1+6  )+Delt*.5 

IF  I=Iend  THEN 
Delt=T(  10)-T(9) 
Tavg=T(10)-T<9)/2 

END  IF 

Rhoeg=FNRhoeg( Tavg,Egrat ) 

Nueg=FNNueg( Tavg , Egrai  ) 

Mueg=Nueg*Rhoeg 

Cpeg=FNCpeg( Tavg,Egrat ) 

Keg=FNKeg(Tavg,Egrat  ) 

Preg=Cpeg*Mueg/Keg 

Mdot=FNFmcal( I ,T( 1  +  6  ),Pc( I  )  ) 

Ueg=Mdot/(Rhoeg*A>:  ) 

Reeg=Ueg*Di/Nueg 


=    "",DD";Jset 
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2040 

2045 

2050 

2055 

2060 

2065 

2070 

2075 

2080 

2081 

2085 

2086 

2088 

2089 

2090 

2091 

2092 

2094 

2095 

2096 

2097 

2098 

2099 

2100 

2101 

2102 

2104 

2105 

2106 

2110 

2115 

2120 

2125 

2130 

2131 

2132 

2134 

2135 

2136 

21371 

) 

2138! 

2140 

2145 

2146! 

2147 

2148 

2150 

21511 

21521 


Res=4*Mdot/(PI«Mueg*(Di-4*Delta) ) 

Qdot=Mdot*Cpeg*Delt 

Qdpd  )=Qdot/Ao 

IF  1=0  OR  I=Iact  THEN 

Qdpl=Qdp( I  ) 

CALL  Nusselt 
END  IF 

Lmtd=Delt/LOG< (Tsat-T< 1+6  )  >/<Tsat-T< 1+6  >-Delt  )  ) 
Uo(I )=Qdp( I  >/Lmtd 
Cfeg=l 

IF  Reeg<4000  THEN 
IF  lnsert=0  THEN 

Bbl=.0668*<Di/L  )*Reeg*Preg 

Bb2=l  +  .04*( (Di/L  )*Reeg*Preg  )" .6666 

Nueg=3.66+<Bbl/Bb2  ) 
END  IF 

IF  Insert=l  THEN 

Nueg=1.86*<Preg*ReegK . 333333* ( Di/L )" .33333*Cfeg 

Nueg=Cia< Insert, I  >*< 1+5. 484E-3*Preg* . 7*(  Res/Hod >"1 

END  IF 

IF  Insert=2  THEN 

Nueg=Cia( Insert , I  )*<Reeg" .65)*Preg" .46 

END  IF 

IF  Insert=3  THEN 

Nueg=Cia( Insert , I  )*<Reeg".76 )*Preg".46 

END  IF 
ELSE 

BEEP 

PRINT  USING  "10X, ""TURBULENT  CORRELATION 

Nueg=.027*Reeg's  .8*Preg"  .3333*Cfeg 
END  IF 

Hi=Nueg*Keg/Di 

Hod  >  =  l/<  l/Uo(  I  )-Do/(Di*Hi  >-Rm*Al  ) 
IF  1=0  OR  I=Iact  THEN 

PRINT 


"RESULTS  FOR  TUBE  NO . " " , 3D" ; 1+1 


!5)A.5 


PRINT 

PRINT 

USING 

"10X, 

PRINT 

PRINT 

USING 

"10X, 

PRINT 

USING 

"10X, 

PRINT 

USING 

"10X, 

PRINT 

USING 

"10X, 

PRINT 

USING 

"10X, 

PRINT 

USING 

"10X, 

PRINT 

USING 

"10X, 

PRINT 

USIN6 

"10X, 

PRINT 

USING 

"  10X, 

PRINT 

USING 

"10X, 

PRINT 

USING 

"10X, 

"Mass  flow  rate 
"Inside  Tube  Dia 


(m.  ) 


=  "" ,MZ.3DE" ;Maot 

=  "",MZ.3DE";Dia(Itube 


"180  DEG  OUER  Dia.  (HOD) 
"Inlet  temperature 
"Saturation  temp  < Deg 
"DELTA  Tape  Thickness 
"DELT   temp   Dif. 
"Log.  Mean  temp   Dif. 
"Heat  flux 
"Conductivity  E.G. 
"Conductivity  Tube  Metal 


"" ,MZ.3DE";Hod 
=  "",MZ.3DE";T(I+6) 
C)  =  "" ,MZ.3DE" ;Tsat 
=  "\MZ.3DE";  Delta 
=  "" ,MZ.3DE";Delt 
=  "",MZ.3DE";Lmtd 
=  "" ,MZ.3DE" ;Qdp(I  ) 
=  "",MZ.3DE";Keg 

"" ,MZ.3DE" ;Km 
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.MZ 

3DE" 

Rrt 

,MZ 

.3DE" 

Preg 

,MZ 

3DE" 

Reeg 

,mz 

.3DE" 

Res 

,MZ 

3DE" 

Hi 

,MZ 

.3DE" 

Nueg 

,MZ 

3DE" 

Uo(  I  ) 

,MZ 

3DE" 

Kf 

.MZ 

.3DE" 

Rhof 

,MZ 

3DE" 

Hfg 

,MZ 

.3DE" 

Muf 

2153  PRINT  USING  " 10X ," "Resistance  of  Tube  Metal- 

2155  PRINT  USING  " 10X , " "Prandt 1  number 

2160  PRINT  USING  " 10X ," "Reynolds  number 

2161  PRINT  USING  " 10X , " "Reynol ds  number  H&B  S 

2165  PRINT  USING  " 10X , " " Insi de  h.t.c. 

2166  PRINT  USING  " 10X ,"" Inside  NUSSULT  NO. 

2167  PRINT  USING  " 10X , " "OUERALL  H.t.c.  < Uo  ) 
2170  PRINT 

2171!  PRINT  USING  "10X,""KF 

2172!  PRINT  USING  "10X,""RHOF 

2173!  PRINT  USING  "10X,""HFG 

2174!  PRINT  USING  "10X,""MUF 

2175  PRINT 

2177  PRINT  USING  "10X,""TUBE      FM     UEG      DELT     Uo 

HO         HNUS 

2180  PRINT  USING  "10X,""   *      <%>     (m/s)      (K>          (W/wA2. 

K  ) 

2185  PRINT  USING  " 10X ,3D , 4X , 3D .DD , 3X , Z . DD , 4X , DD . 3D , 2X , MZ . 3DE , 3X , MZ . 3DE 

,3X,MZ.3DE,3X" ;I  +  l,Pc( I  ),Ueg,Delt ,Uo( I  ),Ho( I  ),Hnus 

2190  ELSE 

2195  PRINT  USING  " 10X , 3D , 4X , 3D .DD , 3X , Z . DD , 4X , DD . 3D , 2X , MZ . 3DE , 3X ,MZ . 3DE 

,3X";I  +  l,Pc(I  ),Ueg,Delt ,Uo( I  ),Ho< I  ) 

2200  END  IF 

2205  NEXT  I 

2210  IF  ln=0  THEN 

2215  Okaccpt=l 

2220  BEEP 

2225  INPUT  "OK  TO  ACCEPT  THIS  SET  ( 1=DEFAULT=YES ,  0=NO  )? " , Okaccpt 

2230  IF  Okaccpt=l  THEN  OUTPUT  @Fi le ; Pc( *  ) , Emf ( *  ) , Tp< *  ) 

2235  END  IF 

2240  IF  <0kaccpt=l  OR  Irn=l  )  AND  Iout  =  l  THEN 

2245  FOR  1=0  TO  2 

2250  OUTPUT  @Fout;Ho( I  ),Qdp(  I  ) 

2255  NEXT  I 

2260  END  IF 

2265  IF  ln=0  THEN 

2270  Okrpt-1 

2275  BEEP 

2280  INPUT  "UILL  THERE  BE  ANOTHER  RUN  (  1  =  YES=DEFAULT , 0=NO  )" , Okrpt 

2285  IF  Okrpt=l  THEN  1683 

2290  .  ELSE 

2295  IF  Jset<Nsets  THEN  1683 

2300  END  IF 

2305  ASSIGN  SFile  TO  * 

2310  IF  Iout-1  THEN  ASSIGN  ©File  TO  * 

2315  SUBEND 

2320  DEF  FNGrad(T) 

2325  Grad=-3.877857E-5-2*4.7142857E-8*T 

2330  RETURN  Grad 

2335  FNEND 
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2340  DEF  FNKcu(T  ) 

2345  !   OFHC  COPPER  250  TO  300  K 

2350  Tk=T+273.15 

2355  K=434-. 1 12*Tk 

2360  RETURN  K 

2365  FNEND 

2370  DEF  FNNueg(Tc,Egr> 

2375  !   RANGE  OF  VALIDITY:   -20  TO  20  DE6  C 

2380  Tk=Tc+273.15 

2385  Nul=7.U96507E-3-Tk*(7.4863347E-5-Tk*(2.6294943E-7-Tk*3.083332  9E-10)  ) 

2390  Nu2=4.9237638E-3-Tk*(4.9213912E-5-Tk*(1.6437534E-7-Tk*1.8333331E-10)) 

2395  Nu3=8.6586293E-3-Tk*(8.8837902E-5-Tk*(3.0495032E-7-Tk*3.4999996E-10)) 

2400  A2=<Nu3-2*Nu2+Nul >/200 

2405  Al=(Nu2-Nul-940*A2 >/10 

2410  A0=Nul-42*Al-1764*A2 

2415  Nu=A0+Egr*< Al+Egr»A2 > 

2420  RETURN  Nu 

2425  FNEND 

2430  DEF  FNCpeg(Tc,Egr) 

2435  !   RANGE  OF  VALIDITY:  0  TO  20  DEG  C 

2440  Tk=Tc+273.15 

2445  Cpl=1.6701550E+3+Tk*6.3 

2450  Cp2=1.4748125E+3+Tk*6.25 

2455  Cp3=9.5800500E+2+Tk»7.3 

2460  A2=(Cp3-2*Cp2+Cpl)/200 

2465  Al  =  (Cp2-Cpl-900*A2  )/10 

2470  A0=Cpl-40*Al-1600*A2 

2475  Cp=A0+Egr*(Al+Egr*A2) 

2480  RETURN  Cp 

2485  FNEND 

2490  DEF  FNRhoeg(T,Egr  ) 

2495  Rol=1.0607093E+3-T*(3.7031283E-l+T*4.0837183E-3) 

2500  Ro2=1.0748272E+3-T*<4.4266195E-l+T*4.0939706E-3 ) 

2505  Ro3=1.0885934E+3-T*(5.7355653E-l+T*6.1281405E-3) 

2510  A2=(Ro3-2*Ro2+Rol )/200 

2515  Al  =  (Ro2-Rol-900*A2  )/10 

2520  A0=Rol-40*Al-1600*A2 

2525  Ro=A0+Egr*(Al+Egr*A2  ) 

2530  RETURN  Ro 

2535  FNEND 

2540  .   DEF  FNPreg(T,Egr  ) 

2545  Pr=FNCpeg(T,Egr)*FNNueg(T,Egr  )«FNRoeg( T , Egr  )/FNKeg( T , Egr ) 

2550  RETURN  Pr 

2555  FNEND 

2560  DEF  FNKeg(Tc,Egr  ) 

2565  !   RANGE  OF  VALIDITY:   -20  TO  20  DEG  C 

2570  Tk=Tc+273.15 

2575  Kl=2.2824708E-l+Tk*(5.5989286E-4+Tk*3.5714286E-7 ) 

2580  K2=2.5846616E-l+Tk*(2.3978571E-4+Tk«7.1428571E-7 ) 

2585  K3=3.2138932E-l-Tk*(3.0042857E-4-Tk*l .4285714E-6 ) 
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2590 

A2=(K3-2*K2+K1  )/200 

2595 

A1  =  (K2-K1-900*A2  )/10 

2600 

A0=K1-40*A1-1600*A2 

2605 

K=A0+Egr*<Al+Egr*A2 ) 

2610 

RETURN  K 

2615 

FNEND 

2620 

DEF  FNTanh(X ) 

2625 

P=EXP(X ) 

2630 

Q=l/P 

2635 

Tanh=(P-Q)/<P+Q> 

2640 

RETURN  Tanh 

2645 

FNEND 

2650 

DEF  FNTvsv(U) 

2655 

COM  /Cc/  C(7) 

2660 

T=C(0) 

2665 

FOR  1=1  TO  7 

2670 

T  =  T  +  C(  i  >*<ri 

2675 

NEXT  I 

2680 

RETURN  T 

2685 

FNEND 

2690 

DEF  FNBetaU) 

2695 

Rop=FNRho(T+.l  ) 

2700 

Ron=FNRho<T-.l  ) 

2705 

Beta  =  -2/(Rop+Rom  )*(Rop-Rom  )/ .2 

2710 

RETURN  Beta 

2715 

FNEND 

2720 

DEF  FNPsat(Tc) 

2725  i 

0  TO  80  deg  F  CURUE  FIT  OF  Psat 

2730 

Tf=1.8*Tc+32 

2735 

Pa=5.945525+Tf»<  . 15352082+Tf *< 1 . 4840963E-3+Tf *9 

6150671E-6)) 

2740 

Pg=Pa-14.7 

2745 

IF  Pg>0  THEN       !  +=PSIG,-=in  Hg 

2750 

Psat=Pg 

2755 

ELSE 

2760 

Psat=Pg*29.92/14.7 

2765 

END  IF 

2770 

RETURN  Psat 

2775 

FNEND 

2850 

SUB  Xyread 

2855 

BEEP 

2860 

INPUT  "ENTER  FILE  NAME",File$ 

2865 

BEEP 

2870 

INPUT  "ENTER  NUMBER  OF  X , Y  PAIRS" ,N   • 

2875 

ASSIGN  @File  TO  File$ 

2880 

FOR  1=1  TO  N 

2885 

ENTER  @File;X,Y 

2890 

PRINT  X,Y 

2895 

NEXT  I 

2900 

SUBEND 

2905 

SUB  Purge 
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2910 
2915 
2920 
2925 
2930 
2935 
2940 
2945 
2950 
2955 
2960 
2955 
2970 
2975 
2976 
2977 
2980 
2985 
2990 
2995 
3000 
3005 
3010 
3015 
3020 
3025 
3030 
3035 
3040 
3045 
3046 
3050 
3051 
3052 
3053 
3054 
3055 
3056 
3057 
3058 
3059 
3060 
3062 
3063 
3065 
3070 
3075 
3080 
3085 
3090 


BEEP 

INPUT  "ENTER  FILE  NAME  TO  BE  DELETED" , Fi le$ 

PURGE  File$ 

GOTO  2910 

SUBEND 

SUB  Wilson 

COM  /Will/  Doa(4),Dia<4),K(na(4>,Iact,Droot<4) 

COM  /Wil2/  Delta, Isat , Nsets ,Hod,Cia< 3,3 > , Alpaa< 3  ) 

DIM  Enf(  17),T<  17),Xa<  100>,Ya(  100>,Pc(4),Tp<4> 

BEEP 

INPUT  "PLEASE  RE-ENTER  NAME  OF  FILE",File$ 

ASSIGN  ©File  TO  File$ 

INPUT  "0=SMOOTH  COPPER  TUBE,  1=CUNI 

BEEP 

1=  TWISTED  TAPE 


INPUT  "0=NO  INSERT 

BEEP 

INPUT  "GIVE  A  NAME  FOR  XY  FILE",Xy$ 

CREATE  BOAT  Xy$,20 

ASSIGN  @Xy  TO  Xy$ 

L=1.2192 

Do=Doa< Itube) 

Di=Dia< Itube) 

Km=Kma( Itube  ) 

Ax=PI*Di "2/4   !  Cross-sectional  area 

Ao=PI*Do*L 

Rm=Do*LOG(Do/Di  )/<2*Km) 

Initial  values 

Tf=Tsat 

Alpa=.655 

IF  lnsert=0  THEN 

Ci=1.0 

END  IF 

IF  Insert=l  THEN 

Ci=5.172 

END  IF 

IF  Insert=2  THEN 

Ci=.226 

END  IF 

IF  Insert=3  THEN 

Ci=.063 

END  IF 

G=9.81 

lbeg=0 

lend=0   ! CHANGE  TO  4 ,  IF  FIVE  TUBES  IN  BUNDLE 


FIN,  2=K0R0DENSE,  3=FIN  TITAN", Itube 
INSERT,  2=HITRAN1,  3=HITRAN2" , Insert 


IF  Iacta0  THEN 
Ibeg=Iact 
Iend=Iact 

END  IF 
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3095 

FOR  I=Ibeg  TO  lend 

3100 

Sx=0 

3105 

Sy=0 

3110 

5.X5=0 

3115 

Sxy=0 

3120 

Jset-0 

3125 

ASSIGN  ©File  TO  File* 

3130 

ENTER  @Fi le;Dtg$, I  tube, Egrat, Ddl ,Dd2,Dd3 

Dd4,DdS 

3135 

ENTER  eFile;Pc<*),Emf<*>,Tp<#) 

3140 

FOR  J=0  TO  17 

3145 

T< J  )=FNTvsv(Enf( J  )  ) 

3150 

NEXT  J 

3155 

Tvap  =  (T(0)  +  T(  1  )  +  T(2)  )/3 

3160 

Tliq  =  (T(3)  +  T(4)  )/2 

31B1 

Tliql=T(5) 

3165 

IF  lsat«0  THEN 

3170 

Tsat=Tl lql 

3175 

ELSE 

3180 

Tsat=Tvap 

3185 

END  IF 

3190 

Grad=FNGrad(T< 1  +  6  )  ) 

3191 

Tout»<T< 2*1  +  10 )  +  T( 2*1  +  11  ))/2 

3195 

Delt»Tout-T(I+6) 

3200 

Tavg=T(I+6  )+Delt*.5 

3201 

IF  I=Iend  THEN 

3202 

Delt-T(10)-T<9) 

3204 

Tavg=(T<10)-T(9))/2 

3205 

END  IF 

3206 

3210 

Water/Et hy lene  Glycol  Mixture  Froperties 

3215 

Rhoeg=FNRhoeg( Tavg , Egrat  ) 

3220 

Nueg=FNNueg( Tavg , Egrat  ) 

3225 

Mueg=Nueg*Rhoeg 

3230 

Cpeg=FNCpeg( Tavg, Egrat  ) 

3235 

Keg=FNKeg( Tavg, Egrat  ) 

3240 

Preg=Cpeg*Mueg/Keg 

3245 

l 

3250 

Mdot=FNFncal( I , T( 1+6  ) , Pc ( I  )  ) 

3255 

Veg=Mdot/<Rhoeg*Ax  ) 

3260 

Reeg=Veg*Di /Nueg 

3265 

Re5  =  4*Ndot/<Mueg*(PI*Di-4*Delta  )  ) 

3270 

Qdot=Mdot*Cpeg*Delt 

3275 

Qdp=Qdot/Ao 

3280 

Lmtd=Delt/LOG(  (Tsat-KI+B)  )/<Tsat-T<  1  +  6  >- 

-Delt  )  ) 

3285 

Uo=Qdp/Lmtd 

3286 

IF  lnsert=0  THEN 

3287 

Bbl=.0668*(Di/L  )*Reeg*Preg 

3288 

Bb2  =  l  +  .04*(  <Di/L)*Reeg*Preg)A.6666 

3289 

0Piega  =  3.B6+(Bbl/Bb2  ) 

3290 

END  IF 
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3291 
3293 
3294 
3295 
329G 
3297 
3298 
3299 
3300 
3302 
3303 
3305 
3310 
3315 
3320 
3325 
3330 
3335 
3340 
3345 
3350 
3355 
3360 
33G5 
3370 
3375 
3380 
3385 
3390 
3395 
3400 
3405 
3410 
3415 
3420 
3425 
3430 
3435 
3440 
3445 
3450 
3455 
3460 
3465 
3470 
3475 
3480 
3485 
3490 
3495 


IF  Insert=l  THEN 

Omega=( 1+5. 484E-3*PregA . 7*< Res/Hod ) A 1 .25 )" . 5 

END  IF 

IF  Insert=2  THEN 

Omega=(ReegA .65>*<Preg" .46 ) 

END  IF 

IF  Insert=3  THEN 

Onega  =  ( ReegA . 76  )*( Preg" . 46  ) 

END  IF 

R-114/R113  Properties 
Hfg=FNHfg(Tsat  ) 
Kf=FNK(Tf ) 
Rhof=FNRho<Tf ) 
Muf=FNMu<Tf > 

F=(KfA3*RhofA2*G*Hfg/(Muf*Do«Qdp)  >A.  33333 

Ho=!°ilpa*F 

Two-Tsat-Qdp/Ho 

Tf=Tsat/3+2*Two/3 

Y=( l/Uo-Rm)»F 

X=Do*F/(Keg*Omega ) 

PRINT  "OMEGA=" ; Omega ; "F=" ;F; "X=" ;X; "Y=" ;Y 

Xa<Jset)  =  X    !   INEFFICIENT  (MODIFY  LATER) 

Ya( Jset  )=Y 

5x=Sx+X 

Sy=Sy+Y 

Sxs=Sxs+X*X 

Sxy=Sxy+X*Y 

Jset=Jset+l 

IF  JseUNsets  THEN  3135 

ASSIGN  ©File  TO  * 

Slope=(Nsets*Sxy-Sx*Sy  )/ ( Nset s*Sxs-Sx A2 ) 

Intcpt=( Sy-Slope*Sx )/Nsets 

Cic»l/Slope 

Alpac=l/Intcpt 

Cerr=ABS( (Ci-Cic )/Cic ) 

Aerr=ABS( < Alpac-Alpa )/Alpac ) 

IF  Cerr>.001  OR  Aerr>.001  THEN 

Alpa=(Alpa+Alpac  >*.5 

Ci  =  (Ci+Cic  )*.5 

PRINT  "CIC=" ;Cic; "ALPA=" ; AlPa 

GOTO  3100 
END  IF 
BEEP 
BEEP 

PRINTER  IS  1 

PRINT  "CIC=";Cic; "ALPA=" ;Alpa 
Cia( Insert , I  )=Cic 
Alpaa< I  )  =  Alpac 
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3500        PRINTER  IS  701 

3505        FOR  J=0  TO  Nsets-1 

3510  OUTPUT  @Xy;Xa< J  ),Ya< J  ) 

3515        NEXT  J 

3520        PRINTER  IS  1 

3525  NEXT  I 

3530  ASSIGN  @Xy  TO  * 

3535  SUBEND 

3540  SUB  Nusselt 

3545  COM  /Nus/  I  in , Tsat ,Qdp , Hoc ,K f , Rhof , Hf g , Muf , Do , 1 1 ube 

3550  i   Do=.0159 

3555  Ho=1000 

35B0  IF  lin=0  THEN 

35G5       BEEP 

3570       INPUT  "ENTER  TSAT  AND  HEAT  FLUX ", Tsat , Qdp 

3575  END  IF 

3580  Hfg=FNHfg(Tsat  ) 

3585  Two=Tsat-Qdp/Ho 

3590  Tf=Tsat/3+2*Two/3 

3595  Kf=FNK(Tf) 

3G00  Rhof=FNRho(Tf  ) 

3G05  Muf=FNMu(Tf) 

36 10  Hoc=.655*(Kf  A3*Rhof'v2*9.81*Hfg/<Muf*Do*Qdp  )  )A.  333333 

3B15  IF  ABS((Ho-Hoc  )/Hoc  )>.001  THEN 

3G20       Ho=(Ho+Hoc )*.5 

3625       GOTO  3585 

3G30  END  IF 

3635  IF  lin=0  THEN  PRINT  "HO=";Hoc 

3640  SUBEND 

3645  SUB  Plotl 

3650  DIM  Yaa(4  ) 

3655  PRINTER  IS  705 

3660  Idv=i 

3665  BEEP 

3670  INPUT  "OK  TO  USE  DEFAULT  UALUES  (  1=DEF  =  Y , 0  =  N ) " , Idv 

3675  IF  Idv-1  THEN 

3680       Itn-2 

3685       Xnin=l 

3G90       Xmax=5 

3695       Xstep=l 

3700       Ymin-0 

3705       Ymax=2.0 

3710       Ystep=.5 

3715  ELSE 

3720       BEEP 

3725       INPUT  "ENTER  MINIMUM  AND  MAXIMUM  X-UALUES" , Xmin , Xmax 

3730       BEEP 

3735       INPUT  "ENTER  MINIMUM  AND  MAXIMUM  Y-UALUES"  ,  Yrcin  ,  Ymax 

3740       BEEP 

3745       INPUT  "ENTER  STEP  SIZE  FOR  X-AXIS ",  Xstep 
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3750       BEEP 

3755       INPUT  "ENTER  STEP  SIZE  FOR  Y-AXIS" , Yst ep 

3760  END  IF 

3765  PRINT  "IN;SP1;IP  2300,1800,8300,6800;" 

3770  PRINT  "SC  0 , 100 , 0 , 100 ; TL  2,0:" 

3775  Sfx=100/<Xmax-Xmin ) 

3780  Sfy=100/( Ymax-Ymin) 

3785  PRINT  "PU  0,0  PD" 

3790  FOR  Xa-Xmin  TO  Xmax  STEP  Xstep 

3795        X  =  <Xa-Xmn)*Sfx 

3800        PRINT  "PA";X," ,0;  XT; " 

3805  NEXT  Xa 

3810  PRINT  "PA  100,0;PU; " 

3815  PRINT  "PU  PA  0,0  PD" 

3820  FOR  Ya  =  Ymm  TO  Ymax  STEP  Ystep 

3825        Y=(Ya-Ymin)*Sfy 

3830        PRINT  "PA  0, ";Y, "YT" 

3835  NEXT  Ya 

3840  PRINT  "PA  0,100  TL  0  2" 

3845  FOR  Xa=Xmin  TO  Xmax  STEP  Xstep 

3850        X=<Xa-Xmin)*Sfx 

3855        PRINT  "PA" ; X , " , 100 ;  XT" 

3860  NEXT  Xa 

3865  PRINT  "PA  100,100  PU  PA  100,0  PD" 

3870  FOR  Ya=Ymin  TO  Ymax  STEP  Ystep 

3875        Y=<Ya-Ymin)»Sfy 

3880        PRINT  "PD  PA  100,",Y,"YT" 

3885  NEXT  Ya 

3890  PRINT  "PA  100,100  PU" 

3895  PRINT  "PA  0,-2  SR  1.5,2" 

3900  FOR  Xa=Xmin  TO  Xmax  STEP  Xstep 

3905        X=< Xa-Xmin >*Sfx 

3910         PRINT  "PA";X,",0; " 

3915        PRINT  "CP  -2,-l;LB";Xa; "" 

3920  NEXT  Xa 

3925  PRINT  "PU  PA  0,0" 

3930  FOR  Ya=Ymin  TO  Ymax  STEP  Ystep 

3935        Y=< Ya-Ymin>*Sfy 

3940         PRINT  "PA  0 , " ; Y , "" 

3945        PRINT  "CP  -4 , - . 25; LB" ; Ya ; " " 

3950  .   NEXT  Ya 

3955  BEEP 

3960  INPUT  "SELECT  MODE  ( 0=HN/H1 , 1=HN( avg  )/Hl )" , I sm 

3965  Ism=Ism+l 

3970  IF  Idv=l  THEN 

3975       IF  Ism=l  THEN  Ylabel$="HN/Hl" 

3980       IF  Ism=2  THEN  Ylabel$="HN( avg  ) /HI " 

3985       Xlabel$="Tube  Number" 

3990  ELSE 

3995       BEEP 
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4000  INPUT  "ENTER  X-LABEL " , Xlabel$ 

4005  BEEP 

4010  INPUT  "ENTER  Y-LABEL" , Ylabel$ 

4015  END  IF 

4020  PRINT  "SR  1.5,2;PU  PA  50,-10  CP" ; -LEN( Xlabel$  )/2 ; "0 ; LB " ; Xlabel$; " " 

4025  PRINT  "PA  -11,50  CP  0 , " ; -LEN< Ylabel$  )/2*5/G ; "DI  0 , 1 ; LB" ; Ylabel $ • " " 

4030  PRINT  "CP  0,0  DI" 

4035  Okp=l 

4040  BEEP 

4045  INPUT  "WANT  TO  PLOT  DATA  FROM  A  FILE  ( 1=DEF=Y , 0=N  )?" , Okp 

4050  IF  Okp=l  THEN 

4055  BEEP 

40B0  INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE" , Df l le$ 

4065  ASSI6N  @File  TO  Dfile$ 

4070  BEEP 

4075  INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER", Md 

4080  BEEP 

4085  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED" , Nsets 

4090  BEEP 

4095  INPUT  "SELECT  A  SYMBOL  FOR  THE  PLOTTER  < 1  =  * , 2=+  ,  3-c , 4=o , 5="  )" , Sy 

4100  PRINT  "PU  DI" 

4105  IF  Sy=l  THEN  PRINT  "SM*" 

4110  IF  Sy=2  THEN  PRINT  "SM+" 

4115  IF  Sy=3  THEN  PRINT  "SMc" 

4120  IF  Sy=4  THEN  PRINT  "SMo" 

4125  IF  Sy=5  THEN  PRINT  "SM"" 

4130  FOR  1=1  TO  Nsets 

4135  FOR  J=0  TO  3 

4140  ENTER  ©File;Yaa< J  ),D 

4145  IF  J=0  THEN  Ytop=Yaa<0) 

4150  Yaa<  J)  =  Yaa<  J>/Ytop 

4155  NEXT  J 

4160  FOR  J=0  TO  3 

41G5  X  =  (  J+l-Xmin)*Sfx 

4170  Y=(Yaa( J >-Ymin >*Sf y 

4175  PRINT  "PA",X,Y,"" 

4180  NEXT  J 

4185  NEXT  I 

4190  BEEP 

4195  ASSIGN  GFile  TO  * 

4200  GOTO  4040 

4205  END  IF 

4210  BEEP 

4215  INPUT  "LIKE  TO  PLOT  THE  NUSSELT  RELATION  <  1=Y  ,  0  =  N  )? "  ,  Oknus 

4220  PRINT  "PU;SM" 

4225  IF  Oknus=l  THEN 

4230  FOR  Xa=Xmin  TO  Xmax  STEP  Xstep/50 

4235  X-<Xa-Xmin)*Sfx 

4240  IF  Ism=l  AND  Xa>Xmin  THEN  Ya-XaA .75-< Xa-1  )" .75 

4245  IF  Ism=2  AND  Xa>Xmin  THEN  Ya=XaA<-.25^ 
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4250  IF  Xa-Xmin  THEN  Ya=l 

4255  Y=< Ya-Ymin)*Sfy 

4260  PRINT  "PA" ,X,Y,"PD" 

4265  NEXT  Xa 

4270  BEEP 

4275  PRINT  "PU" 

4280  INPUT  "MOVE  THE  PEN  TO  LABEL  THE  NUSSELT  LINE", Ok 

4285  PRINT  "LBNusselt " 

4290  END  IF 

4295  IF  Isn=2  THEN 

4300  BEEP 

4305  INPUT  "LIKE  TO  PLOT  EXPTL  CURVE  < 1  =  Y ,0=N  ) " ,Okex 

4310  Nq=0 

4315  IF  Okex-1  THEN 

4320  BEEP 

4325  INPUT  "ENTER  THE  EXPONENT", Ex 

4330  FOR  Xa=Xmin  TO  Xnax  STEP  Xstep/10 

4335  Nq=Nq+l 

4340  Ya  =  Xa/,(-Ex) 

4345  X-<Xa-Xmin)*Sfx 

4350  Y=(Ya-Ymin)*Sfy 

4355  IF  Nq  MOD  2=0  THEN 

4360  PRINT  "PA" ,X,Y, "PD" 

4365  ELSE 

4370  PRINT  "PA",X,Y, "PU" 

4375  END  IF 

4380  NEXT  Xa 

4385  PRINT  "PU" 

4390  BEEP 

4395  INPUT  "MOVE  PEN  TO  LABEL  AND  HIT  ENTER", Ok 

4400  PRINT  "LBs=0" 

4405  PRINT  "PR  -1  0" 

4410  PRINT  "LB";Ex; "" 

4415  GOTO  4300 

4420  END  IF 

4425  END  IF 

4430  GOTO  4510 

4435  BEEP 

4440  INPUT  "LIKE  TO  PLOT  KERN  RELATIONSHIP  ( 1=Y ,0»N  )?" , Yes 

4445  IF  Yes=l  THEN 

4450  FOR  Xa=Xmin  TO  Xmax  STEP  Xstep/20 

4455  Ya=Xa/,(-l/6  ) 

4460  X  =  (Xa-Xrnin)*Sfx 

4465  Y=(Ya-Ynin)*Sfy 

4470  PRINT  "PA" ,X,Y, "PD" 

4475  NEXT  Xa 

4480  PRINT  "PU" 

4485  BEEP 

4490  INPUT  "MOVE  THE  PEN  TO  LABEL  KERN  RELATIONSHIP" ,Ok 

4495  PRINT  "LBKern;PU" 
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4500  END  IF 

4505  PRINT  "PU  PA  0,0" 

4510  BEEP 

4515  INPUT  "LIKE  TO  PLOT  EISSENBERG  RELATION  ( 1  =  Y ,0=N  )? " , Okei 

4520  IF  Okei=l  THEN 

4525       FOR  Xa=Xmin  TO  Xnax  STEP  Xstep/10 

4530  Ya=.6+.42*Xa"<-.25) 

4535  X  =  (Xa-Xnin  )»Sfx 

4540  Y=(Ya-Ymin>*Sfy 

4545  PRINT  "PA" , X , Y , "PD" 

4550       NEXT  Xa 

4555       PRINT  "PU" 

4560       BEEP 

4565       INPUT  "MOUE  THE  PEN  TO  LABEL  THE  EIS5ENBERG  LINE" ,Ok 

4570       PRINT  "LBEissenberg;PU" 

4575  END  IF 

4580  PRINT  "PU  SP0" 

4585  SUBEND 

4590  DEF  FNPvst(Tc  ) 

4595  COM  /Fid/  Ift 

4600  DIM  K<8) 

4605  IF  lfi=0  THEN 

4610       BEEP 

4615       PRINT  "PUST  CORRELATION  NOT  AVAILABLE  FOR  R-114" 

4620       STOP 

4625  END  IF 

4630  IF  Ifi=l  THEN 

4635       DATA  -7.691234564,-26,08023696,-168.1706546,64.23285504,-118.9646225 

4640       DATA  4 .  1671 1732 , 20 . 9750676 , 1 .E9, 6 

4645       READ  K<* ) 

4650       T=(Tc+273.15)/647.3 

4655       Sum=0 

4660       FOR  N=0  TO  4 

4665  Sum=Sum+K(N)*< 1-T  )"(N+1  ) 

4670       NEXT  N 

4675       Br  =  Sum/(T*(  1+K(  5  )*  (  1-T  )+K  (  6  >*<  1-T  T2  )  )-<  1-T  ) /<  K  (  7  >*  (  1-T  )  "2+K(  8  )  ) 

4680       Pr=EXP(Br) 

4685       P=22120000*Pr 

4690  END  IF 

4695  IF  Ift-2  THEN 

4700       Tf=Tc*l. 8+32+459. 6 

4705       P=10'(33. 0655-4330. 98/Tf-9 . 2635*L6T( Tf )+2 . 0539E-3*Tf ) 

4710       P=P*101325/14.696 

4715  END  IF 

4720  IF  IH=3  THEN 

4725       A=9.394685-3066.1/(Tc+273.15> 

4730       P=133.32*10AA 

4735  END  IF 

4740  RETURN  P 

4745  FNEND 


197 


4750  DEF  FNHfg(T) 

4755  COM  /Fid/  I  ft 

4760  IF  If t=0  THEN 

4765       Tf=T*1.8+32 

4770       Hfg=6.1451558E+l-Tf*(6.951079E-2+Tf*( 1 . 3988688E-4+1 . 9607843E-7*Tf ) ) 

4775       Hfg=Hfg*2326 

4780  END  IF 

4785  IF  I f t=l  THEN 

4790       Hfg=2477200-2450*<T-10> 

4795  END  IF 

4800  IF  Ift-2  THEN 

4805       Tf=T*1.8+32 

4810       Hfg=7.0557857E+l-Tf*(4.838052E-2+1.2619048E-4*Tf ) 

4815       Hfg=Hfg*2326.1 

48161       Hfg=( 1.611-.0031*T  )»1.0E+5 

4820  END  IF 

4825  IF  Ift=3  THEN 

4830       Tk=T+273.15 

4835       Hfg=1.35264E+6-Tk»(6.38263E  +  2  +  Tk*. 747462  ) 

4840  END  IF 

4845  RETURN  Hfg 

4850  FNEND 

4855  DEF  FNMu(T) 

4860  COM  /Fid/  I  ft 

4865  IF  If t=0  THEN 

4870       Tk=T+273.15 

4875       Mu=EXP(-4.4636  +  1011.47/Tk  >*l.E-3 

4880  END  IF 

4885  IF  I f t =1  THEN 

4890       A=247.8/(T+133.15) 

4895       Mu=2.4E-5*10"A 

4900  END  IF 

4905  IF  I f t=2  THEN 

4910  i      Mu=8.9629819E-4-T*( 1 . 1094609E-5-T*5 .566829E-8 ) 

4911  Tk=T+273.15 

4913       Mu=1.34E-5*10.0/,<503.0/(Tk-2.15)  ) 

4915  END  IF 

4920  IF  I f t=3  THEN 

4925       Tk=l/(T+273.15) 

4930       Mu=EXP(-11.0179+Tk*<  1  .  744E+3-Tk*<  2  .  80335E+5-TM1 .  1266 1E+8  )  )  ) 

4935  END  IF 

4940  RETURN  Mu 

4945  FNEND 

4950  DEF  FNVvst(Tt ) 

4955  COM  /Fid/  I  ft 

4960  IF  lft=0  THEN 

4965       BEEP 

4970  PRINT    "UUST   CORRELATION   NOT    MAILABLE   FOR   R-114" 

4975  STOP 

4980  END    IF 
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4985  IF  Ift-1  THEN 

4990       P=FNPvst(Tt) 

4995       T=Tt+273.15 

5000       X=1500/T 

5005       F 1  =  1 / ( l  +  T*l.E-4  ) 

5010       F2  =  ( 1-EXP(-X  ))A2.5*EXP(X  )/XA.5 

5015       B=. 00 15*F1-. 000942 *F2-. 0004882 *X 

5020       K=2*P/(4G1.52*T ) 

5025       y»(l+<l+2*B*K)*.5)/K 

5030  END  IF 

5035  IF  If -t  =2  THEN 

5040       Tf=Tt*1.8+32 

5045       V=13.955357-Tf*(  .  16  127262-Tf *5 . 1726190E-4 ) 

5050       V=V/1B.018 

5055  END  IF 

50B0  IF  Ift=3  THEN 

50G5       Tk=Tt+273.15 

5070       P=FNPvst(Tt) 

5075       V-133.95*Tk/P 

5080  END  IF 

5085  RETURN  V 

5090  FNEND 

5095  DEF  FNCp(T) 

5100  COM  /Fid/  Ift 

5105  IF  Ift-0  THEN 

5110       Tk=T+273.15 

5115       Cp=.40118+Tk*( 1.65007E-3+Tk«(1.51494E-6-Tk*6.67853E-10) ) 

5120  END  IF 

5125  IF  I f i  =  1  THEN 

5130       Cp  =  4.2U20858-T*(2.26826E-3-T*<4.42361E-5+2.71428E-7*T}> 

5135  END  IF 

5140  IF  Ift=2  THEN 

5145!       Cp=9.2507273E-l+T*(9.3400433E-4+1.7207792E-6*T) 

5146       Cp=(929.0+1.03*T)/1000. 

5150  END  IF 

5155  IF  I  f t=3  THEN 

51B0       Tk=T+273.15 

5165       Cp=4.1868*(1.6884E-2+Tk*(3.35083E-3-Tk*(7.224E-6-Tk*7.61748E-9)) ) 

5170  END  IF 

5175  RETURN  Cp*1000 

5180  FNEND 

5185  DEF  FNP.ho(T) 

5190  COn  /Fid/  Ift 

5195  IF  If t=0  THEN 

5200       Tk=T+273.15 

5205       X=l-(1.8*Tk/753.95> 

5210       Ro=36.32+61.146414*X'M  1/3  )+ 16 . 418015*X  +  17 .476838**" .5+1 . 119828*XA2 

5215       Ro=Ro/. 062428 

5220  END  IF 

5225  IP  Ift=l  THEN 
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5230       Ro-999.S2946+T*<  . 012G9-T*( 5 . 482513E-3-T* 1 .234147E-5  )  ) 

5235  END  IF 

5240  IF  Ift-2  THEN 

5245  Ro=1.6207479E+3-T*(2.2186346+T*2.3578291E-3) 

5246  1        Ro-l./< <  ,G17+.000E47*T"1 .1  )»1.0E-3) 
5250  END  IF 

5255  IF  If t=3  THEN 

5260       Tk-T+273. 15-338. 15 

5265       Uf=9.24848E-4+Tk*(6.2796E-7+Tk*<  9  .  2444E-10+Tk*3  .057E-12  ) 

5270       Ro-l/Vf 

5275  END  IF 

5280  RETURN  Ro 

5285  FNEND 

5290  DEF  FNPr(T  ) 

5295  Pr=FNCp(T)*FNMu(T  )/FNK(T  ) 

5300  RETURN  Pr 

5305  FNEND 

5310  DEF  FNK(T) 

5315  COM  /Fid/  Ift 

5320  IF  Ift-0  THEN  K= . 071- . 000261 *T 

5325  IF  Ift-1  THEN 

5330       X=(T+273.15)/273.15 

5335       K=-.92247+X*(2.8395-X*< 1 .8007-X*(  .52577- .07344*X  )  )  ) 

5340  END  IF 

5345  IF  Ift-2  THEN 

5350!      K  =  8.2095238E-2-T*(2.2214286E-4+T*2.3809524E-8  ) 

5351!       K=.0802-.000203*T 

53521      K=.071-.000261*T 

5353  Tf=T*1.8+32 

5354  K=4.8461905E-2-Tf*6.5714286E-5 

5355  K=K*1.7308 
5357  END  IF 

5360  IF  Ift-3  THEN 

5365       Tk=T+273.15 

5370       K  =  4.1868E-4*(519.442+.320920*Tk  ) 

5375  END  IF 

5380  RETURN  K 

5385  FNEND 

5390  DEF  FNHf(T) 

5395  COM  /Fid/  Ift 

5400  IF  Ift-0  THEN 

5405       BEEP 

5410       PRINT  MHF  CORRELATION  NOT  FOR  R-114" 

5415       STOP 

5420  END  IF 

5425  IF  Ift=l  THEN 

5430       Hf=T*(4.203849-T*(5.88132E-4-T*4.55160317E-6) ) 

5435  END  IF 

5440  IF  Ift-2  THEN 

5445       Tf-T*1.8+32 
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5450 
5455 
54B0 
5465 
5470 
5475 
5480 
5485 
5490 
5495 
5500 
5505 
5510 
5515 
5520 
5525 
5530 
5535 
5540 
5545 
5550 
5555 
5560 
5565 
5570 
5575 
5580 
5585 
5590 
5595 
5600 
5605 
5610 
5615 
5620 
5625 
5630 
5635 
5640 
5645 
5650 
5655 
5560 
5665 
5570 
5675 
5680 
5685 
5690 
5695 


"Select  Option  X-Y  Limits 

"0  Use  default  values 

"1  Use  new  values 


CODE" , Icode 


Hf =8. 2078571 +Tf*<  . 19467857+Tf *1 . 3214286E-4 > 
Hf=Hf*2.326 
END  IF 
IF  Ift-3  THEN 

Hf=250   !  TO  6E  VERIFIED 
END  IF 

RETURN  Hf*1000 
FNEND 
SUB  PlotZ 

COM  /Drl/  Star, Sym, Icon 
COM  /Fid/  Ift 
DIM  C(9),Xya(7),Doa(4) 
DATA  0.0158,0.0158,0.0158,0.0158 
READ  Doa<*  ) 
Fw»l 

PRINTER  IS  1 
BEEP 

PRINT  USING  "4X, " 
PRINT  USING  "6X," 
PRINT  USING  "6X, " 
INPUT  Okd 
BEEP 

INPUT  "ENTER  TUBE 
Do=Doa< Itube) 
Iht=2 
BEEP 

PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
INPUT  Iht 
PRINTER  IS  705 
IF  Okd=0  THEN  IAXIS 
IF  lht=0  THEN  !<h 

Ymin=0 

Ynax=60 

Ystep=10 

Xmin= .2 

Xmax=1.4 

Xstep=.2 
END  IF 
IF  Iht-1  THEN 

Xnin=0 

Ymin=0 

Ymax=.5 

Xmax«15 

Xstep=3 

Ystep=. 1 
END  IF 
IF  Iht=2  THEN  !(h  vs  t  ) 


"4X, 
"6X, 
"6X, 
"6X, 


Select  opt  ion; ' 


versus  q 

versus  Delta-T 

versus  Delta-T  (default) 


DEFAULT  VALUES 
vs  q  ) 


!  (q  vs  t  ) 
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ENTER  MINIMUM  AND  MAXIMUM  X-UALUES" , Xmin ,Xmax 
ENTER  MINIMUM  AND  MAXIMUM  Y-VALUES" , Ymin , Ymax 
ENTER  STEP  SIZE  FOR  X-AXIS" , Xstep 
ENTER  STEP  SIZE  FOR  Y-AXIS" , Ystep 


5700  Xmin=0 

5705  Ymin=0 

5710  Xmax=50 

5715  Ymax=6 

5720  Xstep=10 

5725  Ystep=l 

5730       END  IF 

5735  END  IF 

5740  IF  Okd=l  THEN 

5745       BEEP 

5750       INPUT 

5755       BEEP 

5760       INPUT 

5765       BEEP 

5770       INPUT 

5775       BEEP 

5780       INPUT 

5785  END  IF 

5790  BEEP 

5795  PRINT  "IN;SP1;IP  2300,1800,8300,6800;" 

5800  PRINT  "SC  0, 100,0, 100;TL  2,0;" 

5805  Sfx=100/(Xmax-Xmin> 

5810  Sfy=100/(Ynax-Ymin> 

5815  BEEP 

5820  lcg=0 

5825  INPUT  "LIKE  TO  BY-PASS  CAGE  ( 1=Y , 0  =  N=DEFAULT  )? " , leg 

5830  IF  Icg=l  THEN  6175 

5835  PRINT  "PU  0,0  PD" 

5840  FOR  Xa=Xmin  TO  Xmax  STEP' Xstep 

5845        X  =  (Xa-Xnin)*Sfx. 

5850        PRINT  "PA";X,U ,0;  XT; " 

5855  NEXT  Xa 

5860  PRINT  "PA  100,0;PU; " 

5865  PRINT  "PU  PA  0,0  PD" 

5870  FOR  Ya=Ymin  TO  Ymax  STEP  Ystep 

5875        Y=(Ya-Ymin)*Sfy 

5880         PRINT  "PA  0, ";Y, "YT" 

5885  NEXT  Ya 

5890  PRINT  "PA  0,100  TL  0  2" 

5895  FOR  Xa=Xmin  TO  Xmax  STEP  Xstep 

5900  .       X=<Xa-Xmin>*Sfx 

5905        PRINT  "PA";X," ,100;  XT" 

5910  NEXT  Xa 

5915  PRINT  "PA  100,100  PU  PA  100,0  PD" 

5920  FOR  Ya=Ymin  TO  Ynax  STEP  Ystep 

5925        Y  =  <Ya.-Ynin)*Sfy 

5930        PRINT  "PD  PA  100,",Y,"YT" 

5935  NEXT  Ya 

5940  PRINT  "PA  100,100  PU" 

5945  PRINT  "PA  0,-2  SR  1,5,2" 
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5950  FOR  Xa=Xmin  TO  Xmax  STEP  Xstep 

5955  X=(Xa-Xmin)*Sfx 

5960  PRINT  "PA";X," ,0; " 

59B5  IF  Xa<l  AND  Xa<>0  THEN  PRINT  "CP  -1 . 5, -1 ; LB0; PR  -1 , 0; LB " ; Xa ; " " 

5970  IF  Xa=0  THEN  PRINT  "CP  -.5,-1 ;LB0"  . 

5975  Xin=0 

5980  IF  Xa  MOD  1=0  THEN  Xin=l 

5985  IF  Xa>=10  THEN  PRINT  "CP  -2 ,-1 ; LB" ; Xa ; " " 

5990  IF  Xa>l  AND  Xa<10  AND  Xin=l  THEN  PRINT  "CP  -1 . 25 , -1 ; LB" ; Xa ; " " 

5995  IF  Xa>l  AND  Xin=0  THEN  PRINT  "CP  -2 , -1 ; LB" ; Xa ;  " " 

G000  IF  Xa=l  THEN  PRINT  "CP  -1,-1; LB1 .0" 

6005  NEXT  Xa 

B010  PRINT  "PU  PA  0,0" 

G015  FOR  Ya=Ymin  TO  Ymax  STEP  Ystep 

B020  Y=(Ya-Yrun)*Sfy 

G025  PRINT  "PA  0,";Y,"" 

6030  IF  lht=0  AND  Ya>0  THEN  PRINT  "PR  2,0" 

G035  IF  Ya<l  AND  Ya<>0  THEN  PRINT  "CP  -4 , - . 25 ; LB0 ; PR  -2 , 0 ; LB" ; Ya; " " 

G040  IF  Ya=0  THEN  PRINT  "CP  -2,-.25;LB0" 

G045  IF  Ya>l  AND  IhU2  THEN  PRINT  "CP  -5  ,  -  .  25  ;  LB"  ;  Ya  ;  " " 

6050  IF  Ya>9  AND  Iht=2  THEN  PRINT  "CP  -4 , - .25 ; LB" ; Ya ; " " 

6055  IF  Ya>0  AND  Ya<10  AND  Iht=2  THEN  PRINT  "CP  -3 , - . 25 ; LB" ; Ya ; " " 

G0G0  i  IF  Ya=l  THEN  PRINT  "CP  -4 ,- . 25 ; LB1 . 0 " 

G0G5  NEXT  Ya 

6070  IF  Okd=l  THEN 

G075  BEEP 

6080  INPUT  "ENTER  X-LABEL" , Xlabel$ 

6085  BEEP 

6090  INPUT  "ENTER  Y-LABEL" , Ylabel$ 

6095  END  IF 

G100  IF  IhtOl  THEN 

6105  PRINT  "SR  1.5,2;PU  PA  -12,35;DI  0,l;LBh;PR  1 , 0 . 5 ; LBo ; PR  -l,0.5;LB/< 

kW/rn1' 

6110  PRINT  "PR  -1,0.5;SR  1,1.5;LB2;SR  1.5,2;PR  .5,.5;LB.;PR  .5,0;LBK>" 

6115  ELSE 

6120  PRINT  "PA  -12,39;DI  0 , 1 ; LBq/( MU/m ; SR  1,1.5;PP  -1 , 0 . 5; LB2 ; SR  1.5,2;P 

R  1.0, 0.5; LB)" 

6125  END  IF 

6130  IF  lht=0  THEN 

6135  PRINT  "SR  1.5,2;PU  PA  40 , -10 ; DI ; LBq/( MW/m; PR  0 . 5 , 1 ; SR  '  1 , 1 . 5 ; LB2 ; SR 

1.5, 2; PR  .5,-1 ; LB)" 

6140  ELSE 

6145  PRINT  "DI  PA  38 , -10 ; LB( T ; PR  .5,-l;LBs;PR  .5,1;LB-T;PR  -2.4,3  PD  PR 

2,0  PU;PR  .5,-4;" 

6150  PRINT  "LBwo;PR  .5,1;LB)/K" 

G155  END  IF 

6160  PRINT  "CP  0,0  DI" 

6165  Xlg-l.E+G 

6170  Xug=-l.E+6 

6175  Xal=50 
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6180 

Yal=95 

6185 

Nrun=0 

6190 

BEEP 

6195 

INPUT  "WANT  TO  PLOT  DATA  FROM  A  FILE  <  1-Y,0=»N )?"  ,Ok 

6200 

XI 1=1 .E+G 

6205 

Xul  — l.E  +  6 

6210 

Okp=0 

6215 

IF  Ok=l  THEN 

6220 

BEEP 

6225 

INPUT  "ENTER  THE  NAME  OF  THE  PLOT  DATA  FILE" ,D_f l le$ 

6230 

ASSIGN  GFile  TO  D_file$ 

6235 

IF  IcombOO  THEN  6265 

6240 

Sx=0 

6245 

Sy=0 

6250 

Sx2=0 

6255 

Sxy=0 

6260 

Md=l 

6265 

BEEP 

6270 

INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER  <DEF=l)",Md 

6275 

Npairs=9 

6280 

BEEP 

6285 

INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED  (DEF  =  9)",Np( 

airs 

6290 

Nrun=Nrun+Npairs 

6295 

PRINTER  IS  1 

6300 

BEEP 

6305 

PRINT  USING  "4X," "Select  a  symbol: 

6310 

PRINT  USING  "4X,""  1   Star    2   Plus  sign 

6315 

PRINT  USING  "4X,""  3   Circle   4   Square 

6320 

PRINT  USING  "4X,""  5   Rorcbus 

6325 

PRINT  USING  "4X,H"  6   Right-side-up  triangle 

6330 

PRINT  USING  "4X,""  7   Up- si  de-down  triangle""" 

6335 

INPUT  Sym 

6340 

BEEP 

6345 

INPUT  "ENTER  TUBE  NUMBER  FOR  PLOTTING  < 0=TOP , 1=SEC0ND , . 

.  .  )" 

6350 

PRINTER  IS  705 

6355 

IF  Sym=l  THEN  PRINT  "SM*" 

6360 

IF  Syrc  =  2  THEN  PRINT  "SM+" 

6365 

IF  Sym=3  THEN  PRINT  "SMo" 

6370 

IF  Md>l  THEN 

6375 

FOR  1  =  1  TO  (Md-1  ) 

6380 

ENTER  @File;Xya<*  ) 

6385 

NEXT  I 

6390 

END  IF 

6395 

FOR  1=1  TO  Npairs 

6400 

ENTER  @File;Xya(* ) 

6405 

Ya=Xya< Itube*2 ) 

6410 

Xa=Xya( Itube*2+1 ) 

6415 

Yc=LOG(Xa ) 

6420 

Xc=LOG(Xa/Ya ) 

6425 

Sx=Sx+Xc 

,  I  tube 
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6430  Sy=Sy+Yc 

6435  Sx2=Sx2+Xc"2 

6440  Sxy=Sxy+Xc*Yc 

6445  IF  Iht-0  THEN 

6450  Xt=Ya 

6455  Ya=Ya/Xa 

6460  Xa=Xt 

6465  IF  Xa/l.E+6>Xul  THEN  Xu-l-Xa/1  .E+6 

6470  IF  Xa/l.E+B<Xll  THEN  Xll=Xa/l.E+6 

6475  END  IF 

6480  IF  Iht-0  THEN 

6485  X=<Xa*l.E-6-Xmin)*Sfx 

6490  Y-(Ya*l.E-3-Ynin)*Sfy 

6495  END  IF 

6500  IF  Iht-1  THEN 

6505  X=(Xa-Xmin)*Sfx 

6510  Y=(Ya*l .E-6-Ymin)*Sfy 

6515  END  IF 

6520  IF  Iht=2  THEN 

6525  X«<Xa/Ya-Xnin>*Sfx 

6530  Y«<Ya*l.E-3-Ymin)#Sfy 

6535  END  IF 

6540  IF  Y>100  OR  Y<0  THEN  6585 

6545  IF  5ym>3  THEN  PRINT  "SM" 

6550  IF  Syrc<4  THEN  PRINT  "SR  1.4,2.4" 

6555  PRINT  "PA",X,Y,"" 

6560  IF  Syn>3  THEN  PRINT  "SR  1.2,1.6" 

6565  IF  Sym=4  THEN  PRINT  "UC2 ,4 , 99 , 0 , -8 , -4 , 0, 0 , 8 , 4 ,0 ; " 

6570  IF  Sym=S  THEN  PRINT  "UC3 , 0 , 99 , -3 , -5 , -3 , 6 , 3 , 6 , 3 , -6 

6575  IF  Syn=6  THEN  PRINT  "UC0 , 5 . 3 , 99 , 3 , -8 , -6 , 0 , 3 , 8 ; " 

6580  IF  Sym=7  THEN  PRINT  "UC0 , -5 . 3 , 99 , -3 , 8 ,6 , 0 , -3 , -8 ; " 

6585       NEXT  I 

6590       BEEP 

6595       INPUT  "WANT  TO  LABEL  < 1  =  Y , 0=N  )? " , I lbl 

6600       IF  Ilbl=l  THEN 

6605  IF  Sym>3  THEN  PRINT  "SM" 

6610  IF  Syw<4  THEN  PRINT  "SR  1.4,2.4" 

6615  PRINT  "PA" ,Xal,Yal,  "" 

6620  IF  Sym>3  THEN  PRINT  "SR  1.2,1.6" 

6625  IF  Sym=4  THEN  PRINT  "UC2 , 4 , 99 , 0 , -8 , -4 , 0 , 0 , 8 , 4 ,0 ; " 

6630  .  IF  Syi*i=5  THEN  PRINT  "UC3  ,  0,  99  ,  -3  ,  -6  ,  -3  ,  6  ,  3  ,  6  ,  3  ,  -6  ;' 

6635  IF  Sym=6  THEN  PRINT 

6640  IF  Sym=7  THEN  PRINT 

6645  PRINT  "SM" 

6650  IF  Syei<4  THEN  PRINT 

6655  PRINT  "PR  2,-1.0;SR 

6660  Yal=Yal-5 

6665  BEEP 

6670  INPUT  "WANT  TO  ADD  ANOTHER  STRING  < 1=Y , 0=N >? " , las 

6675  IF  Ias=l  THEN 


"UC0, 

3.Z 

.,99,3 

-8 ,  -6 

0 

3 

,8; 

"UC0, 

-5 

3,99,- 

-3,8,6 

0 

~> 

- » 

"PR    2 

,0' 

1.0,1 

.8; 

LB";D. 

.file* 

"  ' 
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6680 
6685 
6690 
6695 
6700 
6705 
6710 
6715 
6720 
6725 
6730 
6735 
6740 
6745 
6750 
6755 
6760 
6765 
6770 
6775 
6780 
6785 
6790 
6795 
6800 
6805 
6810 
6815 
6820 
6825 
6830 
6835 
6840 
6845 
6850 
6855 
6860 
6865 
6870 
6875 
6880 
6885 
6890 
6895 
6900 
6905 
6910 
6915 
6920 
6925 


STRING", Label$ 
1.0,1.8;LB";Label$; 


BEEP 

INPUT  "ENTER  THE 
PRINT  "PR  2,0;5R 
GOTO  6665 
END  IF 
END  IF 
BEEP 

INPUT  "WANT  TO  COMBINE  ANOTHER  FILE?  < 1=Y , 0=N > " , Icomb 
ASSIGN  ©File  TO  * 
Xll=5 
Xul=4S 

IF  lcomb<>0  THEN  6220 
BEEP 

INPUT  "WANT  TO  PLOT  A  LEAST-SQUARES  LINE  <  1  =  Y ,0=N  ) " , I  Is 
IF  Ils=l  THEN 
BEEP 

INPUT  "SELECT  EXPONENT:  0=COMPUTE ,  1=0 .75  ) " , Iexp 
BEEP 

INPUT  "SELECT  CURUE  TYPE  < 0=SOLID , 1  =  DASHED ) " , 1 1 1 
Ilt-Ilt+1 
PRINT  "SM" 
IF  lexp=0  THEN 

Bb=(Nrun«Sxy-Sy»Sx  >/(  Nrun*Sx,2-SxA2  ) 
ELSE 

Bb=.75 
END  IF 

Aa=<Sy-Bb*Sx  )/Nrun 
Aa=EXP<Aa) 
PRINTER  IS  1 

PRINT  USING  "10X,""a  =  "",Z.4DE";Aa 
PRINT  USING  "10X,""n  =  "",Z.4DE";Bb 
PRINTER  IS  705 
ln=0 

IF  Ih+=0  THEN  Xxstep=Xstep/40 

IF  lht>0  THEN  Xxstep=Xstep/10 

FOR  Xa=Xll  TO  Xul  STEP  Xxstep 

IF  Xa>.99*Xmax  THEN  6995 

IF  Iht-1  THEN  Ya=Aa*XaABb 

IF  lht=0  THEN  Ya  =  Aa"< 1/Bb  )*(Xa*l .E+6 )'<  (Bb-1  )/Bb ) 
IF  Iht=2  THEN  Ya=Aa*XaA ( Bb-1 ) 
IF  lht=0  THEN 

Y=( Ya*l.E-3-Ymin)*Sfy 
X-(Xa-XMin)*Sfx 
END  IF 
IF  Iht=l  THEN 

Y»<Ya*l.E-B-Ynin)*Sfy 
X=(Xa-Xnin)*Sfx 
END  IF 

IF  Iht=2  THEN 
.  Y=(Ya*l .E-3-Ymin>*Sfy 
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6930 

X=(Xa-Xmin>»Sfx 

6935 

END  IF 

6940 

IF  Y<0  THEN  Y=0 

6945 

IF  Y>100  THEN  6990 

6950 

IF  Ilt=l  THEN 

6955 

PRINT  "PA",X,Y,"PD" 

6960 

ELSE 

6965 

In=In+l 

6970 

Ir-In  MOD  lit 

6975 

IF  Ir=l  THEN  PRINT  "PA" , X , Y , "PD" 

6980 

IF  lr=0  THEN  PRINT  "PA" , X , Y , "PU" 

6985 

END  IF 

6990 

NEXT  Xa 

6995 

PRINT  "PU" 

7000 

END  IF 

7005 

lcomb=0 

7010 

60T0  6185 

7015 

END  IF 

70Z0 

PRINT  "PU  SM" 

7025 

6EEP 

7030 

INPUT  "WANT  TO  PLOT  NUSSELT  LINE  ( 1=Y ,0=N  )? " , Inp 

7035 

IF  lnp=0  THEN  7125 

7040 

BEEP 

7045 

INPUT  "ENTER  TSAT  <DEFAULT=18  DE6  C)",Tsat 

7050 

Hfg=FNHfg(Tsat  > 

7055 

Xll=5 

7060 

Xul=45 

7065 

FOR  Xa=Xll  TO  Xul  STEP  Xstep/50 

7070 

Tfilm=Tsat-Xa*.5 

7075 

Kf-FNK(Tfiln) 

7080 

Rhof=FNRho(Tfilm) 

7085 

Muf=FNMu(Tfilrc) 

7090 

Ya=. 728* ( K fA 3*Rhof A2* 9.81 *Hf g/( Muf *Do*Xa  )  )". 

25 

7095 

X  =  (Xa-Xmm)»Sfx 

7100 

Y=(Ya*l.E-3-Ymin>*Sfy 

7105 

PRINT  "PA",X,Y,"PD" 

7110 

NEXT  Xa 

7115 

PRINT  "PU  PA  0,0" 

7120 

PRINT  "PU  PA  0,0  SP0" 

7125 

SUBEND 

7130 

SUB  Plo+3 

7135 

CON  /Drl/  Star, Syn, Icon 

7140 

COM  /Fid/  Ift 

7145 

DIM  0(9) 

7150 

Fw-1 

7155 

PRINTER  IS  1 

7160 

BEEP 

7165 

PRINT  USING  "4X, ""Select  Option  X-Y  Limits: 

7170 

PRINT  USING  "6X,""0  Use  default  values 

7175 

PRINT  USING  "6X,""1  Use  new  values 
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7180  INPUT  Okd 

7185  PRINTER  IS  705 

7190  IF  Okd=0  THEN 

7195       Xnin=0 

7200       Ymm=0 

7205       Xmax=15 

7210       Ymax=15 

7215       Xstep=3 

7220       Ystep=3 

7225  ELSE 

7230       BEEP 

7235       INPUT  "ENTER  MINIMUM  AND  MAXIMUM  X-UALUES" , Xmin , Xmax 

7240       BEEP 

7245       INPUT  "ENTER  MINIMUM  AND  MAXIMUM  Y-UALUES" , Ymin , Ymax 

7250       BEEP 

7255       INPUT  "ENTER  STEP  SIZE  FOR  X-AXIS" , Xstep 

7260       BEEP 

7265       INPUT  "ENTER  STEP  SIZE  FOR  Y-AX IS" , Ystep 

7270  END  IF 

7275  BEEP 

7280  PRINT  "IN;SP1;IP  2300,1800,8300,6800;" 

7285  PRINT  "SC  0 , 100 , 0 , 100 ; TL  2,0;" 

7290  Sfx=100/(Xmax-Xmin> 

7295  Sfy=100/( Ymax-Ymin) 

7300  BEEP 

7305  lcg=0 

7310  INPUT  "LIKE  TO  BY-PASS  CAGE  ( 1-Y , 0=N=DEFAULT )? " , leg 

7315  IF  Icg=l  THEN  7620 

7320  PRINT  "PU  0,0  PD" 

7325  FOR  Xa=Xmin  TO  Xmax  STEP  Xstep 

7330        X=(Xa-Xmin)*Sfx 

7335        PRINT  "PA";X,",0;  XT; " 

7340  NEXT  Xa 

7345  PRINT  "PA  100,0;PU; " 

7350  PRINT  "PU  PA  0,0  PD" 

7355  FOR  Ya=Ymin  TO  Ymax  STEP  Ystep 

7360        Y=(Ya-Ymin)*Sfy 

7365        PRINT  "PA  0, ";Y,"YT" 

7370  NEXT  Ya 

7375  PRINT  "PA  0,100  TL  0  2" 

7380  .   FOP  Xa=Xmin  TO  Xrcax  STEP  Xstep 

7385        X=<Xa-Xmin)*Sfx 

7390        PRINT  "PA" ;X," ,100;  XT" 

7395  NEXT  Xa 

7400  PRINT  "PA  100,100  PU  PA  100,0  PD" 

7405  FOR  Ya=Ymin  TO  Ymax  STEP  Ystep 

7410        Y=(Ya-Ynin)*Sfy 

7415        PRINT  "PD  PA  100,",Y,"YT" 

7420  NEXT  Ya 

7425  PRINT  "PA  100,100  PU" 
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7430  PRINT  "PA  0,-2  5R  1.5,2" 

7435  FOR  Xa-Xnin  TO  Xmax  STEP  Xstep 

7440        X-<Xa-Xmin)*Sfx 

7445        PRINT  "PA";X, " ,0; " 

7450        IF  Xa<l  AND  Xa<>0  THEN  PRINT  "CP  -1 . 5 , -1 ; LB0 ; PR  -l,0;LB";Xa; 

7455        IF  Xa=0  THEN  PRINT  "CP  -.5,-1 ;LB0" 

74B0        Xin=0 

74B5        IF  Xa  MOD  1=0  THEN  Xin=l 

7470        IF  Xa>=10  THEN  PRINT  "CP  -2 , -1 ; LB" ; Xa ; " " 

7475        IF  Xa>l  AND  Xa<10  AND  Xin=l  THEN  PRINT  "CP  -1 .25 , -1 ; LB" ; Xa ; ' 

7480        IF  Xa>l  AND  Xin=0  THEN  PRINT  "CP  -2 , -1 ; LB" ; Xa ; " " 

7485        IF  Xa=l  THEN  PRINT  "CP  -1,-1;LB1.0" 

7490  NEXT  Xa 

7495  PRINT  "PU  PA  0,0" 

7500  Iht=2     !  MODIFY 

7505  FOR  Ya-Ymin  TO  Ynax  STEP  Ystep 

7510        Y=(Ya-Ymin)*Sfy 

7515        PRINT  "PA  0, " ; Y , " " 

7520        IF  lht=0  AND  Ya>0  THEN  PRINT  "PR  2,0" 

7525        IF  Ya<:i  AND  Ya<>0  THEN  PRINT  "CP  -4  ,  -  .  25  ;  LB0  ;  PR  -2,0;LB";Ya 

7530        IF  Ya=0  THEN  PRINT  "CP  -2,-.25;LB0" 

7535        IF  Ya>l  AND  I h  t  <  2  THEN  PRINT  "CP  -5 , - . 25 ; LB" ; Ya ; "  " 

7540        IF  Ya>9  AND  Iht=2  THEN  PRINT  "CP  -4 , - . 25 ; LB " ; Ya ; " " 

7545        IF  Ya>0  AND  Ya<10  AND  Iht=2  THEN  PRINT  "CP  -3 , - . 25;LB" ; Ya; "' 

7550  !       IF  Ya=l  THEN  PRINT  "CP  -4  ,  - . 25 ; LB1 . 0" 

7555  NEXT  Ya 

75B0  IF  Okd=0  THEN 

7565       Xlabel$="X" 

7570       Ylabel$="Y" 

7575  ELSE 

7580       BEEP 

7585       INPUT  "ENTER  X-LABEL" ,Xlabel$ 

7590       BEEP 

7595       INPUT  "ENTER  Y-LABEL" , YlabelS 

7B00  END  IF 

7B05  PRINT  "SR  1.5,2  PU  PA  50,-10  CP" ; -LEN( X label$  )/2 ; "0 ; LB" ; Xlabel$; 

7610  PRINT  "PA  -11,50  CP  0, " ;-LEN( Ylabel$ >/2*5/6 ; "DI  0, 1 ; LB" ; YlabelS 

7615  PRINT  "CP  0,0  DI" 

7620  Nrun=0 

7625  BEEP 

7630  INPUT  "WANT  TO  PLOT  DATA  FROM  A  FILE  <  1=Y,0=»N )?"  ,Ok 

7635  Okp=0 

7640  IF  Ok=l  THEN 

7645       BEEP 

7650       INPUT  "ENTER  THE  NAME  OF  THE  PLOT  DATA  FILE" , D_f i le$ 

7655       ASSIGN  ©File  TO  D_file$ 

7660       IF  lcomb<>0  THEN  7690 

7665  S>:=0 

7670  Sy=0 

7675  Sx2=0 
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7680  Sxy=0 

7685  Md=l 

7690  BEEP 

7695  INPUT  "ENTEP  THE  BEGINNING  RUN  NUMBER  <DEF=l>",Md 

7700  Npairs=9 

7705  BEEP 

7710  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED  ( DEF=9  ) " .Npairs 

7715  BEEP 

7720  INPUT  "SELECT  TUBE  NUMBER  ( 0=TOP , l=SECOND )",Itube 

7725  Nrun=Nrun+Npairs 

7730  PRINTER  IS  1 

7735  BEEP 

7740  PRINT  USING  "4X ," "Select  a  symbol: 

7745  PRINT  USING  "4X,""  1   Star    2   Plus  sign 

7750  PRINT  USING  "4X,""  3   Circle   4   Square" "n 

7755  PRINT  USING  "4X,""  5   Rombus 

7760  PRINT  USING  "4X,""  6   Right-side-up  triangle 

7765  PRINT  USING  "4X,""  7   Up-side-down  triangle""" 

7770  INPUT  Syn 

7775  PRINTER  IS  705 

7780  IF  Syi*i=l  THEN  PRINT  "SM*" 

7785  IF  Syrc=2  THEN  PRINT  "SM+" 

7790  IF  Sym=3  THEN  PRINT  "SMo" 

7795  Md=Itube*9 

7800  IF  Md>l  THEN 

7805  FOR  1=0  TO  (Md-1  ) 

7810  ENTER  @File;X,Y 

7815  NEXT  I 

7820  END  IF 

7825  FOR  1=1  TO  Npairs 

7830  ENTER  ©File; Xa,Ya 

7635  Sx=Sx+Xa 

78A0  Sy=Sy+Ya 

7845  Sx2=Sx2  +  Xa'-2 

7850  Sxy=Sxy+Xa*Ya 

7855  X=(Xa-Xmin)*Sfx 

7860  Y=< Ya-Ymin)*Sfy 

7865  IF  Y>100  OR  Y<0  THEN  7910 

7870  IF  Sym>3  THEN  PRINT  "SM" 

7875  IF  Sym<4  THEN  PRINT  "SR  1.4,2.4" 

7880  PRINT  "PA",X,Y,"" 

7885  IF  Sym>3  THEN  PRINT  "SR  1.2,1.6" 

7890  IF  Sym=4  THEN  PRINT  "UC2 ,4 , 99 , 0 , -8 , -4 ,0 , 0 , 8 , 4 , 0 ; " 

7895  IF  Syn=5  THEN  PRINT  "UC3, 0 , 99, -3 , -6 , -3 , 6 , 3 , 6 , 3 , -6 ; " 

7900  IF  Sym-6  THEN  PRINT  "UC0,5. 3 , 99 , 3 , -8 , -6 , 0 , 3 , 8 ; " 

7905  IF  Sym=7  THEN  PRINT  "UC0 ,-5.3 , 99, -3,8 ,6,0 ,-3,-8; " 

7910  NEXT  I 

7915  BEEP 

7920  INPUT  "WANT  TO  LABEL  < 1  =  Y , 0=N  )? " , I lbl 

7925  IF  Ilbl-1  THEN 


210 


7930 
7935 
7940 
7945 
7950 
7955 
7960 
7965 
7970 
7975 
7980 
7985 
7990 
7995 
8000 
8005 
8010 
8015 
8020 
8025 
8030 
8035 
8040 
8045 
8050 
8055 
80B0 
80B5 
8070 
8075 
8080 
8085 
8090 
8095 
8100 
8105 
8110 
8115 
8120 
8125 
8130 
8135 
8140 
8145 
8150 
8155 
81B0 
8165 
8170 
8175 


IF  Sym=5  THEN  PRINT 
IF  Sym=6  THEN  PRINT 


IF  Syn>3  THEN  PRINT  "SM" 

IF  Sym<4  THEN  PRINT  "SR  1.4,2.4" 

PRINT  "PA",Xal,Yal,"" 

IF  Sym>3  THEN  PRINT  "SR  1.2,1.6" 

IF  Sym=4  THEN  PRINT  "UC2  ,  4  ,  99  ,  0  ,  -8  ,  -4  ,  0  ,  0  ,  8  ,  4  ,  0  ;  " 

"UC3, 0,99, -3, -6, -3, 6, 3, 6, 3,-6: " 
"UC0,5. 3, 99,3,-8, -6,0,3, 8;" 
IF  Sym=7  THEN  PRINT  "UC0 , -5 . 3 , 99, -3 , 8 , 6 , 0 , -3 , -8 ; " 
PRINT  "SM" 

IF  Sym<4  THEN  PRINT  "PR  2,0" 
PRINT  "PR  2,-1.0;SR  1 . 0 , 1 . 8 ; LB" ; D_f i le$ ; " " 
Yal=Yal-5 
BEEP 

INPUT  "WANT  TO  ADD  ANOTHER  STRING  < 1-Y.0-N >?" , las 
IF  Ias=l  THEN 
BEEP 

INPUT  "ENTER  THE  STRING" , Label? 
PRINT  "PR  2,0;SR  1 . 0 ,  1  . 8 ; LB" ; Label$ ; " " 
GOTO  7990 
END  IF 
END  IF 
BEEP 

INPUT  "WANT  TO  COMBINE  ANOTHER  FILE?  < 1  =  Y , 0=N  ) " , Icomb 
ASSIGN  @File  TO  * 
IF  IcornbOO  THEN  7B45 
Ils  =  l 
BEEP 

INPUT  "WANT  TO  PLOT  A  LEAST-SQUARES  LINE  < 1=DEF=YES ,0=NO  ) " , 1 1 s 
IF  Ils=l  THEN 
BEEP 
INPUT  "SELECT  CURUE  TYPE  <  0=SOLID  ,  1=DASHED  ) "  -,  1 1 1 

irt-nt+1 

PRINT    "SM" 

IF    lexp=0    THEN 

Bb=(  Nrun*Sxy-Sy*Sx  )/  (  Nrun*Sx2-Sx  "2  ) 
ELSE 


Bb=.75 

END  IF 

Aa=(Sy-Bb*Sx  )/Nrun 

PRINTER  IS  1 

PRINT  USING  "10X,""a  =  "",MZ.3DE' 

;Aa 

PRINT  USING  "10X,""n  =  "".MZ.3DE' 

';Bb 

PRINTER  IS  705 

ln=0 

FOR  Xa=Xmin  TO  Xmax  STEP  ( Xmax-Xmin ) 

Ya=Aa+Xa*Bb 

Y=(Ya-Ynin^*Sfy 

X=(Xa-Xmin)*Sfx 

IF  Y<0  THEN  Y=0 

IF  Y>100  THEN  GOTO  8220 
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8180  IF  Ilt=l  THEN 

8185  PRINT  "PA",X,Y, "PD" 

8190  ELSE 

8195  In=In+l 

8200  Ir-In  MOD  lit 

8205  IF  Ir-1  THEN  PRINT  "PA" , X , Y , "PD1 

8210  IF  lr=0  THEN  PRINT  "PA" , X , Y , "PU' 

8215  END  IF 

8220  NEXT  Xa 

8225  PRINT  "PU" 

8230       END  IF 

8235       lcomb=0 

8240       GOTO  7620 

8245  END  IF 

8250  PRINT  "PU  PA  0,0" 

8255  PRINT  "PU  PA  0,0  5P0" 

8260  SUBEND 

8270  DEF  FNFncaKI.T.Pc  ) 

8280  IF  1=0  THEN 

8290       Mdot0=.00206*Pc-. 00414 

8300       Mdot24=.00225*Pc-. 00329 

8320  END  IF 

8330  IF  1=1  THEN 

8340       Mdot0=.00204*Pc-. 00615 

8350       Mdot24=.00216*Pc-. 00138 

8370  END  IF 

8380  IF  1=2  THEN 

8390       Mdot0=.00208*Pc-.0036 

8400       Mdot24=.00223*Pc-. 00301 

8420  END  IF 

8430  IF  1=3  THEN 

8440       Mdot0=.0028*Pc-. 00122 

8450       Mdot24=.0029»Pc+.0048 

8460  END  IF 

8470  Mdt=<T*<Mdot24-Mdot0)/24 )+Mdot0 

8480  RETURN   Mdt 

8490  FNEND 
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